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ing problems and the chance for greater responsibility 
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You Can't Inspect Quality into a Piece! 


You must put quality into a piece a// along the way, in each opera- 
tion. That's the kind of work you get from the Diversey team... 
extreme quality and precision in every step. 

The missile component you see above was an intricate and diffi- 
cult job. You see it now passing final rigid inspection with flying 
colors because it was handled completely by Diversey. Nowhere 
else can you get this devotion and follow through on your com- MISSILE 
plete missile machining jobs. Diversey has the best men, the best SEND sosemaiinaie 
equipment, and the best team spirit to give you the finest missile FOR 
metal machining quality in the industry. FREE 


When your job must pass rigid inspection, have it done at 
Diversey. BOOKLET 
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LEADERS IN CONTOUR MACHINING: 


ENGINEERING COMPANY 
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FROM NOSE TO NOZZLE, FROM FIN TO FIN, CONTOUR TURNED PARTS—WITH PRECISION BUILT IN 


Juty 1958 437 


x 
| 
: 
Diversey 
7 
f 
f 
| 


74 publication of the 


“—-AMERICAN ROCKET so ¥: 
Research and Development 

IRWIN HERSEY—DIRECTOR OF PUBLICATIONS 

SUMMERFIELD 


EDITOR 
ASSISTANT EDITOR LARKIN JOYNER 
ART EDITOR JOHN CULIN 


ASSOCIATE EDITORS 


AL! BULENT CAMBEL, Northwestern University 
IRVIN GLASSMAN, Princeton University 
M. H. SMITH, Princeton University 


CONTRIBUTORS 


MARSHAL FISHER, Princeton University 
GEORGE F. McLAUGHLIN 


WALTER BRUNKE 


ADVERTISING & PROMOTION MANAGER 
WILLIAM CHENOWETH 


D. C. Emery & Associates James C. Galloway & Co. 
155 East 42 St., New York, N. Y. 6535 Wilshire Blvd., Los Angeles, Calif. 
Telephone: Yukon 6-6855 Telephone: Olive 3-3223 


Jim Summers & Associates R. F. and Larry Pickrell 
35 E. Wacker Dr., Chicago, Ill. 318 Stephenson Bldg., Detroit, Mich. 
Telephone: Andover 3-1154 Telephone: Trinity 1-0790 


Louis J. Bresnick John W. Foster 
304 Washington Ave., Chelsea 50, Mass. 239 Ath Ave., Pittsburgh, Pa. 
Telephone: Chelsea 3-3335 Telephone: Atlantic 1-2977 


AMERICAN ROCKET SOCIETY 


Founded 1930 
OFFICERS 

President George P. Sutton 
Vice-President - John P. Stapp 
Executive Secretary -— James J. Harford 
Secretary | A. C. Slade 
Treasurer Robert M, Lawrence 
General Counsel - Andrew G. Haley 


BOARD OF DIRECTORS 
Terms expiring on dates indicated 
Krafft Ehricke, 1959 H. S. Seifert, 1958 
S. K. Hoffman, 1958 a ~ a K. R. Stehling, 1958 


Simon Ramo, 1960 ee a rte Martin Summerfield, 1959 
H. W. Ritchey, 1959 Wernher von Braun, 1960 


Maurice J. Zucrow, 1960 


TECHNICAL DIVISION CHAIRMEN 


David G. Simons, Human Factors John F. Tormey, Propellants and 
Lawrence S. Brown, Instrumenta- Combustion 

tion and Guidance Brooks T. Morris, Ramjet 
Y. C. Lee, Liquid Rocket Ivan Tuhy, Solid Rocket 


Krafft A. Ehricke, Space Flight 


438 


Scope of JET PROPULSION 


This Journal is a publication of the Ameri- 
can Rocket Society devoted to the advance- 
ment of the field of jet propulsion through 
the dissemination of original 
new knowledge or new deve opments. As 
used herein, the term “jet propulsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-consuming 
engines "and underwater systems as well as 
rockets. JET PROPULSION is open to con- 
tributions dealing not only with propulsion 
but with other aspects of jet-propelled flight, 
such as flight mechanics, guidance, telemeter- 
ing, and research instrumentation. In- 
creasing emphasis will be given to the scientific 
problems of extraterrestrial flight. 
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Manuscripts must be as brief as the proj er 
presentation of the ideas will allow. Ex- 
clusion of dispensable material and concise- 
ness of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-space: typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract, and 12 illustrations’ 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments. 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are usually published within two 
months of the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication according to the above-stuated 
requirements as to subject scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as & 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. [If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to begiven 
as follows: for journal articles: authors first, 
then title, journal, volume, year, page numbers; 
for books: authors first, then title, publisher, 
city, edition, and page or chapter numbers. 
Line drawings must be clearand sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accomplished by 
written assurance as to security. clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Engineering, Prince- 
ton University, Princeton, N. J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 
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Recent Advances in Convective Heat Transfer With 
Dissociation and Atom Recombination 
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The author is a 1957-1958 General Electric Charles A. Coffin 
Fellow, working toward the Ph.D. degree in aeronautical engi- 
neering at Princeton University. A Daniel and Florence 
Guggenheim Jet Propulsion Fellow during 1955-1957, Mr. Rosner 
received his BME summa cum laude at the City College School 
of Technology, New York, N. Y., in June 1955. He is now en- 
gaged in research on the problem of heat transfer in dissociated 
gases. 


I. Introduction 

LTHOUGH energy transfer processes in partially disso- 
{“& ciated gases had been investigated as early as 1904 by 
Nernst (36),? it is only in the last few years that the aerody- 
namics of propulsion has attained a state of development 
which obliges the engineer to consider the effects due to dis- 
sociation and recombination reactions in the fluid surrounding 
a vehicle (the external problem) and in the working fluid used 
in a propulsion or cooling system (the internal problem). 
These are effects by which the molecular nature of “real’’ 
fluid media manifests itself, particularly at high temperature. 
In particular, it is the purpose of this paper to briefly describe 
recent progress in our understanding of heat transfer phe- 
nomena connected with dissociation and recombination proc- 
esses, with the discussion being limited to continuum flows in 
the absence of ionization, surface meiting, vaporization or 
sublimation, or subsurface reactions. Because of the rate at 
which the literature in this field is expanding, no claim can be 
made as to the completeness of either this review or the ac- 
companying bibliography. It is hoped, however, that at least 
an over-all picture is conveyed, and that the bibliography will 
be useful as an introduction to this class of problems for the 
interested reader. 

If one considers high Mach number flight, or entry into 
planetary atmospheres in which there are polyatomic con- 
stituents, then local deceleration of the medium relative to 
the body gives rise to temperatures which are high enough to 
cause fragmentation of even the most stable molecules. This 
may occur either behind the bow shock which stands off 
blunt-nosed bodies or in the shear layer near the lateral sur- 
faces of vehicles as a result of locally intense viscous heating. 
In either event, energy transfer to the surface of a vehicle de- 
pends critically upon whether the atoms or molecular frag- 
ments (“radicals’’) recombine exothermically in its immediate 
vicinity, at the surface of the body itself, or in the wake left 
by the vehicle. Thus, the net energy transfer depends upon 
flight conditions, surface conditions and physical scale in a 
way which could not be anticipated without making a de- 
tailed study of energy transfer processes in nonuniform, react- 
ing gaseous mixtures. As to the internal problem, there are 
now rocket propellant combinations characterized by high 
combustion temperatures and substantial atom and radical 
concentrations in the products. Again, chamber and nozzle 
cooling requirements cannot be understood (50) without an 
investigation of the effects of chemical rate processes. Of 
relevance also is the fact that gases which readily dissociate, 
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such as the halogens, may be competitive with more conven- 
tional coolants for certain applications. 

With these examples in mind, we proceed to discuss recent 
contributions to our present understanding of heat transfer 
phenomena with dissociation and recombination. 


Il. Effects of Gas-Phase and Surface Reaction © 


The effects to be discussed are those characteristic of the 
more general problem of nonequilibrium flows of reactive 
gaseous mixtures. This general subject is known by the rather 
ominous name: aerothermochemistry. Since dissociation and 
recombination reactions form only a small subclass of the 
class of all “reactions,’’ much of the discussion may be taken 
to include not only dissociation and recombination reactions 
but also excitation, de-excitation (32) and re-arrangement 
“reactions.”’ “Composition”? would then be considered in the 
broad sense of Ubbelohde (49), who first regarded each 
molecular quantum state as a distinct chemical specie for 
purposes of energy transfer calculations in pure polyatomic 
gases. In fact, if this point of view is adopted, the relation 
between the Eucken conductivity of pure polyatomic gases 
and the “effective” heat conductivity (Part IV) of reacting 
or nonuniform gas mixtures will be seen to be a close one. 

Several surveys of the fundamental equations of aerother- 
mochemistry are available (21, 25, 37), with the emphasis 
being, therefore, upon rate processes within the solution do- 
main, i.e., the gas phase. However, boundary conditions (44) 
for these equations must, in general, include the rates of chemi- 
cal processes at surfaces, i.e., heterogeneous combustion or 
catalysis rates. Heat transfer problems with dissociation and 
subsequent recombination owe most of their complexity and 
interest to the fact that both homogeneous and heterogeneous 
rate processes play an important role in determining the 
character of solutions. For example, it is possible for the 
chemical composition of a gas mixture near a surface to vary 
from point to point in the absence of gas-phase chemical rate 
processes (i.e., in chemically “frozen” flow). This will be the 
case when the bounding surface acts as a sink for several 
species present in the “free stream.”’ The diffusion processes — 
which accompany such nonuniformities of chemical composi- | 
tion are of fundamental importance in aerothermochemical 
problems regardless of their cause; i.e., they may be due to — 
heterogeneous and/or homogeneous chemical rate processes. 
In fact, this steady-state diffusion constitutes the energy trans-_ 
fer mechanism which sets this class of heat transfer problems — 
apart from those involving flow fields of constant composition. 
This is a result of the fact that, in general, each chemical — 
specie has a distinguished energy content relative to the ele-— 
ments in some arbitrarily chosen state. Thus, selective local — 
diffusion is itself a mechanism for energy transport which — 
must be dealt with if subsequent chemical change occurs in | 
the gas phase or at the bounding surface. We say, in a loose. 
sense, that “chemical” as well as “thermodynamic” energy — 
‘an be transported through such gas mixtures, as discussed in — 
Part IV. In addition, exothermic or endothermic chemical — 
reactions in the gas phase strongly modify the steady-state — 
temperature distributions (see Part V) and thus directly in- 
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fluence the magnitude of ordinary conductive energy transfer, 
as well as the steady-state distributions of all fluid properties 
which are temperature dependent. 


III. Macroscopic Approach to 
Aerothermochemistry 


It has not been fruitful, as yet, to adopt a completely micro- 
scopic approach to problems in aerothermochemistry. This 
‘an be attributed to the present state of development of the 
kinetic theory of dilute gas mixtures in which chemical reac- 
tions are included, and the fact that the cross sections for 
elementary reactions have not been accurately predicted by 
quantum-mechanical methods. Instead gaseous mixtures are 
regarded as continua in which there is normal transfer of 
momentum (10) between particles of the different components. 
The conservation equations of the mechanics of continuous 
fluid media are written; these are the well-known equations of 
continuity (for the individual species), momentum (for the 
mixture) and energy (for the mixture), modified so as to in- 
clude the effects of chemical reaction explicitly in the con- 
tinuity equations, and implicitly in the momentum and energy 
equations. It is significant, and perhaps unsatisfactory from 
the point of view of consistency, that a curious mixture of 
disciplines must be drawn upon to make this system of partial 
differential equations a determined one. The “flux vectors’ 
which appear in the equations must be written in terms of a 
set of phenomenological coefficients, the macroscopic proper- 
ties of the mixture, and their spatial gradients. For this pur- 
pose either recourse is made to the linear theory of irreversible 
thermodynamics, or one descends to the microscopic level, 
i.e., the kinetic theory of dilute monatomic znert gas mixtures; 
the latter provides not only the form of the flux vectors for 
inert gas mixtures but also expressions for the necessary phe- 
nomenological coefficients in terms of ‘fundamental’ molecu- 
lar properties. In either event, only linear transport terms 
are considered. However, as a result of the large departures 
from equilibrium often caused by chemical reactions, it is 
found (24) that the chemistry does not strictly fit into a linear 
transport theory. For this reason recourse is made to classical 
chemical kinetics to provide independently phenomenological 
relations for the rates of elementary reactions in terms of 
both the macroscopic variables themselves, and a number of 
new chemical parameters. The assumption is made that 
these kinetic laws, determined largely under uniform condi- 
tions, hold locally in the nonuniform velocity, temperature 
and composition fields to be determined. The resulting set of 
equations, which then includes parameters associated with 
the thermodynamic behavior of the gas mixture, “‘transport”’ 
properties and chemical kinetics, forms the basis of our present 
understanding of reacting continua, and, in particular, heat 
transfer in dissociating gas mixtures. 

Confronted with this formidable array of equations, bound- 
ary conditions and possible “special cases,” one can pose an 
enormous number of problems, some more tractable than 
others, some more realistic or technically important than 
others. On the one hand, one can construct simple mathe- 
matical models and seek closed-form results for a class of 
flows applicable to a “general” gas. Such models are intended 
to illustrate the dependence of the solutions upon the parame- 
ters pertinent to the problem and each study attempts to in- 
clude at least one of the couplings which characterizes prob- 
lems in aerothermochemistry. On the other hand, with the 
advent of machine computation, one can apply comparatively 
brute force techniques to rather comprehensive problems in- 
volving specific gas mixtures and later sort out the parameter 
dependences for cases of current technical interest. Both ap- 
proaches contribute greatly to our understanding of this 
general class of phenomena, and both appear in the literature. 
Simplifications can and have been introduced in many direc- 
tions, through treatment of steady-state problems charac- 
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1. Reduced spatial dimensionality. 
2. Simple hydrodynamic flow configurations; e.g., Couette 
flow, two-dimensional and axi-symmetric boundary layer 


flows, ete. 


3. Consideration of simple gas models; e.g., incompressi- 


ble “frozen” fluids, dilute “binary” gas mixtures, etc. 


4. Extremes of gas phase reaction rates; e.g., chemically 
frozen, quasi-equilibrium flows; or the adoption of a simple 
phenomenological rate law. 

5. Extremes of wall catalytic activity; e.g., noncatalytic, 
“fully” catalytic; or the adoption of a simple phenomeno- 
logical heterogeneous rate law. 

Because attention is here focused principally on heat trans- 
fer problems of engineering interest, emphasis will be upon 
“boundary layer” solutions to the macroscopic conservation 
equations, and their experimental verification. 


IV. The Energy Flux Vector—** Effective” 
Thermal Conductivity 


The simplest “flow field” in which to study energy trans- 
port in dissociating mixtures is that for which the mean mass 
motion vanishes everywhere, i.e., the ideal ‘conductivity 
cell.’’. Such studies and their local counterparts have been 
carried out in the past half-century (36, 11, 22, 23, 33, 39, 7) 
in an effort to predict the large “effective’’ thermal conduc- 
tivity of reacting gaseous mixtures as well as the observed 
variation of this conductivity with temperature and, particu- 
larly, pressure. 

The flux vector Jg which appears in the energy equation 
represents the directed local energy transport per unit time 
and area relative to the center of gravity of a moving fluid 
element. For the special case of the conductivity cell, J is 
simply constant across the cell; it is the eigenvalue corre- 
sponding to a steady-state solution of the energy equation for 
prescribed wall conditions. Locally, it can be anticipated that 
this energy flux will be the sum of contributions from: 

(a) Molecular transfer of translational, rotational and 
vibrational energy. 

(b) Transfer of chemical energy (i.e., energy relative to the 
elements in their standard states) owing to concentration 
(Fick) diffusion of distinct chemical species. 

(c) Thermo-diffusion; i.e., specie and, hence, energy trans- 
port due to the presence of a temperature gradient (Soret 
effect). 

(d) Diffusion thermo-effect; molecular transport of trans- 
lational, rotational and vibrational energy as a result of the 
presence of concentration gradients (Dufour effect). 

(e) Pressure diffusion; specie and, hence, energy trans- 
port due to the presence of a pressure gradient. This will be 
absent in the “conductivity cell.” 

(f) Thermal or chemi-luminescent radiant energy flux. 

Ordinarily, mechanisms (c) to (f) are neglected as being of 
secondary importance, although the effects of thermal diffu- 
sion have been studied (33, 39). It has been shown that if the 
gas phase reaction rates are fast (see Part E) in either the for- 
ward or reverse direction, the assumption of local thermo- 
chemical equilibrium is a good one, i.e., the actual com- 
position-temperature profiles differ only slightly from those 
corresponding to local thermochemical equilibrium, except in 
the immediate vicinity of walls. In this case, and only in 
this case, the reacting gas behaves as though it were a pure 
substance with an enhanced “effective” thermal conductivity 
Aert, in the sense that locally the energy flux vector may be 
written: Jo = —AetVT. This effective conductivity is de- 
composable into a “frozen” and “reactive” contribution, i.e., 
ers = Ay + A,; the “frozen” conductivity being the conduc- 
tivity one would compute on the basis of Chapman-Enskog 
theory (21) for the local mixture, including the appropriate 
Eucken corrections; the ‘reactive’ contribution being that 
due to local transport of chemical energy through the mix- 
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ture. The effective heat conductivity is thus seen to be the 
Eucken conductivity for polyatomic molecules generalized so 
as to include true mixtures. 

Several investigators have given equivalent expressions for 
the reactive contribution A, in terms of the thermodynamic 
coordinates of the mixture, enthalpy changes across the par- 
ticipating chemical reactions, and the binary and self diffusion 
coefficients. For the simple dissociation: A, @ 2A character- 
ized by the enthalpy change AH these results reduce to 


where mechanisms (c) to (f) have been neglected, a is the de- 
gree of dissociation (weight fraction of atoms), and Dj» is the 
atom-molecule binary diffusion coefficient. The relative im- 
portance of transport of chemical energy in nonuniform mix- 
tures is reflected in the value of the ratio ),/A,, which may be 
computed using Equation [1] or the equivalent result 


A; Cer | 

Where Le, is the “frozen” Lewis number (pD.».Cp,)/d,,Cp, 
the “frozen” specific heat, and Cp, the equilibrium specific 
heat of the binary mixture. This ratio is found to be of the 
order of 10 for the dissociation of oxygen, hydrogen and the 
halogens at temperatures eorresponding to a@ = 0.50 (7). 
These treatments have been verified through studies of the 
observed pressure and temperature dependence of the con- 
ductivity of hydrogen fluoride vapor (7) and the N.O;-NO, 
system (9). 

In case local thermochemical equilibrium is a poor assump- 
tion, the mixture cannot be characterized by an “effective’’ 
thermal conductivity and each of contributions (b) to (f) must 
be treated individually. There is, however, a special case 
worth noting which is of practical importance. For binary 
mixtures in which only mechanisms 1 and 2 are considered, 
the nonequilibrium energy flux vector takes a particularly 
simple form provided the frozen Lewis number for the mixture 
is exactly unity. This is seen from the approximate expres- 
sion for the energy flux*® vector 


where h; and he are the individual enthalpies,* @ and (1 — a) 
are the respective weight fractions, h is the specific enthalpy of 
the mixture, and Cp, is again the frozen specific heat of the 
mixture. For atom-molecule binary mixtures a@ represents 
the “degree of dissociation” and D,. becomes the atom-mole- 
cule binary diffusion coefficient. Unfortunately, however, 
the frozen Lewis number for atom-molecule mixtures is not, 
in general, unity (28) so that solutions of the energy equation 
based upon this assumption cannot be expected to reproduce 
all of the features of more accurate solutions for the case of 
“pure” dissociating gases. Solutions to the stagnation point 
problem for frozen Lewis numbers different from unity have 
recently been obtained, and are briefly discussed in Part V. 

The local equilibrium results for the heat conductivity bear 
directly upon many heat transfer problems of engineering in- 
terest. Recent work of Schotte (48) and Brokaw (6) suggests 
that, provided the appropriate mean values of the true 
molecular transport properties are introduced, the well-known 
empirical relations which correlate heat transfer measurements 
for the fully developed, laminar and turbulent internal flows of 
nonreacting fluids may also be used to compute heat transfer 

5’ The assumption Le; = 1 simplifies the energy flux vector to 
—[(A;/Ce,;) grad h] for multicomponent reacting mixtures as well, 
provided the gas may be replaced by an “effective’’ binary mixture 
with regard to diffusion processes. The implication and utility of 
these assumptions for treating chemically reacting laminar flows 
have been thoroughly discussed by Lees (26). 

‘Unless otherwise specified all enthalpies should be taken to 
include the enthalpy of formation at the reference temperature. 


rates in reversibly dissociating mixtures. The choice of mean 
values of the transport properties is expected to be more 
critical than in the corresponding nonreactive case since a 
prescribed wall to fluid (bulk) temperature difference is now 
associated with much larger changes in these properties. On 
the other hand, if the rates of diffusion and chemical processes 
are such that local thermochemical equilibrium cannot be 
maintained everywhere in the boundary layer, and, for ex- 
ample, atom recombination occurs at a cooled wall, then the 
“effective” conductivity loses its meaning on a molecular 
level, and the heat transfer depends both upon gas phase re- 
combination rates and the efficiency with which the wall re- 
combines incident atoms. In many such cases, the concen- 
tration, and, hence, enthalpy of the mixture at the wall, will 
not be known a priori (1, 8, 45) and, unlike those cases which 
may be treated with the conventional heat transfer correla- 
tions, the heat transfer rates become specific to the atom,‘sur- 
face combination under investigation. This is briefly dis- 
cussed in Part V-2. 


V. Chemical Reactions With Convection 


Parallel to the recent “conductivity cell” and local energy 
flux vector studies several solutions to problems involving 
chemical reactions in shear flows have recently been given. 
The effects of reactions on such flows and the important limit- 
ing cases may best be discussed in terms of dimensionless 
numbers which involve the “chemistry” of the process. For 
this purpose, it is convenient to think in terms of the charac- 
teristic times: Throws Teoll) T’chem chem. The first of 
these gives the time scale which characterizes any forced 
convection problem, viz., the time it would take a fluid a 
of the “undisturbed”’ stream to traverse the characteristic 
length of the problem For developing flows, this is always | 
of the same order as the time Tair: it takes a labile specie to. 
diffuse across the concentration boundary layer. The second — 
of these, introduced by the molecular nature of the medium, is — 
the average time, 7.11, between molecular encounters. The — 
presence of gas phase and surface reactions further introduces — 
the time scales 7%chem and 7" chem, Which may be thought of in | 
terms of the average time between those encounters which — 
lead to chemical reaction. Important regimes are then de-_ 
termined by the magnitude of the ratios of these times; e.g., — 
continuum flows are those for which >> ‘fast’ 
actions are those for which Tenem ~ Teo; Chemically “frozen” 
flows are those for which 7% >> Triow, local thermochemical 
equilibrium flows are those for which Triow >> chem, ete. 
Furthermore, these ratios and their energy counterparts de-_ 
termine the conditions under which, in particular, wrote il 
flows are aerothermochemically similar. The discussion of | 
convective heat transfer with dissociation and recombination | 
may now be cast in terms of these time scales. 


1. Laminar Flow 


Insight into some of the effects of surface reactions in 
laminar shear flows has been obtained from studies of the — 
properties of solutions to the boundary layer equations 
in cases for which the mass transfer and/or energy 
release does not affect the hydrodynamic field; constant 
density and transport properties usually being assumed. The 
limiting case of no volume reaction but pure surface catalysis | 
has been investigated by Chambré and Acrivos (1, 8) and this — 
writer (45), the object being to predict the streamwise dis- 
tribution of surface reaction rate in terms of the surface shear 
stress distribution and surface catalytic condition. The sig- 
nificant catalytic parameter has been displayed (1, 8, 45) and 
interpreted as a ratio of characteristic times (45). This — 


“catalytic parameter” thus becomes an additional sealing or _ 


similarity criterion for flows in which surface reactions (e.g., 
atom recombination) are of importance. Studies which in- 
clude the strong exothermic effects of gas-phase atom re- 
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combination (or, the strong endothermic effects in the case of 
dissociation) on the temperature profiles in compressible 
laminar boundary layers have also appeared. The relevant 
external heat transfer problems of greatest interest are those 
for which either the “free stream” seen by the boundary layer 
is partially dissociated, or dissociation occurs within the 
boundary layer itself by virtue of intense viscous heating. 
In the atmosphere of the earth both of these cases are aerody- 
namically encountered at hypersonic Mach numbers since 
thermal dissociation of stable diatomic molecules such as 
oxygen and nitrogen only becomes important for tempera- 
tures in excess of 2500 K. If the vehicle walls are necessarily 
maintained at lower temperatures then, aside from radiation, 
energy is transferred to the vehicle both by ordinary molecu- 
lar conduction and by exothermic recombination of atoms on 
the surface. The object of these analyses has been to predict 
the contribution of each mechanism to the local heat transfer 
rate for prescribed wall temperatures and free stream condi- 
tions. 

For practical as well as mathematical reasons much at- 
tention has recently been directed toward the blunt-nosed 
body. The large aerodynamic heating rate and low material 
heat capacity in the tip region of the sharp-nosed body pre- 
cludes its use at hypersonic speeds. As a result, there has 
been a renaissance of interest in the blunt-nosed body and, 
hence, interest in the local heating rate at the forward stag- 
nation point. A detached, locally normal shock is presumed 
to produce the dissociated “free stream” seen by the boundary 
layer, the assumption being that the division of the flow field 
behind the shock into a viscid and more-or-less inviscid region 
is a valid one. Mathematically it can be shown that the 
boundary layer conservation equations can be reduced to a 
system of exact ordinary differential equations for arbitrary 
gas phase reaction rates only near the stagnation point (11). 
Furthermore, an approximate similar solution with nonzero 
atom fraction at the wall can be found (16) for arbitrary first 
order wall kinetics in the case of stagnation point flow with no 
gas phase recombination. Recent analytical work on the 
stagnation point heat transfer problem with atom recombina- 
tion is classified below according to the comparative magni- 
tude of gas phase and wall reaction rates considered by each 
investigator. Wall rates are specified in terms of the “cata- 
lytic constant” ky defined by the first-order phenomenological 
rate law 


giving the mass rate R4 w of atom consumption in terms of the 
partial density, pC4, of atoms at the wall. A rate law of this 
form will adequately represent the recombination process in 
cases for which the equilibrium atom concentration cor- 
responding to the wall temperature is negligible compared to 
the actual local atom concentration (44,46) i.e., if the rate of 
dissociation of incident molecules is negligible. See Table 1.5 

The local equilibrium solutions numerically obtained by 
Fay and Riddell are based on the thermodynamic properties 
of argon-free air. On the other hand, solutions for the case of 
arbitrary gas phase recombination rates pertain to “air-like’”’ 
mixtures, i.e., those mixtures which may be treated as binary 
mixtures of “atoms” and “molecules” both as to transport 
properties and chemical kinetic behavior. Their numerical 
results for both the “frozen” and “equilibrium” cases have 
been expressed as simple engineering relations of the Nusselt 
form 


Nuz Delle CA,eE 
—— = {| (Pr,) <1 — —1) —— >}... [5] 
V Re, (= 

where the constants C, a, n and 8 are given for each case, and 
Ey is the molecular dissociation energy, c4,. and h? being, 
respectively, the weight fraction of atoms and stagnation 


5 Numbers in Table represent references appropriate to row 
and column headings. 
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Table 1 Laminar stagnation point heat transfer with dissociated 
free streem—a guide to recent analytical studies 


raitt/Tchem ———> 
\ sas phase} 

\ rates! infinite 
wall zero finite (quasi- 
rates (‘‘frozen’’) equil. ) 

_Tdiff 

zero 12,47 
(noncatalytic) 

finite 16, 46 

(catalytic) 12, 28 
“infinite’’ 12, 28 12,47 
(catalytic) | 16, 46 


enthalpy of the mixture at the outer edge of the boundary 
layer. As regards experimental verification, for realisti: 
choices of transport properties, the equilibrium relations of 
Lees (28) and Fay and Riddell (12) have successfully corre- 
lated recent shock-tube stagnation point heat transfer meas- 
urements in high temperature air over a wide range of simu- 
lated flight conditions (40, 42). 

The “fully” catalytic wall considered both by Lees and 
Fay and Riddell corresponds to the ability of the wall to 
maintain zero atom fraction at the surface but yet recombine 
atoms at a finite rate, i.e., kw must be infinite. That this is 
incompatible with the microscopic requirements of collision 
theory was independently recognized by Scala (46) and 
Goulard (16) each of whom considered the case of arbitrary 
wall catalytic activity. Ifa “fully” catalytic wall is taken to 
be one which recombines every incident atom then kw cannot 
be infinite (16), and the local atom fraction at the wall will be 
small but never zero. For the case of arbitrary catalytic ac- 
tivity the wall atom fraction is an essential unknown (1, 8, 45), 
i.e., its knowledge is equivalent, through (4), to a knowledge 
of the contribution of surface recombination to the heat 
transfer rate. 

The stagnation point problem for intermediate values of 
Tilow/T chem (Damkohler’s first ratio) has been treated (11) 
by introducing a phenomological reaction rate law appropriate 
to the three-body-collision recombination mechanism; again, 
with the gas treated as a binary mixture of “atoms” and 
“molecules.”” These solutions display the transition from 
the fully equilibrium to the chemically “frozen” layer, and 
the way in which the conductive and surface recombination 
contributions to the net heat transfer may be expected to vary 
with flight conditions and physical scale. 

The above treatments of the laminar stagnation point 
problem have lead to the following general conclusions: 

1. For “binary” gaseous mixtures, if recombination does 
take place, either entirely within the gas phase at local equi- 
librium or entirely at the wall, the net heat transfer rate, 
which depends principally upon the enthalpy difference h,? — 
h,, is nearly the same; the ratio of the two rates being of the 
order of (Le,;)w’/* (28). 

2. However, if the wall is noncatalytic, and, in addition, 
flight conditions are such that no appreciable recombination 
takes place in the gas phase, then the net heat flux can be less 
than half of that of case 1 (12, 28). 

This last conclusion, which reveals the sensitivity of the 
heat transfer rate to the nature of the vehicle surface, is con- 
sidered to be of potential importance in the solution of the 
“entry” problem. 

While the usual boundary layer similarity treatments, 
which lead to an exact system of ordinary differential equa- 
tions, do not strictly apply away from the vicinity of the 
stagnation point, these solutions have been extended down- 
stream along the body surface by means of the assumption of 
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“local similarity” (28). This assumption enables the heat 
transfer distribution to be determined from a knowledge of 
only the stagnation point heat transfer rate and the pressure 
and surface enthalpy distributions, and is expected to be good 
in case the variation of wall enthalpy is everywhere small com- 
pared to the enthalpy difference ho — hw across the layer, 
i.e., for “highly cooled” layers. The result has been verified 
over the forward portion of hemisphere cylinders by Hartwig 
(19) on the basis of steady-state heat transfer measurements. 
For a detailed review of frozen and equilibrium boundary 
lavers including a discussion of other approximate methods, 
the reader is referred to the forthcoming book of Hayes and 
Probstein (17). 

[n contrast to the laminar boundary layer (with dissociated 
free stream) which develops along blunt-nosed bodies, and in 
which viscous dissipation may be shown to be of secondary 
iniportance, there is an important class of hypersonic com- 
pressible shear flows in which viscous dissipation is the 
primary cause of dissociation; the free stream itself being 
undissociated. 

In order to gain insight into the relevant boundary layer prob- 
lem the simplest laminar shear flow with chemical reaction, 
i., Couette flow, has been considered. In particular, Liep- 
mann and Bleviss (29) have directed attention at the effect: of 
“ree-stream Mach number’’ (i.e., the Mach number corre- 
sponding to the linear velocity of the moving plate) on the en- 
thalpy and velocity profiles, recovery temperature, skin friction 
and heat transfer, assuming both local thermochemical equilib- 
rium and constant equilibrium Prandtl number. On the other 
hand, Broadwell (5) drops the assumption of local thermochem- 
ical equilibrium but introduces the simplifications of constant 
fluid properties, extremes of wall catalytic activity, and a 
“linearized”’ gas phase chemical source term, in order to in- 
vestigate the character of solutions for varous values of a reac- 
tion rate parameter which is essentially Taitt/T%chem- 

As in the case of Couette flow, with increases of upstream 
Mach number, the temperature rise due to viscous heating in 
laminar boundary layers may eventually be expected to lead 
to thermal dissociation of the diatomic species initially present 
in the free stream. With the onset of this strongly endother- 
mic process, the recovery temperature becomes temporarily 
insensitive to further increases in free stream Mach number, 
or, in the nonadiabatic case, maximum temperatures within 
the layer become lower than would be anticipated for flows of 
constant composition. Several investigators have addressed 
the question of the effects of this dissociation on boundary 
layer heat transfer for prescribed free stream and surface 
temperatures. Early numerical results of Moore (35) for flat 
plates in air indicated that, at reasonable wall temperatures, 
the heat transfer is essentially that which would be computed 
without taking dissociation into account.6 However, the 
equilibrium Prandtl numbers used by Moore were shown to be 
highly inaccurate in the temperature range of dissociation 
(18) and further calculations of the same type, but using more 
realistic Prandtl numbers, have since been reported (41, 34). 
More recently, Kuo (26) has treated the laminar boundary 
layer equations for a reversibly dissociating binary mixture 
and concludes that the change in the flat-plate heat transfer 
rate due to dissociation is expressible in terms of the cor- 
responding change in py, and in any case is not large. This 
conclusion is reported to be unaltered by the introduction of a 
more consistent Prandtl number for reacting mixtures than 
that used by Kuo (51). 


2. Turbulent Flow 


Chemically reacting turbulent shear flows, of importance in 
both external as well as internal aerodynamics, have recently 
received attention. Following Reynolds, the field equations 
governing globally steady reacting turbulent shear flows can 
be considered to be derivable, by a suitable time averaging 


® Moore implicitly assumed a Lewis number of unity. 
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process, from the nonsteady laminar flow equations. How- 
ever, as in the nonreactive case, the appearance of unknown 
“correlations” necessitates recourse to a semi-empirical ap- 
proach to turbulent momentum, energy and mass transfer 
processes. Lees (28) has discussed the extent to which exist- 
ing treatments for pure gases in turbulent flow are applicable 
to mixtures with chemical reaction. For fluids of arbitrary 
laminar Prandtl and Lewis numbers, the effects of chemical 
reactions in regions of the flow where laminar transport proc- 
esses are important must be studied in detail. On the other 
hand, extensions of treatments based on the Reynolds analogy 
(i.e., equal “exchange coefficients” ¢ for heat, momentum and 
mass transfer) to cover the case in which reaction takes place 
either only at the wall or entirely in the fully turbulent regions 
is relatively straightforward (2, 27). For fully developed 
internal flows or boundary layers with slowly varying free 
stream and wall properties, simple approximate expressions 
for the total heat transfer rate are obtained in terms of (a) the 
local skin friction coefficient, (b) local free stream and surface 
enthalpies and (c) laminar frozen Prandtl and Lewis num- 
bers. Corrections for the effect of variation of fluid proper- 
ties can be made, for example, by the reference enthalpy 
method. However, in the case of purely surface reaction with 
arbitrary wall catalytic activity, the wall enthalpy distribu- 
tion cannot be considered to be known a priori, since the con- 
centration distribution of atoms along the wall will be de- 
termined in part by the true kinetics of the surface reaction 
(1, 8, 45). The contributions to the heat transfer may then 
best be treated individually, with the mass-momentum analogy 
being used to obtain the approximate surface atom concen- 
tration distribution through application of the rate law for 
the surface. If the laminar frozen Lewis number is unity, or 
all reaction takes place in the fully turbulent portion of the 
layer, then the principal effect of reaction on the heat transfer 
can be accounted for by a modification of the effective “driv-_ 
ing force” (2). 
treatment of developing boundary layers with nonzero pres- 
sure gradient. 


Bartz (3) has utilized this approach in on 


Alternatively, for local equilibrium fully de- | 


veloped flows, the work of Schotte (48) and Brokaw (6) has — 


indicated that only a modification of laminar transport prop- 
erties and heat capacities is necessary to make existing heat 


transfer results applicable to reacting systems. For highly | 
‘cooled turbulent boundary layers in dissociation equilibrium — 


on blunt bodies, Rose, Probstein and Adams (42) choose to | 


semi-empirically modify results for laminar (frozen) Prandtl — 


and Lewis numbers of unity and treat the layer as incom-_ 
pressible as far as skin friction and velocity profile are con-— 


cerned. The boundary layer growth is determined from the | 


integrated momentum equation, and a modified Reynolds | 


analogy is then used to compute the local heat transfer in | 


terms of the skin friction and boundary layer thickness. — 


Their result has been used to correlate shock tube measure-— 
ments on a hemisphere cylinder. 


At the present time this writer knows of no theoretical study _ 


of the case in which reactions occur at arbitrary rates in the — 
partially laminar zone near the surface, although possible 
approaches to this problem have been indicated (2, 27). 
There is also an increasing need for experimental measure- _ 
ments under controlled conditions both to verify the trends 


indicated in recent analytical work and to suggest areas for — 


future development. 
VI. 
The extent to which analytical solutions to the equations of 
aerothermochemistry will approximate reality may be ex- 
pected to depend upon the reliability of the basic thermody- 
namic, transport property and chemical kinetic data used in 
actual computations. A brief statement of the state of affairs 
in each category is therefore given here. 


Input Data 


1. Thermodynamic Properties 


So long as pressures and temperatures are such that the 
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gas in question may be regarded as a mixture of perfect gases, 
it is possible to make accurate calculations of the over-all 
equation of state, heat capacities, enthalpy, etc., of equi- 
librium mixtures, or nonequilibrium mixtures of prescribed 
composition, using the methods of equilibrium statistical 
mechanics. Lighthill (30), for example, has considered the 
statistical mechanics of dissociating pure diatomic gases, and 
has further introduced the approximation of a three parameter 
family of “ideal” dissociating gases for purposes of investigat- 
ing the effects of dissociation on the dynamics of internal and 
external flows. For exact calculations involving air in ther- 
mochemical equilibrium, extensive tables and Mollier dia- 
grams are now available (14), and, with increased interest in 
other planetary atmospheres, similar diagrams may be ex- 


pected to appear. 
2. Transport Properties 


Because of the temperature levels involved, experimental 
data on the viscosity coefficient, thermal conductivity and 
diffusion coefficients for dissociating gas mixtures of current 
interest are virtually nonexistent. As a result, extensive use 
of Chapman-Enskog theory (21) has been made (4) in order 
to apply low temperature data for pure nonpolar gases to the 
calculation of high temperature transport properties of gas 
mixtures. However, the need for fundamental studies of 
atom-molecule interactions, of particular interest in com- 
puting transport properties of dissociated gases (28), promises 
to remain with us for some time. 


3. Chemical Kinetic Parameters 


a. Gas Phase Reactions 


With the recent development of aerodynamic methods (13, 
38) for determining chemical reaction rates, there has been 
increased activity in obtaining quantitative data on the rate 
constants for elementary dissociation and recombination reac- 
tions at high temperatures. In the absence of data for a par- 
ticular gas, or data in the appropriate temperature range, ex- 
isting rate theories have been used to estimate rate constants. 
For recombination reactions these include “collision” theory, 
and the theories of Eyring, and Wigner, a comparison of which 
has recently been given by Heims (20). It may be said that 
there is still a need for both a great deal of further experimental 
rate data, and a reliable theory. 

The behavior of inviscid flows with chemical reaction has 
recently received much attention (15, 20). These studies 
provide methods of estimating the free stream conditions to 


be imposed on the boundary layer equations for aerodynamic 
heat transfer calculations. 

b. Surface Reactions 


If for any reason a surface provides a favorable site for the 
recombination of atoms but is not itself chemically “‘attacked”’ 
it is said to act as a catalyst. A few quantitative studies of 
the catalytic activity of specific metallic and ceramic surfaces 
for the recombination of specific atoms have been carried out 
(31). Asa result it has been established that (a) virtually all 
recombinations of atoms on surfaces are first-order over a wide 
range of surface temperature, and (b) metals are invariably 
“better” catalysts than glass, Pyrex and silica surfaces. 
However, quantitative data on the magnitude and tempera- 
ture dependence of the rate constant for atom/surface com- 
binations of current aerodynamic interest are scant. In this 
case theoretical calculation of either absolute reaction rates or 
temperature dependences is less reliable than in the case of gas 
phase rates. 
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A Preliminary Study of the Application of Steady- — 


- State Detonative Combustion to a Reaction Engine 
R. DUNLAP,? R. L. BREHM? and J. A. NICHOLLS* | 


University of Michigan, Ann Arbor, Mich. 


The purpose of this study is to determine the feasibility 
of a reaction engine employing a continuous detonation 
process at the combustion chamber. A _ reaction-type 
engine employing steady-state detonative combustion is 
considered. A simplified analysis treats the supersonic 
mixing of fuel and air together with the requirements 
necessary to achieve steady-state detonative combustion. 
Calculations of specific thrust and specific fuel consump- 
tion as functions of flight Mach number are made for hy- 
drogen and acetylene fuels. The results of this study indi- 
cate that some supersonic diffusion of the air is necessary 
even though supersonic combustion exists. It is concluded 
that the speed range of air-breathing engines may be 
materially extended. 
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Nomenclature 
a = speed of sound 
A = area 
Co = specific heat at constant pressure 
F = thrust 
g = acceleration of gravity 
h = enthalpy per unit mass 
Cr = specific thrust 
= defined by Equation [20] 
Ke = defined by Equation [21] | 
m = molecular weight 
M = Mach number 
= pressure 
Q = heat release per unit mass | 
Ro = universal gas constant 
SFC = specific fuel consumption 
T = absolute temperature 
= velocity 
w = mass flow rate 
x,y = coordinate axes 
= ratio of specific heats 
p = density 
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é(M) = defined by Equation [13] 
n(M) = defined by Equation [15] 


Subscripts 
a, f = gases a and f (alsoa = air and f = fuel) 
= Chapman-Jouguet detonation 


e = exit conditions 
ig = ignition value 
n = normal component 
s = stagnation conditions 
t = tangential component 
= undisturbed air conditions 
1, 2, 3, 4 = stations in the engine 


theoretical work devoted to the study of detonation 
waves. It has been found that upon ignition of various com- 
bustible gas mixtures under certain conditions of pressure, 
volume and composition, a flame front will propagate through 
the mixture at speeds ranging from 3000 to 12,000 fps. This 
front (shock-initiated combustion) is called a detonation wave 
and the speed at which it propagates is the detonation ve- 
locity. So far this phenomenon has been observed mainly in 
shock tubes in which the wave is in a transient state. At 
present, however, attempts are being made to accelerate a 
combustible gas at high total temperature and pressure to 
the local detonation velocity and then cause it to ignite, thus 
establishing a standing (steady-state) detonation wave (1).4 

One is led to wonder if such a burning process could be 
applied to a propulsive system. It seems that a jet device 
utilizing this mode of combustion, as opposed to deflagration 
burning, would offer several advantages. For example, the 
supersonic inlet diffuser could be greatly simplified since the 
burning would occur at supersonic speeds. Thus the in- 
coming air need not be diffused subsonically and hence shock- 
swallowing problems, as well as large total pressure losses, 
could be eliminated. The static pressure rise, necessary for 
expelling the exhaust gases, would now be realized, at least 
in part, across the detonation wave rather than entirely 
through the diffuser. Furthermore, since the detonation 
process occurs at high velocities and total temperatures, it 
seems to be a natural means of extending the speed range of 
air-breathing vehicles. Other advantages would include a 
shortened combustion chamber with no need for an ignition 
device. An obvious disadvantage is the lack of static thrust. 

In view of possible applications, a preliminary analysis was 
made in an effort to determine some of the performance char- 
acteristics of an engine in which heat is added by means of a 
standing detonation. Special attention is given to the super- 
sonic mixing of fuel and air, and solutions are presented for 
the two cases when either the pressure or the area remain 
constant throughout the mixing zone. A general discussion 
of the detonation process, together with the method in which 
it could be applied to a steady-flow engine, is also given. 
Finally, problems associated with the matching of the vari- 
ous flow processes to produce an efficient thrust-producing 
mechanism are discussed and some over-all characteristics, 
such as specific thrust and specific fuel consumption, are 
calculated. 

It should be pointed out that the analysis contained herein 
is limited to the following general assumptions: The working 
fluid is assumed to be an inviscid perfect gas mixture. Aver- 
age values of the specific heats, applicable to the temper- 
ature range and mixture considered, are used throughout. 
Heat and frictional losses are assumed absent, and any total 
pressure losses due to the formation of oblique shock waves 
while compressing a supersonic gas are neglected. Finally, 
it is supposed that there exists an “effective ignition temper- 
ature” below which the fuel-air mixture must be kept prior to 
detonating. 


4 Numbers in parentheses indicate References at end of papers. 
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Fig. 2 Normal detonation wave 


In view of these assumptions, it should be clear that the 
results presented in this paper can only indicate approximate 
values of performance characteristics comparable to other 
idealized engines, which is all that was intended. 


General Considerations 


The following is a brief discussion of the flow system as- 
sumed in the analysis of this engine. Fig. 1 is a sketch of this 
flow system. 

The desirability of an engine requiring no diffusion of the 
incoming air has already been mentioned. At first one 
might think that no diffusion would be required in this engine 
since both mixing and combustion are occurring at super- 
sonic velocities. However, it will be shown that for maxi- 
mum pressure recovery through the engine there must be 
supersonic diffusion to some extent previous to mixing the 
fuel and air. This diffusion occurs between stations © and 
1. 

The fuel is injected into the supersonic stream at station 
1 and is assumed fully mixed with the air at station 2. 

Detonation occurs at station 3. Under certain conditions, 
namely, at low flight speeds, the back pressure should be suffi- 
cient to initiate the detonation normal to the flow at some 
station 3. However, this type of initiation may possibly be 
unstable, and in practice some form of stabilization, such as 
a thin body placed in the stream, will probably be required. 
At high flight speeds, certain mixing requirements and effi- 
ciency considerations dictate that a body must definitely be 
placed in the stream. The resulting oblique shock wave 
initiates the chemical reaction and the body serves as the 
means of stabilization. For the oblique wave portion of this 
analysis, the detonation is assumed to be stabilized on a two- 
dimensional wedge. 

Chemical equilibrium conditions are assumed to exist at 
station 4 immediately downstream of the detonation. Fi- 
nally, frozen equilibrium flow is assumed in the isentropic ex- 
pansion to atmospheric pressure at the exit. Downstream 
reflections of the stabilized wave are neglected and of course 
variable geometry is implied. 

It should be pointed out that the question of stability of 
the assumed mode of combustion is as yet an unresolved one. 
In fact, this is the primary objective of the Air Force contract 
supporting the study reported here. This problem is being 
attacked both theoretically and experimentally. For pur- 
poses of this paper, stability has been implicitly assumed. 
It is felt that this is highly probable for wedge stabilization 
although it is admitted that a normal detonation wave exist- 
ing between stations 2 and 3 would in all likelihood be un- 
stable and some means of stabilization would be required. 
Such a wave could not move all the way upstream to position 
1, however, due to incomplete mixing. 

In the over-all consideration of the analysis, two areas of 
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special interest and importance stand out—the supersonic 
mixing of fuel and air and the detonation process. Each will 
be considered separately in the following two sections. 


Considerations of the Detonative Process 


For purposes of analysis a detonation wave is usually 
treated as a discontinuity with heat addition.’ First, con- 
sider the one-dimensional model of a detonation in which the 
reference system is attached to the wave and the gases move 
relative to it (Fig. 2). 

The conservation equations for this system can be written 


4 


a3 a4 


mass: p3V3 = or 


momentum: P; + p3;V3? = Py + or 
+ = Py + [2] 


energy: + > +Q = + or 
2 = % 2 1 
—1 2 2 


Next, consider an extension of the above model to the two- 
dimensional case wherein an oblique detonation is stabilized 
on a wedge (Fig. 3). In this model the incoming velocity is 
resolved into components normal and tangential to the wave. 
Since there is no net pressure force in the tangential direction, 
there can be no change in the tangential velocity across the 
wave. The normal components may be treated as in the one- 
dimensional case. Hence the two-dimensional analysis re- 
duces to the vector addition of a (tangential) velocity to the 
one-dimensional normal wave model. Note that the con- 
servation equations for this model are identical to the previous 
one-dimensional equations except that V3; and Vy are now 
understood to be the normal components of the total veloci- 
ties. Note also that in the energy equation the tangential 
components of the velocity will cancel. 

One finds upon analyzing the above conservation equations 
that, for a given upstream temperature and heat release, solu- 
tions to the equations exist only for upstream velocities 
above a certain minimum. For all velocities above this 
minimum, the equations predict two possible solutions for 
the downstream conditions; however, one solution is ruled 
out by entropy considerations.® 

It is interesting to note that the minimum velocity is the 
velocity always measured in shock-tube investigations of det- 
onation. Associated with this minimum velocity are down- 
stream conditions such that the Mach number of the burned 
gases relative to the wave is always unity. Also, the total 
pressure loss across such a detonation is a minimum. This 
type of wave is called a Chapman-Jouguet detonation. 

Since total pressure is of prime importance in the efficient 
conversion of thermal energy into kinetic energy, the total 
pressure loss across the detonation in this engine was mini- 
mized by imposing the conditions necessary for a Chapman- 
Jouguet type wave. This proves to be a fortunate choice 
from still another standpoint since the heat release computa- 
tion can be eliminated. Thus, by using experimental values 
for the detonation velocity, the momentum and mass equa- 
tions can be solved without considering the energy equation, 
and hence the laborious chemical equilibrium calculations 
behind the wave are avoided. 


Mixing Analysis 


An exact solution to the supersonic mixing of fuel and air 


> It may be thought of as a shock wave with combustion. 
6 See (1) for a comprehensive discussion of the oblique wave 
solution. 
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would be prohibitively complicated. Consider, for example, — 
a scheme for introducing gaseous fuel into a supersonic stream — 
of air as shown in Fig. 4. Air, at a high stagnation temper- 
ature and pressure, would flow around the nozzles at super- 
sonic speeds and proceed to mix with the fuel through molecu- 
lar and turbulent diffusion. The resulting shear flow would 
be nonuniform and nonsteady. The flow would be further 
complicated by the presence of shock waves generated “ie 
the nozzle region, as well as by possible local burning in the 
boundary layer in the vicinity of the nozzle exit. Since it — 
would clearly be quite difficult to make a detailed investiga- 
tion, the following analysis will be directed toward a simpli- | 
fied solution of the mixing problem. 

Consider the mixing of two inviscid perfect gases (denoted — 
by subscripts a and f) as shown in Fig. 5. It is assumed that 
the rate of change of any flow parameter in the x or y direc- 
tion is zero at stations 1 and 2 (except at the initial interface 
of the gases where a discontinuity may exist in the y direc-— 
tion in the temperature and velocity), and that the gases are’ 
fully mixed at station 2. 

With these assumptions the conservation equations may be — 
written in difference form as 


mass: PiaViaA ta + pp Vip Ary = at 


Vie! vy" 
energy: ia hia + + wis + 9 
V.? 


momentum in .r-direction: 
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prA»V2? PiaAaV 10? pr Ay [6] 


Using the equation of state 


and approximating the enthalpy by 
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where C;, is an average value of the specific heat between 7’ = 
0 and T = T,, the conservation equations are rewritten in the 
form 
ViaM aAjaP; + yyMyAyPi [7] =constont [Mig 
Qe ay a» = constont 
2 4 
Cra Wia 
0.7 
M,*2.0 | P=constont 
A we | 
0.5 
PdA = (14 [9] | IN 
where 
0.4 
Ay = Ata + Ajy 
Notice that the momentum equation contains an integral Tp +2000 °R 
term that is dependent upon the variation of pressure during 0.3 
the mixing, while the other conservation equations express a 
relationship between the end points only. It is apparent 
that a solution can be found for constant pressure mixing, in 02 
which case 
PdA = P,(Az — A,) = — Aj).. {10} 
Ws 
or for constant area mixing where ee 7 
As 
PdA =0.. 2 3 4 5 6 ? 
Ai Mig 
In general the integral can be approximated by using an Fig. 6 Total pressure ratio across the mixing zone 


average constant value of the pressure during mixing although 
this case will not be considered. 

By combining Equations [7, 8, 9] and [10 or 11], the fol- 
lowing solutions for the downstream Mach number are found 
for constant pressure or constant area mixing 


It should be understood that the subscript for y in (1/) 
and $(./) should correspond to the subscript for 1. 
Of particular interest in the engine analysis are the down- 


constant pressure mixing: a, - = : stream total temperature, given directly by Equation [8], and 
Cpe tie the downstream total pressure. The total pressure ratios 
l across the mixing zone for the two cases are given for P; = 
(Mya) + — (My) by 
Y2 Mia Cora Wia Wia 
wher and for A; = A, 
1+ x 


constant area mixing: 7(M/.) = 


Are Ay 
1 — Mie? + Yy — Mi? 
A; A, 


us Cra 

| = Fig. 6 shows a plot of the total pressure ratio as a function 
Ons wis ~ (Mia)... [14] of M for a stoichiometric hydrogen-air mixture. For this 
' Yia Wis Aas (Mra) case, and in the remainder of the report, the subscripts a and 
Vis Wia Ana N(M;) f will refer to the air and fuel, respectively. The downstream 
yee _y. = boundary condition (7, = 2000 R) was chosen so that the eal- 
7 eos culations would correspond to a practical case wherein the 

combustible must be kept below an agzition | 
"ee es temperature. Curves for two different fuel-inlet conditions 

V soe 15] are shown. It can be seen that as the Mach number of the | 

air increases the total ratio continually decreases 
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and that in the Mach number range between M = 2 and M = 
4 the least loss in total pressure is obtained by injecting the eel 
fuel at the low energy level. It is physically plausible that — 

there should be a better total pressure recovery (referred to — 
the air total pressure) when the fuel is at a low energy level 
since more thermal energy from the air must be transferred 
to the fuel.7 As the Mach number of the air increases, the 
air total temperature becomes so large that the effect of the 
difference in fuel-inlet conditions becomes less significant. 

The results shown in Fig. 6 indicate the desirability of dif- 
= fusing the incoming air somewhat before injecting the fuel 
because of the importance of total pressure recovery in engine | 
performance. | | 
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Stoichiometric Ho -air | 
Tig = 2000 °R 
| T= 520°R 


Detailed Examination of the Flow System Mao 
f Fig. 7 Specific thrust vs. flight Mach number 

The integration of the previous work into an over-all engine : 
system will now be considered. 

An important limitation peculiar to this type of engine is | 
that the static temperature of the unreacted fuel-air mixture 
“a must be kept below an “effective ignition temperature.” 
The “effective ignition temperature”’ is herein defined as that _ 
temperature at which a moving combustible gas mixture will — 
— spontaneously ignite. Clearly, this temperature will de-— 
pend on the dynamic and thermodynamic state of the gas, _ 
composition, confining boundaries, the flow process involved, 
etc. To the authors’ knowledge, no studies involving all 
these parameters have been made. To effect the calcula- 

tions for this engine, it was necessary to assume some value of 
ignition temperature which appeared reasonable on the basis 
of known values found in the literature for stagnant mixtures. 

In applying the ignition-temperature condition, certain 
ramifications must be considered. The total pressure loss 
— across the detonation can be minimized by the attainment peren 
of a Chapman-Jouguet detonation at the lowest possible oa 
detonation Mach number (3). Since this Mach number | 
varies inversely as the square root of the temperature,’ it is 1 ! 
obvious that the detonation should occur at the highest pos- Moc 
sible temperature. This is advantageous from still another Fig. 8 Specific thrust vs. flight Mach number — 
M) standpoint since, as can be seen from Fig. 6, the total pres- 
sure losses through the mixing zone are decreased when the 
wh- Mach number of mixing is decreased. Hence, for best per- 
und formance, the static temperature of the unreacted gas must 
Lios be allowed to approach its limiting value, the ignition temper- 

P, } ature. Also, from this reasoning, it is concluded that super- 
sonic diffusion of the air prior to mixing with the fuel is is ae 
necessary. 

In locating the detonation in this engine, the above dis- 
cussion would indicate that it should coincide with the point 
of maximum static temperature, namely, at the end of mix- 
ing. Recall, however, that to achieve detonation, a certain CoMesir| 
minimum velocity is required. At low flight speeds this 
minimum velocity will not be realized at the end of mixing, 5 
and the flow must be expanded. The detonation will then 
occur at some downstream position (station 3, Fig. 1), as a 
Chapman-Jouguet normal wave. As the flight speed is in- 
creased, the conditions sufficient to generate this type of det- 
onation will occur further upstream until finally they will 
exist just at the end of mixing. For any higher flight speeds, 
the wave must be stabilized on a wedge so that the normal 
17 | component of velocity remains a minimum (Chapman-Jouguet 
oblique wave). 

With the general engine configuration thus determined, it is _ 2 
possible to compute over-all performance. By specifying the 
flight Mach number, ambient air temperature, fuel, fuel-air 
ratio and fuel-inlet conditions, the computations are made — = 
al- as follows: 
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T,g72000°R 
Tem 2520°R 
220 


All mixtures stoichiometric 
Tig = 2000 °R 
Te =520°R 
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Mi,=2.0 
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SFC ( 


on 7 Recall that in a frictionless heat subtraction process the total 
pressure of the gas will increase. = ; 
8’ The Chapman-Jouguet detonation velocities used in this re- Meo 


he port were assumed to be independent of initial temperature and 
eS pressure and were obtained from (2). Fig. 9 Specific fuel consumption vs. flight Mach number 
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(a) Setting T, = 7;,, determine M, and 7.2 from Equation 
8}. 
(b) Knowing M, and T.2, compute M,, from Equation [12 
or 14). 


(e) Calculate from Equation [16 or 17]. 
(d) Find = Mz and compare with V p. 
: 1. If V2 is less than Vp, use the normal wave solution 

7 below (Equation [18]) to compute V,. 
_ 2. If V. is greater than Vp, use the oblique wave solu- 


ae tion (Equation [19]). 

The following equations were derived by use of the conser- 

vation equations between the end of mixing and the exit. The 
solution across the detonation wave was determined so that 

= Chapman-Jouguet (normal or oblique) wave always 


occurred. 


normal detonation: 
is & 
= Ye x 
Vo 
2 
2 
2] _y-1 
2 2 


oblique detonation: 
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<= 


vit 1) . [19] 


where 


2 (as? + Vp?)? 
= 
k, 
2 bee, V ve [21] 


Performance characteristics such as specific impulse and 
specific fuel consumption can now be computed by the usual 


formulas 
Wia 


T => 

SFC = 3600 — wy 

Cr 


Results and Conclusions 


Figs. 7, 8 and 9 show the results of calculations of specific 
thrust and specific fuel consumption for hydrogen and acety]- 
ene fuels. In each case the fuel-air mixture was stoichiomet- 


ric and the effective ignition temperature was assumed to 
be 2000 R. The Chapman-Jouguet detonation velocities 
were approximately equal at a value of 5900 fps. The con- 
stant pressure mixing solution was used throughout. 

From these graphs, it can be seen that there is no thrust 
below a flight Mach number of 4 and that the specific thrust 
reaches a maximum between Mach 6 and 7, then decreases, 
The shape of the specific thrust curve can be explained as 
follows: First, below a certain flight speed the total energy 
of the fuel-air mixture is insufficient to achieve steady deto- 
nation. At slightly higher flight speeds a detonation will 
occur if the gases are expanded to a very high Mach number. 
However, for very high detonation Mach numbers, the total 
pressure loss is excessive and no thrust is realized. As the 
flight speed is increased, the Mach number of detonation 
decreases, and hence the specific thrust begins to increase. 
Finally the specific thrust reaches a maximum when the det- 
onation is stabilized at its limiting position (the end of 
a. For still higher flight speeds the Chapman-Jouguct 
oblique wave solution exists at a constant Mach number of 
detonation at 7’;,, and hence the total pressure ratio across the 
wave remains constant. Thus the specific thrust decreases 
because of the rising total pressure loss associated with the 
mixing process. 

The importance of fuel inlet conditions can be seen froin 
Fig. 7. The upper curve corresponds to M,,; = 2 and T,,;,; = 
2000 R, while the other curve represents the lower limiting 
values for these parameters. The fuel-inlet conditions are 
seen to be quite important in determining the operating range 
for this type of engine. 

The curves in Figs. 8 and 9 compare two different fuels. 
Since the detonation velocities were the same for both mix- 
tures and the effective ignition temperatures assumed equal, 
the differences in these curves are due to the unequal molecu- 
lar weights, fuel-air ratios, and specific heat ratios. 

In deciding upon the feasibility of the proposed engine, 
the specific thrust and specific fuel consumption may be com- 
pared to an ideal ramjet. The present engine offers com- 
parable performance at much higher flight Mach numbers, 
and indicates a possible means of extending the speed range 
of air-breathing vehicles. Finally, it must be stressed that 
these performance calculations are highly dependent on the 
effective ignition temperature assumed and are thus limited 
by the accuracy of this assumption. Furthermore, the re- 
sults have not been optimized from the standpoint of mixing 
process and fuel-inlet conditions. 
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Solutions to Stoolman’s External Diffusion Equation for = 
Instability of a Normal Shock Inlet Diffuser = 


CHIEH-CHIEN CHANG? and CHENG-TING HSU: 7 


Simple approximate solutions for the external diffusion ¢/v = yJ\M¢ = modified acoustic impedance 
regime of a normal shock inlet diffuser in subcritical oper- p(z,t) = p + p’ = instantaneous density 
ation are obtained in analytical forms. These solutions p(x) = time average of p(z,t) . 
correspond to Stoolman’s original differential equations = M*/M2* = u/u velocity ratio 
é = = shock oscillating velocity, toward 


in the external diffusion regime with an additional as- 
sumption of linear axial velocity variation. The solu- 
tions satisfy the upstream boundary conditions of flow Superscripts 


perturbations immediately behind the oscillating shock. 


downstream; £ is the amplitude v" 


Deviation of Stoolman’s numerical result from the present ’ = nonsteady small perturbation, function of x and ¢ 
one is explained. Acoustic impedance for the external 

diffusion region is also obtained, which indicates that flow Subscripts 

in the supersonic inlet-diffuser is intrinsically stable in : = section upstream of normal shock 


the absence of viscous dissipation. = section immediately behind the normal shock 
Nomenclature Introduction x 
A(x) = time-average flow area in external diffusion regime _ problem of transmission of ‘asada waves through © 
= local sound velocity _ a normal shock was first treated by Burgers (1),4 who — 
B,C = — to be determined indicated that the shock would oscillate under the interac-_ 
b = 7— = constant for gas property (b = } for air) tion of sound waves. Now, if this normal shock is detached — 
y+1 in front of the inlet-diffuser during subcritical operation (as — 
Cp, = specific heat at constant pressure and volume, respec- P 
wii gece —_— dais shown in Fig. 1), we are interested in knowing the amplitude | 
H(z,t) = eS/c: = h +h’ = instantaneous entropy ‘unction ratio of the reflected to the incident wave in the external 7 
h(x) = time average of H (z, t) diffusion regime when it impinges on the shock. The magni- — 
I(z) = function of z as specified in Equation [11] tude of this amplitude ratio (2, 3) is given as a function of — 
M = u/a = Mach number acoustic impedance 
M* = u/a* = dimensionless velocity 
P(z,t) = p+ p’ = instantaneous pressure 1+¢ 
p(x) = time average of P (z, t) 1-¢ 
R(z) = radius of cylindrical flow area A (x) : 
- = radial coordinate normal to diffuser axis Therefore the problem is reduced to finding the acoustic 
S = entropy impedance or other flow perturbation ratios® in the external 
t = time 
U(z,t) = u + u’ = instantaneous velocity component in z cicieatine é 
direction "DETACHED SHOCK 
u(x) = time average of U (z, t) ; 
w(2z 
w(y) = new dependent variable for v, see Equations [11, 13] i 
- = coordinate parallel to diffuser axis (positive toward i 
GOwnewenm) 
y = V/2bz = new independent variable for z +7 
2 2 
= = = (4%) = = new independent vari- 
a* leu | 
able for z ‘Ty | 
7 = cp/cy = ratio of gas specific heats (y = 1.40 for air) ty | 
¢(x) = time-independent parts of u’/u, p’/p, and p’/p, 
respectively; i.e., u’/u = vet!, p’/p = and 
p’/p = gett TIME 
w = angular frequency | position oF \ \ \ 
= ¢(x) — yo(x) = time-independent part” h'/h; 
i.e., h’/h = FLow AREA \ \EXTERNAL DIFFUSION | 
= — — = acoustic impedance 
yp u’ ; ss Fig. 1s External diffusion regime of a normal shock inlet diffuser 
Received Sept. 24, 1957. 
1 This investigation is part of a project sponsored by the Office ‘ Numbers in parentheses indicate References at end of paper. 
of Scientific Research, Air Research — Command, 5 The flow perturbation ratios in the forms of (p’ /p) )/(u'/u ), 
USAF, under contract No. AF 18(600)-1553. (h’/h)/(u’/u) were designated as density and entropy ‘‘transfer 
? Professor of Fluid Mechanics, Aeronautical Engineering functions,’’ respectively, by Stoolman (2) as followed Tsien’s 
Department. Mem. ARS. terminology (5, 6) used in ‘rocket nozzles and were later modi- 
* Research Associate, Aeronautical Engineering Department. fied as acoustic and entropy impedance in the forms of (a/yp)- 
Mem. ARS. (p’/u’) and (a/e,)-(S’/u’), respectively, by Mirels (3). 
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diffusion regim 

Stoolman was one of the first interested in the above prob- 
lem which is directly related to the instability of supersonic 
inlet-diffusers (2, 3). Stoolman has formulated linearized 
flow equations in the external diffusion regime with assump- 
tions of one-dimensional* quasi-steady, nonviscous, thermally 
-nonconductive, but compressible flow. Stoolman also 
imposed that the transverse velocity perturbation ratio 
v’/v is assumed to be small as compared to the longitudinal 
one u’/u on the control surface. Finally, he derived the 
linearized equations for the external diffusion regime as fol- 


 Jows ((2), p. 45) 


d 
+ Do - = @.... [2] 
u dz 


4s dv 


at: ou + 


where o and v are the two dependent variables while u, du/dz, 
a and dA/Adzx are known functions of x for a given geometry 
of the time averaging flow boundary. Stoolman obtained 
the solution by numerical integrations, first starting at the 
section immediately behind the normal shock where the 
boundary conditions of Burgers (1) were applied. This 
integration was extended downstream from the shock into 
the external diffusion regime, the flow area of which was ex- 
pressed in terms of Mach number distribution. In the pres- 
ent note, it is found that some simple analytical solutions 
of the above equations can be obtained approximately if 
linear velocity distribution in the external diffusion regime 
is assumed. 


Linear Velocity Distribution 


In one-dimensional steady flow, 
classical theory that 


it is well known from 


dA du 
ax —() — M*) 3 
Adz ) udx 
__ Introducing Equation [3] into Equation [1], we obtain 
do dv du 
dx dx 


In view of Stoolman’s experiments ((2), Fig. 10), the axial 
_ velocity gradient in the external diffusion regime is usually 
so great just preceding the flow instability that it is probably 
reasonable to assume that the velocity u decreases linearly 
with the distance x. Therefore we choose 


du u a* 


dz 
«a It can be further simplified if we use a new independent 
variable as introduced by Tsien (5, 6) 


\* u \? 
(2) - (3) [6] 
x a 


in a terms of which we « 


2 

Substituting Equations [5, 6, 7] into Equations [4, 2], we 
obtain, respectively 


‘an express 


dz dz 1 — 


bo  w «& 


® The justification of the assumption of quasi one-dimensional 
flow has been discussed by Stoolman (2) and Trimpi (4). 
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where b = (y — 1)/( 
Equations [8, 9] are two first-order equations in o¢ and », 
and can be combined to arrive at one second-order equation 


4 2 dp 4 1 — bz? 0 (10) 
—— — — | 
dz? 2z z—1/ dz 222(z — 1)(1 — bz) 


in v in which there are four singularities z = 0, 1, 1/6, and © 
corresponding respectively to J = 0, 1, and = in actual flow 
field. However, for the subsonic diffusion process behind a 
normal shock, only z = 0 and 1 of these singularities are the 
two limiting points, other singularities being outside of the 


range of interest in the variable z ae 
Solutions 


Equation [10] can be reduced to its normal form 


3bz + 5 5 8 8 
1622(1 — bz) ( + ) 


yp = —1)]..... [12] 


Since b = }< landz < 1, we can retain only up to the second- 
order term of z in the series expansion of J(z) without losing 
much accuracy. In this manner, we obtain the well-known 
Whittaker equation 


series rae (7) can A expressed as 


where B and C are two arbitrary constants, to be determined 
from the boundary conditions immediately behind the normal 
shock. 

It is felt that the above series solution of v is not con- 
venient for obtaining the solution for o in Equations [8 or 9]. 
By properly neglecting the high order terms of z in Equation 
[14], we take v in the following simple form 


The solutions for ¢ found from Equations [8 and 9] are, re- 
spectively 


B 1- 
o = — + — [1 — log(1 — + + const. 
7 ... [16] 
jeg 2b(1 — bz) 
(on 
) (2B 
2b 


The maximum difference between the second terms of the 
above two solutions is less than 2 per cent for air. There- 
fore it can be concluded that the first independent solutions 
of vy and o (containing the arbitrary constant B) as shown in 
Equations [15, 16, 17] are in close approximation and the 
second independent solutions (containing the arbitrary con- 
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[15] 


stant C) in the same equations are actually the exact solu- sat T 
tions. For later calculations, we will choose Equation [17] \\ 
for o. An interesting observation is that the solution of v “PRESENT “SOLUTION 
still retains two singularities at z = 0 and 1, while the solu- 24 \\ —-~ STOOLEAN'S MUMERICAL SOLUTION 
tion of o contains only one singularity at z = 1. 
4 
50 
The arbitrary constants B and C can be determined from A — i 
the boundary conditions immediately behind the oscillating & MA 
normal shock. Consider a normal shock oscillating slightly. ene ‘Ss : 
The upstream flow conditions (p, p, 7’) remain constant ex- E we WN 
cept for the magnitude of relative velocity, while the flow o 2.07— mse ee 
conditions behind the shock are disturbed and vary with the 
shock velocity ratio Assuming &/u. << 1, Stoolman == 
(2) obtained the following relations from Burgers’ general | 
expressions (1) by retaining only the first order terms in per- me ae aa aa aa = 
VELOCITY RATIO 


turbation 


1 
(“) /2 = (1 ... [18a] 


g 
¢ Me [18b] 
h’ 
= € = b)- 
(; Ue al... 
2 
bM, (M, 2)+ 6 [18¢ | 


YUM? — bf 


From the above boundary conditions, it is seen that the 
pairs y, and &, o2 and e are mutually in phase, but the former 
pair has a phase shift of 180 deg from the latter. For this 
reason, o»/(£/Us) and €/(£/us) are always expressed as nega- 
tive quantities. It is important to note that entropy will 
increase if the normal shock is moving upstream due to any 
disturbance such as isentropic compression waves originating 


in the plenum chamber. 
MIS 


Equations [15, 17] can be expressed in terms of the velocity — ie 

ratio \ and free-stream dimensionless velocity ./,* by mak- 

ing use of both Equation [6] and Prandtl’s normal shock 

relation = 1 


Final Results 


BM CM, *2(M,*? — bd?) [19] 
M,*? —\ A(M,*? — d?) F 
M,*? M,*? 
b 
&......... [20] 
M,*? — 2b 


Since for \ = 1, vy = vg and o = o», Band C ean be deter- 


mined from Equations [19, 20], i.e. 


2b 
B= — + (M,*2 — + 
—b 
(M,*? 
(22) 
2 
M,* \ M,*?— — bM,*)f 


where v2, ¢2 and & are given in Equation [18] and are fune- 
tions of only. 

Equations [19, 20] are plotted in Figs. 2 and 3, respectively, 
with A as the variable and ./, as the parameter. Stoolman’s 
numerical examples are also included in these figures for com- 
parison. Deviation of Stoolman’s numerical results from the 
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Fig. 2 Velocity perturbations in the external diffusion regime 
(quasi-steady ) 
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Fig. 3 Density perturbations in the external diffusion regime 


present analytical ones may be due to the following reasons 

1 Different entropy boundary condition used in Stool- 
man’s analysis.’ Stoolman considered that the entropy in 
the external diffusion regime is equal to €, multiplied by an 
area ratio of captured free stream flow area A, to the inlet 
diffuser area. 

2 Accumulated numerical errors possible due to succes- 
sive integrations in Stoolman’s results. The velocity per- 
turbation curve does not approach naturally to infinity at the 
singularity \ = 0 for Stoolman’s original curve and the J; = 
2 curve is above the 1/,; = 4 curve for Stoolman’s result as 
shown in Fig. 2. The above phenomena are probably due 
to numerical inaccuracy. 

3 Assumption of linear velocity distribution in the pres- 
ent analysis. The flow area A(x) used here may slightly 
deviate from Stoolman’s model. 

4 Approximation of the first independent solution in the 
present analysis. It is estimated within 2 per cent of error. 

Furthermore, the pressure perturbation ratio can be ob- 
tained from the simple relationship 


7 This is the main discrepancy from Stoolman’s analysis (2), 
whose conclusion leads to instability for supersonic inlet-diffusers 
even in the absence of viscous dissipation. 
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Therefore the modified acoustic impedance ¢/v or yM¢ is 
plotted in Fig. 4. 


Conclusion 


Tt is shown in the above analysis that the acoustic imped- 
ance is always a negative quantity. In other words, the 
amplitude of the reflected wave from the shock is always 
damped. Since the inlet is the only mechanism responsible 
for the instability of the diffuser (3), it is obvious that the 
flow is intrinsically stable without taking consideration of vis- 
cosity and unsteadiness. Further work will be published to 
show that the viscous dissipation is responsible for the in- 
stability of ae inlet diffusers even in a quasi-steady 
analysis. 
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Fig. 4 Acoustic impedance in the external diffusion regime 
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~The Optimization n of the N-Step Rocket With Different 


Construction Parameters and Propellant _ — 


The April 1957 issue of JET PROPULSION published the 
paper of M. Goldsmith “On the Optimization of Two- 
Stage Rockets,”’ giving the derivation of equations for de- 
termining the optimum weight distribution for two-stage 
tandem for the case of different construction parameters 
and propellant specific impulses in each stage. In the 
paper, the results were generalized for the case of the n- 
stage tandems, arbitrary number of stages, with different 
construction parameters and propellant specific impulses 
in each stage. All principal equations for the optimized 
two-stage rocket are easy to obtain from our equations for 
an optimized n-step rocket. The principal results pre- 
sented in this paper were worked out in 1955 and were 
presented by the author in May 1957 at the Conference 
on Rocket Techniques and Astronautics in Warsaw, 
Poland. The criterion of optimization used here is the 
minimization of gross weight for a given required burnout 


Received Aug. 26, 1957. 
1 Member of the PTA—Polish Astronautical Society. Address: 
ul. Weteranow Nr. 34-8, Lublin, Poland. 
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in Each Stage 


M. SUBOTOWICZ! 
Nomenclature 
n = number of stages 
Mou, . - -;Mon = mass of the Ist,..., nth stage 


m = useful payload; payload of the nth stage 

Mel, . - -» Men = empty mass of the Ist, , nth stage 

= An = Moing1)/Mon = = 
parameters of the Ist, . .., nth stage 


A = tI Ai = m,/mMo = total payload parameter 
0; = = Men/Mon = construction 
parameters of the Ist, . . ., nth stage 
Mpn = propellant mass of the Ist, ..., nth stage 
& = Mn/on, —n = Mpn/Mo, propellant 
parameters of the Ist, . . ., nth stage 


Daye) xp ee = burnout velocity of the Ist, ..., nth stage 

Wiy = «ny We = propellant exhaust velocity of the Ist, ..., nth 
stage 
n 


= = = total burnout velocity of the n-stage 


i=l 
rocket 
mor Mon 


Mon — Mpn 
., nth stage 


mor — Mp)’ 
parameters of the Ist, 
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The following general relations exist for the rocket parameters: 


1 1 


Derivation of the Optimization Relationship 
Between the n-Stage Rocket Parameters 


OR flight in a gravitationless and drag-free environment 
the following well-known equation applies 


n n n 
v= wigrn= — — &)]...[2] 
i=1 i=1 


The criterion for optimum or minimum gross mass is found by 
the method of Lagrange of the indefinite multiplicators. The 
total payload parameter A as a function of parameters @ and & 
according to [1] is 


n n 

A= A= & — 8)...... [3] 
We find now such a collection of the values of the propellant 
parameters &,..., € for which A has an extremum with the 
extra-condition 


n 
B= > 
~_ 


n 
[w; lg(1 — é,)] = 0. .... [4] 
1 


According to the Lagrange method 


En) + ., &.).. . 


and 


Using the relations [6] and [4] we can get the values of the 
multiplicator a@ and values of &, . . 
extremum values of A 


oA 1 
—e=—-— ff (1 —& — 6) = —A 
oA 1 
6, 
o(aB) 
B 
o&,, 
Substituting [7,] in [6] we get 
ob 
and hence | 
[91] 
In an analogous manner we get 
1 
1-&-@ 1-& [92] 
1 
= 0.. [Qn] 


To get &,..., &n_1) as a function of &, only, we divide [9,.] by 
[9:], then by [92], and so on to [9(n-1)]. We get finally 
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Which indicate the 


6,w; (1 — 
OnWn + (1 En) (w; Wn) 


(I — fn) 
+ (1 — En) (won-1)) Wn) 


1 — = 


The equations [10,;] through [10,,1)] are the optimization 
relationships? sought between the n-step rocket parameters 
with different construction ratios (, ..., 8, and propellant 
exhaust velocities w), ..., Wr. From these relations for w; 
= W. = ... = wy, we can derive the Malina-Summerfield*® 
conditions for the n-step rockets with equal parameters. | 


Derivation of the Principal Equation of the 
Optimized n-Step Rocket 


To get the value of &, we put [10,]—[10,,_1)] in [4] and find 
£, with numerical methods. Then it is easy to obtain the 
values of &, . . ., &n-1) from relations [10,]-[10,,1)]. But 
this method proposed by Vertregt is rather difficult. 

We show another simple graphical method to get the op- 
timized parameters of an n-step rocket. This method is more 
general but analogous to the Goldsmith method for a two- 
stage rocket. 

At the beginning we shall treat our problem for the case of 
a three-step rocket. But there are no principal limitations to 
extending our method for the case of the n-step rocket. 

The burnout velocity of the three-step rocket according to 
[2] is 


v= — [wilg (1 + we Ig (1 — &) + ws Ig (A &3)]. . . [23] 
The optimization conditions for a three-step rocket according 


to [10,]-[10,,_1)] are 


1-§&= ..... 
&;) 


= 10 
+ (1 &)(we ws) 


1-&= 


Our conditions [10,a], [10.a] can be written according to [1] 
as 
+ 63) 


10,5 
+ (As + 02)(wi — ws) 


6. As + 4 ) 


+ + (As + 6; (we — ws) 
From [10,6] and [10.6] we obtain 
3 2 6. 
= + (wi — ws); ws = We xt (we — ws) 


We shall express \; and from [11] 


ws; 


\ = 
Ww) 


— + (we — 
3 


We put the parameters \; and A» from [12] in Equation [23], 


where the expressions with (1 — &) are replaced by (A + 6) 
according to [1] 


2 An analysis analogous to the present calculation has been 
published by Vertregt (Journal of the British Interplanetary So- 
ciety, March-April, 1956, p. 95) who shows an explicit solution 
may be obtained for the case of equal exhaust velocities. 

3 Malina, J. F., and Summerfield, M., Journal of the Aero- 
nautical Sciences, August 1947, p. 471. 
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: 
ws); 6, we 6, We 
w, — + — ws) 
Wn-i 6, Wn-i 
We lg + 0. + lg (As + 63) [13] Wn n Wn 
3 
wot + (ws — ws) The relations [18,]—[18,_1] can be easily derived from the op- 


timization conditions [10,]-[10,-1]. We can find the other 


After simple calculations we get from [13] for a ve _ parameters of the separate stages of the step rocket by the 


cana ie application of the relationships [1]. We find the starting mass 
—v — w lg 6, — we lg 6 — ws lg & = 7 of the step rocket by using the equation 
+ we + ws) lg (: +%)- w, lg [ + (1 =) | my = mA! 
3 3 > 


we Ig [1 + (1 where A 
We 6; 


By using the derivated optimization conditions [10,]—[10(.-1)] 
for the n-step rocket we obtain in an analogous way the follow- 
ing general equation for the optimal n-step rocket 


( is 8.) (2 ) +3) 
n—1 


where i = 1, 2,...,”. This is the principal equation of the 
optimal n-step rocket. 


LF (5/4) 


Graphical Method of Obtaining the Optimal a 
a 7 Parameters of the n-Step Rocket 
> I 


2quation [15] has a suitable form for a fast determination 
of the optimal parameters of the n-step rocket with different 
parameters of each stage: Ai, &:, Ti, Wi. 

The left-hand side of Equation [15] contains only the pa- 
rameters given earlier: v, w;, 6;. The right-hand side of this 


equation is the function F of (A,/6,); An and @, are the 
parameters of the last stage. We can get the optimal parame- 
ters by using Equation [15] in the following way: The right- Fig 1 Curves drawn for y = 0.5 and different values of » ac- 
hand side of [15] cording to Equation [20] 

F(An/0n) = (2 w.) lg (: + = (3 


is used to draw a curve with values of (X,,/6,) on the abscissa 
axis and the numerical values of’ the function F(A,/6,), ac- 
cording to [16], on the ordinate axis. The values of w; (¢ = 1, 
2, ..., m) are given in practical cases. After drawing the 
curve F = F(X,/6,) that is suitable to the given practical 
problem, we shall calculate the value of the left-hand side of 
Equation [15]. This value, denoted by A,, is an exact defined 
number which depends on the burnout velocity v, on the pro- 
pellant exhaust velocities of each stage w,; and on the con- 
struction parameters 6;. We treat v, w; and 6; as given. We 


put the so-caleulated value 
n 
A, = —v— (wig 9;)............. [17] 
i=1 
on the ordinate axis /'(A,,/6,) and look for the suitable abscissa ; - 
we find it by drawing a line parallel to the point of intersection 7 
of this line with the curve F = F(X,/6,). The abscissa of <= 
this point gives the value (A,,/8,). From this we can quickly =. 
derive the payload parameter \, because the value of 6, is a 
given. 
After the calculation of (A,,/6,) we obtain the other ratios 
(Ai/81), (An-1/8n-1) from the relations 92 os os a 
A _ wi On (2 wt [18)] Fig. 2. Curves drawn for q = 0.8 and different values of p ac- 
Wn An \wn cording to Equation [20] 
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Example of the Three-Stage Rocket 


We shall illustrate this method with the three-step rocket. 
We introduce the notations 


p = q = w2/ws; [19] 


The function F(\,/6,) divided by w, we shall denote by 
I'(X,/9,), and the values of A, divided by w, we shall denote 
by A,. Equation [15] for the case of the three-step rocket has, 
by using [J9], the form 


ie — plg@ — qlg & — lg 0; = A; = 
W3 


1 As 

6; 6; 

q 6; 
The curves drawn for different values of p and q¢ according to 
iSquation [20] are given in Figs. 1 and 2. We can easily ob- 


tain the optimal values of the parameters for the suitable p 
and g according to the method given above. 


Derivation of the Goldsmith Equations 


Goldsmith’s equations for determining the optimum weight 
distribution for a two-stage rocket for the case of different 
construction parameters and propellant specific impulses in 
each stage are derived from the theory presented here for the 
n-stage rocket. Goldsmith’s symbols and ours are different. 
The calculations given by Goldsmith are conducted for the 
other parameters as done here. The optimization condition 


calculated according to the method used by Goldsmith for 


the parameter 6 is written in our notation 


or, using Goldsmith’s notation 
Ma 
19, ( 1 [21a] 


This equation is analogous to Goldsmith’s Equation [7]. 
The condition [21] is a special case of the optimization 
condition [10;]—[10,_1] or [18,]-[18,1]. 
The principal Goldsmith equation written in our notation 
for the parameter 6 (n = 2) is 


—v —w lg — we lg 


(wy + ws) (1 + lg +(1 
w; 


The only difference between this equation and the Goldsmith 
equation is that our equation is written with respect to the 
ratio (A,/0,) of the last step; Goldsmith’s equation is written 
with respect to the ratio (1/6) of the first step " 


. [22] 


v We 
we, wi\.. 
lg |! (: [22a] 
Ww) We 


where X = ),/6;. 
It is evident that Equation [22] is a special case of our 


general principal, Equation [15], forn = 2 


On the Optimization of Multistage Rockets — 


H. H. HALL! and 


Aeronutronic Systems Inc., Glendale, 


The optimum weight distribution for multistage rockets 
having different specific impulses and structural factors in 
each stage is derived. Minimization of gross weight for a 
given required burnout velocity and payload is the criterion 
of optimization used. A method is suggested for including 


in first approximation the effects of gravity, drag and turn. 


Introduction 


ARIOUS authors have treated special cases of the prob- 
lem of optimizing the performance of a multistage rocket. 
Malina and Summerfield (1) treated an arbitrary number of 
stages with the two restrictive assumptions that structural 
factors and specific impulse were the same for all stages. M. 
Goldsmith (2) recently treated the case of a two-stage rocket 
having different structural factors and specific impulses in 
stage.4 
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Another paper by Schurmann (4), appearing very recently, 
contains the solution of a related problem. In it, the maximum 
accelerations to be achieved in each stage are prescribed in ad- 
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The purpose of this paper is to present the solution to the 
general problem of a rocket having an arbitrary number of 
stages, and having different structural factors and specific im- 
pulses in each stage. 


F P 


As a starting point for this analysis, consider: ation is re- 
stricted to the usual highly simplified proble m, with the same 
three idealizations expressed or implied in (2). These are: 
(a) Thrust direction is constant (no turning); (b) thrust is the 
only force acting (no drag or gravity); (c) specific impulse is 
constant throughout any given stage, but may vary from 
stage to stage. 

With these assumptions, the mathematical complications 
ure stripped away. To further increase the practical value of 
the method, the notation is also condensed to the minimum 


vance, ‘and, for a solution, the specific impulse must be taken 
equal in all stages. 

In the treatment presented here and in (3), the specific im- 
pulses for each stage are prescribed, and the maximum accelera- 
tions depend on the mass ratios and the burning time, the latter 
being essentially a free parameter. Cross comparisons of the 
method of (4) with that of (3) and this paper are useful since 
conditional acceleration constraints often exist in real problems, 
i.e., the acceleration must not exceed a certain limit, but may be 


smaller if a rocket size advantage results. 
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number of parameters required to express the problem. 
Stages are numbered by the index j, running from j = 1 to 
j = J, the total number of stages to be considered. 


V = velocity at burnout of last stage 

C; = exhaust velocity of jth stage, often denoted as V,,; 1 

= payload weight 

W,; = weight of propellant burned in the jth stage 

W.; = gross weight of the rocket just before the jth stage is ig- 
nited, and after all expendable material from the pre- 
vious stage has been staged off 

W,; = weight of material staged off following burnout of the 


jth stage 
tj; = ratio of the mass of the missile at ignition of the jth stage 
_ to the mass of the missile at burnout of the jth stage 
Gi 
Woj Woi 
8; = the disposable structure factor, representing the fraction 
of a stage which is not burnable propellant 
é 


+ Wo; 


Most of the notational economy is attained by recognizing 
that the solution to the problem depends on the total weight 
staged off each time, but not at all on the detailed distribu- 
tion of that weight among such contributors as structure, 
motor, fins, unused propellant, telemetry equipment, or any 
of the countless other items, each of which in practice turns 
out to have its own intricate scaling laws. 

It is most convenient to deal directly with payload P, and 
gross initial weight, Wo. Experienced designers often re- 
place W,,, W,; and Wo; (propellant consumed, the weight 
staged off, and the gross stage weight, respectively) by the two 
equivalent ratios, 4; and 8;, for which they may have a great 


deal of intuitive “feel.” 
The gross initial weight of the rocket is = 
Wa = Wa + Wa + We + Wet....+ ht 
Wis + Wor + P...... [3] 


Using Equations [1 and 2], it is a simple matter to rewrite 
Equation [3] as 


(i+ 4 
[4] 


j=l 
where the notation indicates the product of factors of the form 
shown. 

The above expression is to be minimized subject to the 
condition that the rocket attains a prescribed burnout velocity 
V, at the end of the last stage powered flight. 

Velocity V is given by Equation [1] of (2), which in our 


notation takes the form 5 


If expression [4] is a minimum, then it may readily be 
shown that its logarithm is also a minimum 


Thus the most convenient statement of the problem is that 
Equation [6] shall be minimized by suitable choice of the y,;, 
subject to constraint [5], where the C; and 8; are arbitrary 
sets of J constants each. 

Solution of the Problem 
‘r= 


The minimization is accomplished most directly by using 
the method of Lagrange multipliers. 


the Lae 
grange multiplier \ in Equation [5] 


Summing Equations [6 and 7] 


j=l 


j=l 


J J 
j=1 j=1 

Differentiating with respect to u; at the extremum, and 
solving for yp; 


AC;B, 
Equations [9 and 5] provide a system of J + 1 equations in 


the J + 1 unknowns, A and J values of u;. To solve this 
set of equations, take the natural logarithms of Equation [9], 


multiply by C;, and sum over j, obtaining = 


j=l j=l ACB; 
Hence, from [5] 
J AC; +1 


Equation [10] is a transcendental equation in a single un- 
known, A, and is easily solved for by interpolation tech- 
niques. Having evaluated X for a given set of C; and y;, the 
uw; are found by substituting into Equation [9]. 

This solution holds for any number of stages and, when 


J = 2, contains the results of (2). ai 
= 


Two special cases are described below, where it is anal 
to explicitly solve for the optimal mass ratios. S" « 


Discussion of Results 


Special Case 


For the special case where C; does not depend on J, i.e., the 
exhaust velocity is the same in all stages, Equation [9] be- 


comes 


Taking the natural logarithm of Equation [11], summing 
over all 7, and using Equation [5], one obtains 


- Substituting into [11] one has the solution 
B; 


A direct physical interpretation is found by noting that 
Equation [11] indicates that the product u;8; has the same 
value for all stages. Thus the optimum mass ratio yu; is 
large where the structure factor is favorable (i.e., small) and 
vice versa if the exhaust velocity is the same in all stages. 


Special Case 2 


If neither exhaust velocity C, nor structure factor B, depends 
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on j, Equation [12] becomes 


V 


- 
= 
; 
y . 
A 
: 


[7] 


ible 


1s 


This is the result first presented in (1), that the optimal igni- 
tion-to-burnout mass ratios are the same for all stages. 


Inclusion of Gravity, Drag and Turning Corrections 


To the extent that one is able to make rough estimates of 
velocity losses due to gravity, drag and turning, these losses 
can be included as modifications of the required V in Equation 
|5], and the corrections will be reflected in the determination 
of the u;. The resulting improvement is, of course, no better 
than the estimates of the losses; hence this method is not a 
substitute for accurate calculations with a high-speed digital 
computer. Instead, it has proved to be a valuable auxiliary, 
in providing very good slide-rule designs, thereby materially 


reducing the amount of expensive computer time needed to 
locate design optima. 
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The problem considered is that of determining rocket 
thrust programs that use the least amount of fuel to propel 
a given mass in a vertical plane from one point to another 
without aerodynamic lift, where gravitational force and 
aerodynamic drag are known as functions of altitude and 
speed. It is solved by formulating it as a Mayer problem in 
the calculus of variations with consideration given to the 
fact that the rocket mass can never increase during the 
flight. A sample problem for a short-range ground-to-air 
rocket has been worked out using a digital computer to 
show some of the features of these RAN Senjnetorton. 


Introduction 


NODDARD (1),3 Hamel (2), Tsien and Evans (3), and 
Leitmann (4, 5) have considered the optimum thrust 
program for a vertical sounding rocket and have established 
that this program involves a rapid boost at the beginning of 
flight, usually followed by a period of continuous burning (sus- 
tain phase), and ending with a zero thrust period (coasting 
phase). Edwards (6) and Cicala and Miele (7) have formu- 
lated the present problem in considerable detail; the principal 
contribution of the present authors is to add to this formula- 
tion a way to handle the constraint of nonincreasing mass and 
some numerical calculations of a typical case. Actually, 
two constraints were considered: (a) that the rocket mass 
never becomes less than the final (empty) mass and (b) that 
the rate of change of rocket mass be zero or negative. It was 
found that specifying only (a) was sufficient for the present 
problem since the resulting thrust programs also satisfied (b). 
Earlier writers have considered calculus of variations prob-— 
lems with such constraints (8, 9). 


Formulation of the Problem 


We will consider the problem on a flat earth and neglect | 
variations of the gravitational force with altitude, which is an 
adequate approximation for short-range rockets. (The more 
complete problem is a direct extension of the present one.) 
Adopting the coordinates shown in Fig. 1, the conservation of 
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momentum parallel and perpendicular to the flight path re- 
quires that 


i. my = —cm — D(y,v) — mgsiny........... [1] 


vy = —gcosy.... 
where m is the instantaneous mass of the rocket, D the aerody- 
namic drag and v the velocity. The rocket has been assumed 
to have high weathercock stability so that its axis is always 
parallel to the flight path. (The flight path is the path of the 
center of gravity of the rocket.) 


_ We will also need the kinematic relations 


; The thrust of the rocket motor has been given as —cm, 
where c is the effective exhaust velocity of the gases relative 
to the rocket and is assumed to be constant. Note that c/g is 
the specific impulse of the rocket propellant. Referring to the 
introduction we see that we need to require that 


and, further, we should check our solutions to be sure that 
m < 0 throughout the flight since fuel cannot be “sucked 
back in”’ to the rocket. 

Now, let us define an “admissible trajectory” as one that 
satisfies Equations [1 through 5] and in addition has certain 


ox 
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specified final values of m, x and y, and certain specified initial 
values of v, x and y. 
The problem, simply stated, is to find the 


jectory that minimizes the initial mass. 


admissible tra- 


The Equations for a Stationary Trajectory 


We first change from time to horizontal distance, x, as in- 
dependent variable using Equation [4] 


dx 
dt = — 
v cos Y 
Next we introduce dimensionless variables as = yoeeeee) 
gx m v gy gt De” 7 
c my c c msg 


Next we eliminate uw as a variable in favor of the quantity 
o = pe’, which does not change during an impulsive burning 
of fuel (see Equation [1]). Log @ = v + log uw has the proper- 
ties of a “‘potential velocity”’ since it is the velocity that could 
be attained at any point in the flight by burning all the rest of 
the fuel instantaneously, i.e., decreasing uw to 1. Such an in- 
stantaneous “boost,’’ while not possible, is a convenient 
idealization to a very rapid burning of fuel. Rewriting Equa- 
tions [1, 2, 3, 5] in the new variables 


se 


J, =o’ +- sin y + Q) = 0 where Q = Q(v, 7) 
[7] 
1 
y? 


where ()’ = d()/dé and we have introduced a new real vari- 
able ¢ in Equation [10] to insure that 4 > 1. These are four 
equations in five variables, ¢, y, 7, v and ¢ so one of them is 


arbitrary; we shall regard 7 as the arbitrary function. Note 
that the time ¢ can be computed from 
[11] 
v 


Consider now an admissible trajectory, characterized by 
the functions ¢, y, , v and ¢, and another admissible trajec- 
tory characterized by @ + 6¢, y + dy, 7 + 6n, v + dv and ¢ 
+ 6¢ where 6() signifies a small variation. Since both are 
admissible trajectories, it follows that 


OJ; = + Og, + Oy, + bn, + dy, + 


Therefore 


1 
é, 
i=1 
where A,(£), 7 = 1, 2, 3, 4, are four more dependent variables 


known as Lagrange multiplier functions, as yet unspecified. 
Integrating Equation [13] by parts in the usual manner, we 


find 
= + Ady + — 
[ Ai sec? Ai sec? Y (@ + Qsin y) + Az sec? dy + 


sin y + Q — 


od’ ov’ on 


Now, we choose As and A, so that the coefficients of 5¢, 
by, dv and 6£ all vanish, i.e. 
tan y 


de’ — Doin y) + As sec? = 0....[16] 


2X2 cos ‘Y 


sin y + Q — vQ,) + — cos-yAye” = .. [17] 


[18] 


Referring to the terms before the integral in [14], we see that 
for an admissible trajectory [dy]: = [dy]:-:, = 0. Now 
(0) and y(£;) are unspecified but we may choose 


Similarly, since v(£,) is unspecified, it follows that $(&,) is 
unspecified, but we may choose 


Equations [19a, 19b] are known as “‘natural boundary con- 
ditions’ of the problem. Because of our choice of \’s through 
Equations [15-19], Equation [14] reduces to 


v 


From [18] it follows that at any given time during the flight 
either A; = Oor ¢ = 0. From [10] it can be seen that ¢ = 0 
corresponds to m = my, i.e., constant mass which means zero 
thrust. Now we have used the variable ¢ for the express pur- 
pose of allowing for an instantaneous boost at & = 0 since 
o(0+) = o(0—), i.e., @ does not change during such a boost. 
Since we wish to find the minimum initial mass u(0—) for a 
given initial velocity and = (0), it will 
suffice to minimize $(0). Referring to Equation [20] we see 


that if we put 
= const 0.............. 
we can obtain a stationary value of $(0) for arbitrary small 
6n if 
hi? — [22] 


Since all of the equations involving the \’s are homogeneous in 


them, we may set 
= 


An admissible trajectory that gives a stationary value of 


initial mass will thus be obtained by the simultaneous solution 
of the seven first-order ordinary differential equations [7, 8, 9, 
10, 16, 17, 22] plus the two equations [15, 18] in the nine de- 
pendent variables ¢, n, Y, 
boundary conditions 


Au, Ae, Az and with the seven 


Equations [15-18] plus [22] are the Euler-Lagrange equa- 
tions corresponding to the integrand of Equation [13]. Since 
the integrand of [13] is not an explicit function of &, a first 
integral of the five Euler-Lagrange equations exists, namely 
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where 
G = Aid 
i=1 


first integral becomes 
de et 
As = etn y + + Qese y) + — 
v 2 
where k, = const. Substituting [26] into [16] gives 
cos Y 


de’ = sec Y ese Y (« - 
y? 


A Second Integral for the Sustain Phase 
If we use an exponential law for the variation of drag with 
altitude, i.e. 
Qy, n) = . [28] 


then it follows that Q, = —@BQ. Using this in Equation [16] 


euation if \; = 0 (as it is during sustain) which upon integra- 
tionyields 


— \; = ke = const............ [29] 


Relations for the Sustain Phase 


During the sustain phase \; = 0, so we have three linear 
simultaneous algebraic equations for dy, A», As, namely [17, 
26, 29]. Solving for A;, we obtain 


2v(ki cos y — ke sin y) 


= 
(1 — 2Bv)d sin y + 24 7Q, 


Substituting this into [15] with Ay = 0, the following first- 
order ordinary differential equation for the velocity v is ob- 
tained 


This last equation, solved simultaneously with the equations 
of motion [7, 8], and the kinematic relation, [9], determines 
the sustain phase of the optimum trajectory. From [17, 24] 
the boundary condition A,(0) = 0 yields 
[Psin y + Q — vQ, = O.. 
Utilizing [32] and \,(0) = 1 in [380] we find that 
((Bv — — BQ + cos y — ke sin y leno = 0 ... [30] 
Substituting [32, 33] into [31] we obtain at & = 0 


Equation [32] can be thought of as determining (0) in terms 
of initial ¢, y and y. Equations [33, 34] can be thought of 
as determining the constants k; and ky in terms of initial @, 7, 
y and v’. 


Boost Velocity 


Equation [17] with A, = 0 provides us with a way to caleu- 
late the amount of the initial boost since u° = poe”, if vo = 0, 
where the superscript refers to conditions before the boost and 
the subscript to conditions after an instantaneous boost. 
Substituting this into Equation [17] we have 


2(0, vo) — vo) + wo sin Yo =O0........ [35] 
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Substituting into [25] and multiplying through by ctn y, the 


. . [26] 


[27] 


and multiplying [7] by \; and [22] by @and adding the result-_ 
ing three equations together, we obtain an exact differential — 


0) k tan y[(Bv — 1)Q, — BQ] + 2vk, ese y — (Q + vQ,) etn 7] 


where Q = Q(n, v). Thus if w°, yo, and the drag function Q 
are given, Equation [35] is a transcendental equation for v. 
Fig. 2 shows a plot of this relation for the particular case of a 
quadratic drag law 


[36 


where Dy = drag at sea level at speed v = c. 

Note: with y = 90°, Equation [17] is the optimizing con- 
dition for the entire sustain phase of a vertical sounding 
rocket, a simple relation apparently missed by Hamel (2) and 

_Tsien-Evans (3), but noticed by Leitmann (4). 


> 


- Conditions at the Juncture of the Sustain Phase and © 
the Coasting Phase 


By considering the conditions at the juncture of a con- 
tinuous burning are and a coasting are (both extremals of 
the problem) it is easily determined that all the \’s as well as 

.v, y, nand ¢ are continuous across this juncture (see Bolza 
(10)). The first derivatives of these quantities may be discon- 
tinuous across the juncture. 


Relations During the Coasting Phase 


During the coasting phase, u = 1, and @ = e’, so that Equa- 
tion [7] becomes simply 


—sec 
= (Qe-’ + sin y) 
v 


Also, eliminating A; between [15, 22] vields 


_ Ai tan y 


Ai’ [y—1 ese ¥(Q — vQ,)] .. 8] 


During the coasting phase then, Equations [37, 38] must be 
solved simultaneously with [8, 9, 27] with initial values of , 7, 
y, A; and ), that are the final values of the sustain phase. 


, Avsin y{Q + vQ, + v2[(By — 1)Q2, — BQ] — — 1) ese + kz cos y{ 2 + vQ, — v? tan? [(By — 1)2, — BQ]} 
— (d@sin y + Q — cos — ke sin 


[31] 


Method of Solution 


It can be shown that the present problem of finding the 
thrust program that minimizes the initial mass of a ballistic 
rocket in traveling between two given points is identical to_ 
the problem of maximizing the range of a rocket with a 
given amount of propellant when the initial point and the final 
altitude are given (or maximizing final altitude with final — 


range given). This latter viewpoint was taken in our sample 
problem. Thus we start with a given initial mass, take an — 
arbitrary y(0), then use [35] to determine = v(0+) and 
my = (the velocity and mass after boost), then proceed 
with the sustain phase by guessing an initial acceleration 
v'(0). The problem is then integrated numerically by finite — 
differences throughout the sustain phase until the fuel is ex-_ 
hausted; at that point, coasting flight is begun (zero thrust) _ 
and it is continued until A. = 0. If A; is not also zero at that — 
point, the sustain phase is started over again with another | 
value of initial acceleration v’(0). It was found that the | 
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Wo= WEIGHT BEFORE BOOST 
Do= CpS = DRAG AT SPEED C 


c = EFFECTIVE EXHAUST VELOCITY 
OF ROCKET GASES RELATIVE TO ROCKET 


Yo = ELEVATION ANGLE AT LAUNCH 
Vo= VELOCITY AT END OF BOOST 


Boost velocity for maxi- 


range ballistic trajectory- 
quadratic drag law 
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Fig. 4 Altitude-range envelope for optimum thrust programming 
of ground-to-air rocket for sample problem 


W,= WEIGHT BEFORE BOOST 
2 
Do =CpS 
4 C = EFFECTIVE ROCKET EXHAUST VELOCITY 
Yo = ELEVATION ANGLE AT LAUNCH 
Wosinye To= INITIAL SUSTAIN THRUST 
Po Fig. 3 Approximate initial sus- 
tain thrust for maximum range 
ballistic trajectory (quadratic 
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Fig. 5 Thrust vs. time curves for optimum thrust programming 
of ground-to-air rocket for sample problem 
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correct value of v’(0) could be estimated rather closely by | 
differentiating the coefficient of A; in Equation [17] (see Fig. 
3). If the final altitude came out below sea level, a higher 
value of y(0) was used on the next try. In this way, by tak-— 
ing the final value of 7 that came out of the problem, we say 
that if we had specified this value as 7; then the thrust pro- 
gram of this problem meets the necessary conditions for a 
stationary value of final range, &,, for that altitude. 


Sample Problem 


A sample problem was calculated using the Harvard 
Univae. The parameters for the problem were 


m'g = 1000 lb; msg = 500 lb; p = poe CpS = 0.47 ft? 


and the drag law was chosen (for simplicity) to be 


For altitudes below about 50,000 ft, 1/a = 30,000 ft approxi- 
mates the density variation with altitude fairly closely. 
The problems calculated were “inverse problems,”’ i.e., the 
maximum range at various altitudes for a rocket with a mass 
ratio of 2 and certain drag characteristics was calculated. 
Fig. 4 shows five optimum trajectories. The envelope of 
these trajectories represents the maximum distance this rocket 
can reach with the best possible thrust programming. Note 
all the trajectories consist of an initial instantaneous boost, 
a sustain phase, and a coasting phase. The initial flight path 
angle for maximum range ground-to-ground appears to be 
about 68 deg, well above the value of 45 deg which it would be 
in the absence of drag. Fig. 5 shows thrust vs. time during 
the sustain phase; note about 32 per cent of the rocket fuel 
is used in the instantaneous initial boost and the remaining 
68 per cent of it is used during the sustain phase in a slightly 
“regressive burning’’ motor; i.e., the thrust drops off with 
time. The burning time is longer for the lower altitude firings. 
Fig. 6 shows velocity vs. time for each flight. For y = 90 
deg the analysis of Tsien and Evans‘ was used and it was 
gratifying to see that their results faired smoothly into the re- 
sults of the numerical calculations. 
Note that for yo = 70 deg, the velocity began to decrease 
near the end of the flight whereas it did not for the other 
cases. This is undoubtedly due to the entry into high density 
air near sea level which the other trajectories did not ex- 
perience. 


* As previously noted by others, there are some typographical 
errors in their Equations [31 and 32], p. 101; the equations (in 
their notation) should read 


as 


tt (1 — — 28 
vo? + (1 — B)m — 28 


B+ 3, 


+ 3 [31] 


= V(1 + 6)? + 48 


‘a av: 
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Fig. 6 Velocity vs. time curves for optimum thrust programming 
of ground-to-air rocket for sample problem 
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NORTH AMERICAN AVIATION, INC. & PHILLIPS PETROLEUM COMPANY 


INTRODUCE 


ASTRODYNE, INC. 


In the formation of Astrodyne, their new jointly owned 
corporation, North American Aviation and Phillips Petro- 
leum Company have brought together two broad back- 
grounds in rocket propellants and rocket hardware. 

From Phillips comes advanced research and manufactur- 
ing experience in solid propellants and JATO units. From 
North American comes demonstrated competence in missile 
systems design and the manufacture and fabrication of 
America’s largest rocket engines. 

And from these two major companies comes Astrodyne’s 
key leadership, teamed to manage America’s major missile 
projects. 

Today Astrodyne has the experience and the facilities to _ 
design, develop, and manufacture complete propulsion syss 
tems, extruded and cast propellants, rocket engines, missile 
controls and hardware, boosters, gas generator charges, and — 7 
auxiliary power units. 

Inquiries are welcomed on any phase of the solid-propel- 
lant field—from preliminary design to quantity production. | 
ASTRODYNE, INc., McGregor, Texas. 


JET PROPULSION 


f 
~ 4 


471 


a 


RAIMOND A. STRUBLE,' CHARLES 


Illinois Institute of 


The trajectory of an aerodynamically stable rocket dur- 
ing its thrust phase is obtained. The rocket or missile is 
interpreted as a point mass (with variable weight) acted 
upon by the three forces, aerodynamic drag, thrust and 
gravity. The principal limitation imposed is that the 
rocket experience positive tangential acceleration as a 
result of its thrust. In practice, the method of solution 
is best suited to those cases wherein the tangential acceler- 
ation is large, and for sufficiently large accelerations the 
results may be construed as closed form solutions. Curves 
have been prepared to facilitate the computations so that 
results of a practical application can be readily obtained. 
The method results in a simple procedure which may be 
applied to rocket and missile design problems. 


Nomenclature 
B- = dimensionless parameter defined by Equation [17] 
B’ = dimensionless parameter defined by Equation [26] 
C = dimensional drag coefficient defined by Equation [9] 
C’ = equivalent drag coefficient defined by Equation [25] 
D = drag force 
g = gravitational acceleration 
k = dimensionless parameter defined by Equation [14] 
k’ = dimensionless parameter defined by Equation [25] 
m = mass rate of discharge of exhaust gases issuing from the 
rocket nozzle 
M = instantaneous mass of the rocket 
M> = initial mass of the rocket 
M> = initial mass of propellant 
r= general parameter for the function Z, 
R = radius of curvature of the trajectory (see Fig. 1) : 
s = arc length along the trajectory - 
t = time 
T = rocket thrust _— = 
T’ = equivalent thrust defined by Equation [11 or 25] 
u = dimensionless velocity (speed) defined by Equation [16] 
u’ = dimensionless velocity (speed ) defined by Equation [26] 
uo = initial value of dimensionless velocity 
uo’ = initial value of dimensionless velocity 
v = magnitude of the velocity vector 
vy» = magnitude of the velocity vector at burnout 
Uv = initial value of the magnitude of the velocity vector 
V = dimensional parameter defined by Equation [13] 
V’ = dimensional parameter defined by Equation [25] 
w = general variable for the function Z, 
x = rectangular coordinate in the horizontal direction (see 
Fig. 1) 
Xx = integral defined by Equation [38] (see Fig. 6) 
y = rectangular coordinate in the vertical direction (see Fig. 1) 
Yx = integral defined by Equation [16] (see Table 1 and Fig. 2) 
Z, = function defined by Equation [35] (see Fig. 5) 
a = constant defined by Equation [24] 
8 = constant defined by Equation [24] 
@ = orientation of the velocity vector or trajectory measured 
from the z-axis (see Fig. 1) 
6 = orientation of the velocity vector at burnout 
9 = initial orientation of the velocity vector Bs: 


1 Introduction 


HIS paper concerns the trajectory of a rocket during its 
thrust phase. The addition of thrust into the trajectory 
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problem introduces an important force term in the equations 
of motion and also additional complications in each term in- 
volving the weight of the rocket. During the thrust phase, 
the weight will normally decrease a significant amount, due 
to the consumption of propellant, and this variation cannot, 
realistically, be neglected. 

The science of exterior ballistics of projectiles (without 
thrust) has undergone extensive development since before 
World War I. In the main, the theories interpret the pro- 
jectile as a point mass in motion under the influence of a 
gravitational field in which the major complications arise in an 
attempt to account for aerodynamic forces. The common as- 
sumption made is that all aerodynamic effects can be charac- 
terized as a drag force, tangential to the trajectory, and in 
opposition to the motion. This assumption is realistic for 
aerodynamically stable projectiles, and to a lesser extent for 
thrusted rockets and missiles in free flight. In this paper the 
rocket or missile is again interpreted as a point mass (with 
variable weight) acted upon by the three forces, aerodynamic 
drag, thrust and gravity. The drag, thrust and velocity 
vectors are assumed to be collinear so that the curvature of 
the trajectory is entirely a result of the normal component of 
gravity. The motion is therefore two-dimensional and con- 
fined to a single vertical plane. 

Recent investigations have been concerned with this rocket 
trajectory problem in limited situations. Culler and Fried 
(1)* have obtained trajectories neglecting aerodynamic drag; 
Struble and Black (2) neglecting gravitational forces; and 
Blatz (3) neglecting trajectory curvature. The first of these 
is a tabulation of trajectories obtained by numerical proce- 
dures, the other two are closed-form solutions relating the 
mass ratio to the burnt velocity. In these latter cases, the 
motion is rectilinear so that the trajectory problem reduces 
simply to considerations of the magnitude of the velocity vec- 
tor. The present analysis treats this trajectory problem 
without essential restrictions. The principal limitation im- 
posed is that the rocket experience positive, tangential ac- 
celeration as a result of its thrust. In practice, the method is 
best suited to those cases wherein the tangential acceleration 
is large. In fact, for sufficiently large accelerations, the re- 
sults may be construed as closed-form solutions. 

The basic equations are introduced in Section 2. Specific 
appligation is then made to a rocket with constant (large) 
thrust and constant drag coefficient in a uniform gravitational 
field. The extension to the general situation is described in 
Section 4. Curves have been prepared to facilitate the com- 
putations so that results of practical application can be ob- 
tained using hand computations involving only the basic 
arithmetic operations. In the general case, the method is 
tantamount to step-by-step numerical integration, in which 
the accuracy obtainable is a function of the labor performed. 
However, in any given case, the labor is not excessive and the 
method results in a simple procedure, which may be applied 


to rocket and missile design. 
= 
2 The Basic Equations 


- 
The equations of motion considered are 


dv 


M T — D — Mgsin@ 


4 Numbers in parentheses indicate References at end of paper. 
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R 
where M is the instantaneous mass of the rocket and v, T and 
D are the magnitudes of the velocity, thrust and drag vectors, 
respectively. The angle of orientation of the velocity vector 
is denoted by @ (see Fig. 1) and the radius of curvature of the 
trajectory by R. 

In Equation [1], 7 and M are generally given functions of 
time, while D varies roughly asv?. The latter term introduces 
an important nonlinearity which is avoided by considering the 
velocity (speed) v as the independent variable. This is prac- 
tical whenever v is a strictly monotonic function of time. In 
the present, analysis it is assumed that the rocket experiences 
only positive acceleration during the thrust phase so that all 
quantities may be considered as functions of v. The time ¢ is 
introduced, whenever desirable, by virtue of its known rela- 
tionship with the instantaneous mass 1/. In particular, the 
mass may be expressed in the form 


where m is the mass rate of discharge of the exhaust gases 
issuing from the rocket nozzle, and Mp) is the mass of the rocket 
at launch. It is customary (4) to interpret the thrust as 
being proportional to m, but this interpretation is not neces- 
sary for the present investigation. 

By introducing are length s along the trajectory, Equation 
[2] may be conveniently reformulated as 


ds 
g co [4] 
where ds/d@ = —R. Using Equation [3], it is found that 
dé 


so that Equations [1, 4] may be written in the differential 
forms 


dM m dv 
M 5 


Integrating these equations from the launch velocity v to the 
instantaneous velocity v, we have 


M m dv 
In — = — as [7] 
My T — D — Mgsin 
tan 
2 4 lg {dM 


In = ; 
tan 

€ 4 


where @ is the launch angle. The integral equations, [7, 8], 
are the basic equations to be solved in the present investiga- 
tion. 


3 A Special Case 


Considerations in this section will be restricted to a rocket 
in a uniform gravitational field for which the thrust is constant 
throughout burning and for which the drag is proportional to v?. 
In fact, we shall assume that the mass rate of discharge m as 
well as the thrust 7 is constant. The drag will be expressed 
in the form a 


where C is defined as the drag coefficient. Here C is a dimen- 
sional constant which depends upon the air density and the 
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Fig. 1 Coordinate system for rocket trajectory 


shape and size of the rocket. Hence, in this case, one is re- 
stricted to trajectories in which only small changes in alti- 
are experienced. 
An additional assumption is made so that the quantity Wg 
sin @ in Equation [7] may be accounted for in an approximate 
manner. Specifically, it is assumed that 7 is so large that 


T — D— Mgsin@>Mgq sin @...........[10] 


throughout burning. This condition characterizes a rocket 
with large tangential acceleration, since Equation [1] implies 
that the quantity (7 — D — Mg sin 6)/Mg is the number of 
g’s acceleration. It is ana priori assumption which depends 
upon the solution itself. However, the quantity Mg sin @ is 
the tangential component of the gravity force and hence is 
never greater than the instantaneous weight of the rocket. 
Further, if the velocity range is anticipated, the minimum 
value of T — D, i.e., the maximum value of D, can be deter- 
mined from Equation [9]. Hence, inequality [10] can, in 
general, be established without a priori knowledge of details 
of the solution. The quantity Mg sin @ is conveniently in- 
terpreted as a variation in the thrust. If, for example, it 
is approximated by its initial value Mog sin 6, the equivalent 
thrust becomes 


Equation [7] may then be integrated using [9] and results 
in 


M V+tu V—v\ 
where 
T’ m 
2V/T'C 


Of 


It is seen that V has the dimensions of a velocity, while k is a 
dimensionless parameter. Equation [12] was obtained in (2) 
for the case of a rocket in rectilinear motion, and in such a case 
completely describes the motion. Here [12] merely relates 
the magnitude of the velocity v to the mass M at each point 
along the trajectory. Since M remains positive, Equation 
[12] requires that v be less than the quantity V. 
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If d//dv is obtained from Equation [12] and substituted 
into Equation [8], there results 


(5 
In = — : Xx 
*) vo J +0 
tan 


4 
l dv 
where the only approximation introduced is in the factor 
dM/dv. However, dM/dv exhibits the same order of ac- 

curacy as .J/ itself since according to Equation [7] 


dM —mM, v m dv 


and the term J/g sin @ occurs only in the combination T — D 
— Mg sin 6. Hence, it follows that Equation [15] exhibits 
the same order of accuracy as Equation [12], when Mg sin 6 
is approximated by Mog sin %. The integral appearing on 
the right-hand side of Equation [15] is expressible in terms of 
elementary or known functions only for special values of the 
parameter k. It has been tabulated in the nondimensional 
form 


u) = exp ——— + — where u = 
0.00667 \1 + u u V 


Tabular values of this function are given in Table 1 for the 
practical range of values of k given by 0 < k < 3. Portions 
of these data have been graphed in Fig. 2 together with the 
function fork = 1. The latter may be obtained from [16] in 
closed form and provides a useful bound to the curves. 


Using [16], Equation [15] may be written 7 a? A 
= 


B Mg (1 
Y,(u) 1 — Uo 


tan (5 *) 
2 


_ Equation [12], and may be written 


Equation [17] describes the trajectory in the hodograph plane 


| oe by giving the orientation 6 of the velocity vector at each point 


along the trajectory as a function of the magnitude v (ie., u = | 
v/V) of the velocity vector. When applicable, i.e., when con- | 
— dition [10] is satisfied, Equation [17] constitutes a closed form 
solution of the trajectory in the hodograph plane. From Fig. 
2 it is seen that away from u = 1, Y;(u) is nearly proportional 
to uso that [Yx(uo)]/[Ys(u)] is nearly v/v. Hence, [17] is 
approximately 


The instant of burnout is determined from Equation [3] 
when If = M, — M,, where M, is the initial mass of the 
propellant. The velocity at burnout, v = »,, is then given by 


1 
V My 
The trajectory orientation at burnout, @ = 6,, is given by 


Equation [17], with u = v,/V or directly in terms of the mass 
ratio M,/M, by the approximate relation [18] in the form 


6, T V + Vo ( vs B 
+{1- k 
¢ 4 V-—w% Mo} 
|V 
tan 9 4 


V —w% 

After burnout, the thrust is zero and the remainder of the 

trajectory is simply that of a projectile (without thrust), 

launched at the orientation 6, and the velocity v,. ; 
‘The hodograph solution as given by [17] can be applied 
— to a wide range of practical cases. However, the accuracy of 

the solutions may be improved in many cases by a simple 

artifice. First the identity 


4) 1+sin 6 
is used to rewrite Equation [17] in the form | 
A us Y,(u) 2B 
tan"\2 + 4) ~ 
. . [22] | 


2 4 Y,(w) 


This form is particularly convenient for direct computations | 
of the quantity sin 6. Using Equations [12, 22], a plot of 
the quantity Mg sin 6 as a function of v can be made and then 
approximated by an expression of the form 


where a and Bare constants. The approximation may be crude 
but in general a and B can be chosen so that 


Mg ain 0 = — [24] 


represents an improvement over the initial estimate Wg sin 6 
= Mg sin 6. The expression [23] is, of course, chosen so 
that if [24] is used in Equation [7], the integration may again 
be performed to yield Equation [12]. The constants a and B 
are interpreted as variations in the thrust 7 and the drag co- 
efficient C, respectively. In fact, using the notation 


ir 
7 T : Cc’ Cc V’ = k’ 
a B \ 


[25] 
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the improved estimate of the quantity .V is given in Equation 28r 
{12] simply by adding the appropriate primes. 
Similarly, Equation [17] (or Equation [22]) yields an im-— 
proved hodograph solution by the addition of primes accord- 
ing to [25] and the additional notation Mog sin@,* 22.8 (Initial Approximation) —> 
v Uo 2Mog [1 + vo 
nical ex: »will trate proce £qs.(31)and (32) 
A numerical example will illustrate the procedure. 7 pa 
Example 1. We begin with the values aan 
Eqs. (28) and (29) 
My, = 1 slug; T = 300 lb; m = 0.0616 slug/sec; C = 85.5 X © o 
10 6 Ib-sec?/ft?; g = 32.2 fps?; vm = 500 fps; 0 = 45° 5 ro 
27 12 2 
= 24-5(~-) (Improved Approximation} 
= 1000 ‘oved Approximation }~ 
and compute the quantities 
7’ =T — Mogsin & = 277.2; k = 0.2; V = 1800 fps; 8 
Uy = 0.278; 2B = 0:260 
The rocket will experience approximately 8 g acceleration at_ 
launch. However, this figure will have been reduced to ap-— F 
proximately 4 g as a result of the increase in drag, when the — 
rocket reaches a velocity of 1500 fps. Using these data, an 
initial estimate of the function J/g sin 6 vs. v is computed R 
using Equations [12, 22] and Table 1. This function is il- 
lustrated by the dashed curve of Fig. 3 and is given by 300 700 900 N00 1300 1500 
1800 — v\%? VELOCITY, v, ft./sec 
M = 1.121 ( [28] 
1800 + v Fig. sin 6 vs. v, Example 1 
= 
Table 1 Yx(u) 
k 
he ee 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 
u 
0.007987 | 1.1972 1.1971 1.1969 1.1967 1.1966 1.1965 1.1963 1.1961 1.1960 1.1958 1.1957 
0.009562 1.4333 1.4329 1.4325 1.4322 1.4318 1.4313 1.4309 1.4305 1.4300 1.4296 1.4292 
0.011447 1.7160 1.7153 1.7145 1.7136 1.7127 1.7121 1.7112 1.7104 1.7095 1.7088 1.7080 
0.013705 2.0546 2.0532 2.0518 2.0503 2.0489 2.0475 2.0460 2.0446 2.0432 2.0417 2.0403 
0.016408 | 2.4598 2.4574 2.4552 2.4527 2.4503 2.4481 2.4456 2.4432 2.4410 2.4385 2.4361 
0.019644 | 2.9453 2.9414 2.9376 2.9338 2.9300 2.9262 2.9224 2.9186 2.9148 2.9113 2.9075 
0.023518 | 3.5265 3.5205 3.5145 3.5085 3.5026 3.4970 3.4910 3.4851 3.4795 3.4736 3.4767 
0.028156 4.2224 4.2131 4.2043 4.1950 4.1862 4.1774 4.1683 4.1595 4.1508 4.1421 4.1334 
0.033709 | 5.0556 5.0420 5.0284 5.0153 5.0018 4.9883 4.9749 4.9615 4.9486 4.9352 4.9224 
| 0.040357 6.0545 6.0339 6.0141 5.9936 5.9739 5.9536 5.9340 5.9145 5.8950 5.8756 5.8562 
| 0.048316 7.2507 7.2210 7.1908 7.1614 7.1321 7.1022 7.0731 7.0442 7.0160 6.9873 6.9588 
| 0.057844 8.6859 8.6408 8.5969 8.5531 8.5105 8.4672 8.4250 8.3830 8.3411 8.2995 8.2581 
| 0.069252 10.4062 10.3408 10.2769 10.2134 10.1503 10.0875 10.0262 9.9642 9.9046 9.8444 9.7845 
0.082910 12.471 12.377 12.283 12.190 12.099 12.008 11.919 11.831 11.744 11.659 11.572 
0.099261 14.953 14.816 14.679 14.544 14.413 14.282 14.154 14.027 13.903 13.780 13.658 
0.118837 17.941 17.740 17.542 17.349 17.159 16.971 16.787 16.607 16.428 16.253 16.080 
0.142274 21.546 21.255 20.968 20.689 20.416 20.146 19.884 19.625 19.373 19.125 18.884 
0.170333 25.912 25.488 25.073 24.670 24.276 23.893 23.519 23.152 22.796 22.448 22.107 
0. 186374 28.437 27.924 27.426 26.942 26.470 26.011 25.564 25.131 24.707 24.293 23.893 
0. 203926 31.224 30.606 30.006 29.424 28.861 28.312 27.780 27.262 26.760 26.274 25.801 
| 0.223130 34.312 33.566 32.842 32.143 31.466 30.812 30.178 29.565 28.968 28.392 27.835 
0.244143 37.739 36.833 35.963 35.121 34.312 33.529 32.773 32.044 31.340 30.661 30.003 
0.267135 41.554 40.455 39.405 38 . 394 37.420 36.485 35.588 34.719 33.889 33.086 32.314 
0.292293 45.815 44.483 43.211 41.989 40.821 39.706 38 . 633 37 . 607 36.624 35.677 34.768 
0.319819 50.597 48 .975 47 .432 45.957 44.558 43.220 41.942 40.723 39.559 38 .444 37.379 
| 0.349938 55.997 54.012 52.133 50.355 48 . 662 47.059 45.536 44.089 42.713 41.401 40.153 
0.382893 62.128 59.698 57.403 55.241 53.202 51.275 49.457 47.737 46.104 44.563 43.099 
0.418952 69.165 66.155 63.345 60.709 58.241 55.919 53.742 51.692 49.764 47.947 46.229 
0.458406 77.316 73.582 70.106 66.880 63.873 61.069 58.458 56.008 53.721 51.578 49.566 
0.501576 86.912 82.212 77.898 73.914 70.232 66.826 63. 669 60.740 58.021 55.484 53.122 
0.548812 98.416 92.444 87.008 82.039 77.494 73.325 69.498 65.970 62.721 59.716 56.940 
0. 600496 112.57 104.85 97.896 91.616 85.936 80.786 76.096 71.823 67.918 64.341 61.057 
0. 637628 124.08 114.78 106.50 99.087 92.435 86.453 81.053 76.165 71.729 67.694 64.008 
0.677057 137.91 126.56 116.56 107.71 99.843 92.833 86.565 80.947 75.876 71.300 67.162 
0.718924 155.04 140.91 128.63 117.88 108.47 100.17 92.823 86.289 80.463 75.249 70.556 
0.763379 177.15 159.06 143.61 130.31 118.81 108.81 100.07 92.398 85.619 79.615 74.277 
0.810584 | 207 .47 183.33 163.16 146.16 131.72 119.39 108.78 99.584 91.588 84.589 78.430 
0.860708 253.41 218.88 190.89 167.96 148.99 133.15 119.82 108.49 98.820 90.477 83.255 
0.913931 337.4 280.6 236.8 202.5 175.3 153.3 135.4 120.6 108.3 98 .00 89.23 
0.941765 | 419.6 337.5 276.9 231.3 196.1 168.7 146.8 129.2 114.8 103. 93.04 
0.970446 | 601. 454. 354. 283. 232 194. 165. 142. 124 110. 98.1 | 
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tan? + r) = 5.828 


180000 
The dotted curve of Fig. 3 is a plot of the function 


E —5 (a5) [30] 


which is chosen as an improved’ approximation to Mg sin @ 
of the form [23]. Using the latter approximation, we have 
a = 24, 8B = 5 X 10-* and obtain the new estimates 


T’=T —a = C’ = C — B = (80.5) 10-8; 
k’ = 0.207; v’ = 1850 fps; uw’ = 0.270 


With these data an improved estimate of Mg sin 6 vs. v may 
be computed according to 


1 — y\0.207 
1850 + v 
39.4 
sin 6 = ee [32] 
5.828 
+{ 39.4 ] 


where u’ = v/ 1850. 

This function is the solid curve of Fig. 3, and differs from 
the improved approximation [30] by not more than 1.3 lb 
throughout the range 500 < v < 1500. Hence, the quantity 
T — D— Mg sin @ has been estimated to within 1.3 lb, and 
since this accuracy is probably greater than one could antici- 
pate in the input data itself, there is no practical need for im- 
proving the solution. 

The initial estimates of the quantities M and sin @ are il- 
lustrated in Fig. 4 by the dashed curves while the improved 
estimates are the solid curves of Fig. 4. The latter are ac- 
curate to approximately one-half of one per cent. Also shown 
in Fig. 4 is a plot of sin @ according to [18] which in this case 
becomes 


4. The General Case 


The hodograph solutions discussed in the preceding section 
are based upon approximations to the quantity Mg sin 86. 
The latter occurs in the combination: T — D — Mg sin 6, 
and may, therefore, be accounted for in an approximate man- 
ner whenever condition [10] holds. In Example 1, we were 
able to approximate the quantity Mg sin @ to within 1.3 lb for 
the range 500 < v < 1500. Although the resulting solution 
was satisfactory, if the range were reduced, the numerical 
accuracy throughout the reduced range could have been im- 
proved. This suggests a procedure that may be applied in 
more general situations. Since Mg sin 6 is known initially, if 
its variation from the initial value is neglected, Equations 
(7, 8] can be integrated as before. The integrations then pro- 
vide a means of determining the variation itself, i.e., Equa- 
tions [12,22]. The range (in v) over which the approximation 

5 This approximation should be compared with the approxima- 
tion Mg sin 9 = Mog sin 6 = 22.8 lb. 
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500 700 900 1100 1300 1500 
VELOCITY, v, ft/sec 
Fig. 4 Sin 6 and vs. v, Example 1 


to a constant is valid can, therefore, be established. At the 
end of this range, the constant value approximating Vg sin 0 
can be changed to its current value and a second portion of the 
trajectory determined. In this way, the hodograph solution is 
constructed piecewise, using successive constant values of Mg 
sin @ computed during the procedure. In principle, arbitrary 
accuracy can be obtained, since the procedure is tantamount 
to step-by-step numerical integration. For each portion of 
the trajectory, the parameter values Mo, T’, m, C, g, vo and 
sin 6) correspond to the initial values of that particular portion 
and are not necessarily launch values. 

In this general procedure, where the trajectory is pieced 
together, variations in the parameters such as thrust, mass rate 
of discharge, drag coefficient and gravitational acceleration 
may be incorporated without additional complications. In 
each portion of the trajectory these parameters are assumed 
constant but may vary from portion to portion. Since the 
parameters may vary with time and/or altitude, in addition 
to velocity, it may be necessary to accompany the step-by- 
step procedure with calculations of elapsed time and current 
altitude. The former is a simple matter of introducing time 
through its known relationship with the instantaneous mass. 
The latter, however, requires an additional integration of the 
hodograph solution, which is desirable in any event. The 
integration of the hodograph solution is treated in the suc- 
ceeding section. 

For many ballistic missiles and rockets, significant changes 
in altitude will occur during flight. The corresponding 
variation of drag coefficient with changes in atmospheric 
density may be incorporated in the step-by-step procedure. 

The effects of discontinuous variations in the thrust and 
drag coefficient are important considerations in trajectory 
calculations. For example, a multistage rocket will exhibit 
sharp changes in thrust and, generally, drag coefficient at suc- 
cessive stages. In addition, if the rocket accelerates through 
sound speed, a sharp increase in drag coefficient will become 
apparent. Major stage changes are, in general, programmed 
as a function of time or possibly altitude, while the discontin- 
uity in drag coefficient will depend upon altitude and flight ( 
speed. Regardless of the source of the discontinuities in these 


parameters, the divisions of the trajectory may be chosen (or y 
determined, as the case may be) so as to coincide with the dis- ( 
continuities. 


Although the computational procedure is straightforward, 
some reduction in labor may be effected by the use of appro- 
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priate curves. For example, the hodograph solution given by 
[quation [17] may be written as 


(; — sin 9\1/28 (; — sin 133] 
1 + sin 6 Ad + sin & 


using the trigonometric identity [21], and Equation [12] as 


(: 
These equations are seen to contain several expressions of the 


form 


1 


which has been plotted in Fig. 5 for various values of r.. This 
compact representation of the various quantities in Equations 
|33, 34] eliminates a large portion of the computational effort. 
A ‘numerical example will illustrate the procedure. 

Example 2. Beginning with the values 


M, = 0.75 slug; 7 = 1001lb; m = 0.0346 slug/sec; 

C = 40 X 10-6 lb-sec?/ft?; g = 32.2 fps?; vy = 400 fps; 
[7 sin 0) = 0.707 

the parameters 


= 0.300; 


T’ = T — Mog sin & = 82.9; k 


V = 1440 fps; w = 0.278 
are computed. Entering Fig. 5, the value of 


is read for w = up = 0.278 andr = k = 0.3, and the constant 


in 2B 2Mog \1 + wo 
computed. It is anticipated the approximation Mg sin @ = 
Mog sin 6 = 17.1 will be satisfactory, say, for the range 
400 < v < 500. With the value of w = u = 500/1440 = 
0.347, Fig. 5 is again entered to obtain 


1+u 


k 
) =0.805 
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Fig. 5 with the quantity 
z, = 
1 — sin 0 
and — sin = w> O read along the abscissa. 
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forr = k = 0.3. The mass .V/ corresponding to this value of v _ 
is then computed by Equation [34] 


- = 0.715 
0.845 
The elapsed time is given by Equation [3] as 
M,—M 0.035 
= — == ~ sec. 
m 0.0347 


The values ¥.(u) = 37 and Y,(u) = 45 are obtained from— 
Fig. 2 corresponding to k = 0.3, up = 0.278 and u = 0.347. 
Entering Fig. 5 a third time with w = sin 6 = 0.707 and r- 
= 1/(2B) = 1.47, the quantity 


1 — sin &\!/28 0.075 
| = 0.075 
Al + sin 
‘is read and then multiplied by 


Ee 7 
> = 1.22. 
} 


to give 


=~ 


| 1 + sin % LY 
With Z, = 0.091, it is seen from Equation [33] that if Fig. 5is — 
entered along the ordinate with this value of Z,, the value of. 
sin 9 = w = 0.670 may be read along the abscissa correspond- _ 
ing to the curve r = 1/(2B) = 1.47. Hence, forv = 500, we © 
have Mg sin 6 = 15.4, which differs from the initial value Mog — 
sin % = 17.1 by 1.7 lb. This variation is approximately 2 — 
per cent of the quantity T — D — Mg sin @ and will lead to — 
about a 1 per cent error in this portion of the trajectory. In | 
some instances it is possible to encounter quantities of the 
form [35] for which Z, > 1. For example, whenever the tra- 
jectory angle 6 becomes negative, one has 


at 1 + sin 0 


In such a case, sin @ is obtained by entering the ordinate of 
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5 Integration of the Hodograph Solution 


‘The differential of are length along the trajectory may be 
expressed in the form: ds = vdt = —(v/m)dM, by virtue of 
Equation [3]. Hence, the are length from launch may be 


where dM /db is given, say, by Equation [12]. In fact, using 
[12], Equation [36] becomes 


V2B fu f1 —u\* udu 


Primes should be added to V, B, k, u and up whenever use is 
made of the quadratic approximation [24] in the hodograph 
solution. The integral 


ufi—u\* udu 
Xx (u) = (; [38] 


has been tabulated and is depicted in Fig. 6 for several 
values of k. Using [38], Equation [37] may be written 


s= [Xe(u) — Xx(uo)] ............ [39] 


It is possible to construct the trajectory by graphical tech- 
niques using [39] and the hodograph solution [17 or 22]. 
However, a numerical procedure for computing the trajectory 
directly in terms of rectangular coordinates is preferable. The 
rectangular coordinates (see Fig. 1) are related to are length 
according to the differential relations 


dx = cos 6 ds 


which in turn may be written 


dx = — cos 0dX;(u) dy =- sin 6 dX;(u). 
g g 


dy = sin 


_ by [89]. Since cos 6 and sin 6 vary quite slowly along the 


ized Closed Form for Burnt Velocity,’’ Jer PRoPpuLsion, vol. 27, 


trajectory, increments in x and y may be computed according 
to [40] using, say, average values for cos 6 and sin 6 obtained 
from [17 or 22] and increments in X;(u) obtained from Fig. 6. 
The increments in x and y are then accumulated to complete 
the trajectory integration. 


6 Concluding Remarks 


In the foregoing, methods for the rapid computation of 
rocket trajectories have been presented. The primary em- 
phasis has been on obtaining solutions in the hodograph plane. 
A step-by-step numerical procedure of general applicability, 
and a closed-form solution limited to rockets with large ac- 
celerations, have been given. Computations in the step-by- 
step procedure are reduced to basic arithmetic operations by 
the use of the curves of Figs. 2 and 5. Since these curves can 
be read to within an accuracy of approximately 1 per cent, 
the step-by-step method is limited to this accuracy. Ex- 
perience dictates that the closed-form solution with 1 per cent 
accuracy or better is available, quite generally, for rockets 
with accelerations of the order of 10 g or more. Finally, the 
trajectory in rectangular coordinates is obtained by a numeri- 
cal procedure using the curves of Fig. 6. In the general case, 
the latter computations are incorporated in the step-by-step 
procedure to allow for variations in the numerous parameters 


with altitude. 
] References 


1 Culler, Glen J., and Fried, Burton D., ‘“‘Universal Gravity 
Turn Trajectories,’’ Journal of Applied Physics, vol. 28, June 


1957, p. 672. 


2 Struble, Raimond A., and Black, Harold D., ‘‘A General- 


Feb. 1957, p. 151. 

3 Blatz, Paul H., ‘“Kinematics of a Vertical Booster,’’ Journal 
of the American Rocket Society, vol. 24, Jan.-Feb. 1954. 

4 Rosser, J. B., Newton, R. R., and Gross, G. L., ‘““Mathe- 
matical Theory of Rocket Flight,’’ McGraw-Hill, New York, 
1947. 


Technical Notes 


The Shock Tube as a Tool for Solid 
Propellant Ignition Research!’ 


> 

MARTIN SUMMERFIELD? and R. F. McALEVY UE 


Guggenheim Jet Propulsion Center, Princeton University, 
Princeton, N. J. 


; _ fundamental problem of the solid propellant igni- 
tion process is to clarify the mechanism by which the solid 
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and radiation absorption are still unknown. By means of 
a shock tube it is possible to ignite a solid propellant under 
controlled conditions to study the purely thermal process. 
This paper describes the design and operation of a shock 
tube being used for this purpose and presents some pre- 
liminary results. 


Introduction 


goup propellant combustion can be initiated by one or 
more types of excitation applied to the exposed surface. 
That is, ignition can be accomplished more or less readily de- 
pending on the circumstances, by input of heat, by surface 
reaction with certain gaseous or condensed substances, by in- 
cident radiation or by concurrent application of any combina- 
tion of these stimuli. Practical syuib igniters generally in- 
troduce all three effects. The ignition research program 
started in our laboratories in 1956 has the general objective of 
clarifying the relative roles of thermal heating, chemical reac- 
tion and radiation in the over-all process of ignition. In 
connection with this objective, it is necessary to determine 
the physical and chemical factors of a solid propellant that 
affect its ignitability. The first phase of our research is to study 


Eprror’s Note: The Technical Notes and Technical Comments sections of Jet PROPULSION are open to short manuscripts describing new 
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experimentally the role of thermal heating alone, and it is 
this aspect of the ignition problem that is the subject of the 
present paper. 

From the standpoint of fundamental combustion theory, a 
solid propellant may be regarded as a chemical substance that 
is capable of decomposing exothermally according to a par- 
ticular chemical kinetic law. The theoretical thermal ignition 
problem is solved if the time required for a specified solid pro- 
pellant to ignite can be obtained as a function of the heat ap- 
plied at a surface. The problem can be formulated one- 
dimensionally by using the unsteady Fourier heat conduction 
equation with an additional term to represent the kinetic heat 
generation 

ot Ox? 


This term is usually represented by an Arrhenius type de- 
pendence of heat release rate on temperature where Q is the 
solid-phase heat of decomposition of the propellant, R the gas 
constant, 7 the absolute temperature, A the pre-exponential 
frequency factor and £ the activation energy. A zero-order 
reaction is implied by this term. The resulting equation is a 
parabolic, second order, nonlinear, partial differential equa- 
tion relating temperature to position and time. 

The rapidity with which various propellants ignite, when 
exposed to the same heating conditions at the boundary, will 
obviously depend strongly on the chemical kinetic factors. 
The reaction rate parameters A and £ of a particular propel- 
lant fix the so-called ignition temperature which that propel- 
lant must exceed before self-sustaining combustion can take 
place. The time required for a propellant to reach this igni- 
tion temperature may be called the ignition time or the igni- 
tion time lag and is measured from the commencement of 
heating. The exact mathematical description of these phe- 
nomena requires the solution of the equation noted above. 
Unfortunately, the equation is so complicated that no 
analytical solution has been found for it. 

A numerical integration method (1)4 has been singularly 
successful in providing the solution to this problem. It was 
obtained, however, for the simplest of boundary conditions: 
A driving temperature difference described by a sharp step 
function suddenly effects the transfer of heat to the propellant 
from a surrounding hot gaseous medium that flows over the 
propellant surface at constant velocity. As a particular re- 
sult, the mathematical solution verified what had been pre- 
viously argued on physical grounds, namely, that for a given 
energy input the time to ignition decreases as the heat transfer 
rate increases. 

The results of this theoretical work have never been put to 
use for two reasons. First, the chemical reaction rate 
parameters are not well known for solid propellants and, 
second, the requirement of a sharp step function as a 
boundary condition is experimentally difficult to fulfill. 

Previous experimental work in the field of solid propellant 


‘ Numbers in parentheses indicate References at end of paper. 
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ignition consisted in transferring heat to a sample of propel- 
lant from some medium, such as a heated wire or gas, and 
measuring the time to ignition. In the hot wire work (2) a 
solution to the heat conduction equation was employed which 
made use of the ignition temperature concept to define the 
transition from pure conductive heating to chemical heat 
generation. This solution correlated the data quite well, 
since the boundary conditions were well defined, but could not 
yield the rate parameters since the kinetic term was excluded 
from the equation. Experiments employing convective heat 
transfer from a flowing hot gas to a propellant sample under 
study suffered from a lack of knowledge of actual boundary 
conditions. The experimental time to ignition when corre- 
lated with the temperature of the flowing gas yields some ap- 
parent rate parameters (3). Such results are qualitative and 
quite useful in determining the relative ease with which 
various propellants will ignite, but are too ill-defined for the 
establishment of an ignition theory or the determination of 
the chemical kinetics of the propellants. 

In designing a suitable experiment primary consideration 
must be given to a means of achieving a sharp step function 
of heat transfer rate and a heating time of the order of a few 
milliseconds. The movement of a mechanical triggering 
mechanism would have to be done in such a short time as to 
make inertial effects almost insurmountable. An easier way 
is to produce the required aerodynamic heating with a shock 
tube. The transition from zero heat transfer to some steady- 
state value after the passage of the shock takes place in a few 
microseconds. 


The Ignition Shock Tube 


The shock tube (4) is undergoing a period of popularity as 
a tool for the study of shock kinematics (reflection, interac- 
tions, etc.), aerodynamics and chemical kinetics. The chemi- 
cal kinetic studies (5) are generally centered on the reactions 
occurring in various gases subjected to shock waves. The 
present work, although dealing with chemical kinetics, will 
study chemical reactions of a solid substance heated by the 
flowing gases generated in a shock tube. 

The operation of a shock tube involves the production of a 
controlled pressure differential. The high pressure section of 
the tube is the pressure chamber, and the low pressure section 
is the expansion channel; the two sections are physically 
separated at the outset by a diaphragm. When the dia- 
phragm is ruptured a pressure discontinuity propagates into 
the low pressure channel as a shock wave which compresses 
the gas in the channel and induces a flow of the channel gas to 
follow it. If the shock propagates with constant strength, — 
then the gas following it will be at a uniform state. This high — 
temperature gas flowing at high velocity heats by convection 
any object (model) placed in the channel. The period of con- 
stant heating is terminated by the passage over the model of 
the shock reflected from the end of the tube, the expansion 
wave reflected from the high pressure end of the tube, or the 
cold expanded driving gas that was initially in the pressure 
chamber. For every combination of chamber length and _ 
channel length, there exists for the model a position of maxi- 
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mum constant heating period, which depends on the operating 
conditions (initial pressure ratio) of the tube (6). This period 
with tube length. 
ae The estimated heating rates (7) for a flow induced in a 
-__ shoek tube and the calculated heating period are theoretically 
sufficient to produce ignition of a solid propellant of reasonabie 
kinetics. A shock tube was constructed in the Flame Labora- 
tory at the James Forrestal Research Center of Princeton 
University to verify this. 

The tube was designed to produce shock waves with Mach 
numbers approaching 6 and heating times up to 7 millisec. 
Physically, it comprises a 15-ft stainless steel cylindrical pres- 
sure chamber or driver section and 56 ft of expansion channel, 
the latter consisting of five 10-ft aluminum sections of tubing 
and two 1-ft-long stainless steel sections separated by a 4-ft 
aluminum section (Fig. 1). The maximum allowable operat- 
ing pressure in the driver section is 1500 psi. The tube can 
be extended to 100 ft by inserting additional aluminum sec- 
tions. It is closed either by a cap which produces a reflected 
shock or by a dump tank which avoids the reflected shock. 

The tube length was selected to provide adequate heating 
time for most propellants, on the basis of reasonable estimates 
of the kinetic factors. The choice of the diameter was dic- 

oe  % tated by the existence of viscous effects in shock tube opera- 
a= »y ii tion characterized by deceleration of the shock wave during 
_ 7 its passage down the channel. Thus, parts of the gas at dif- 
ferent positions along the tube experience different com- 

| pressions and therefore different temperatures and induced 

P E ; velocities. As these gases flow past the model, it is then ex- 
Fig. 2 Propellant cone held in a recess of the stainless steel posed to a varying heat transfer condition that is undesirable. 
re; Gant ehhews shown ” scale This attenuation of shock strength can be minimized by in- 


° creasing the diameter of the tube (8). The tube was therefore 
Sate designed with the rather large ID of 6 in. 
The shock velocity just ahead of the model can be deter- 


mined by measuring the time required for the wave to traverse 
a known distance between two piezo-electric pressure trans- 
ducers situated in the stainless steel sections and separated by 
the 4-ft aluminum section. One of these sections containing 
a sting-mounted solid propellant cone has a viewing window to 
allow an optical system to be used to detect ignition of the 
propellant. 


STING MOUNT 
CEMENT 


20° 
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Ignition Tests 


The first phase of the present work was to demonstrate that 
the shock tube could produce ignition of solid propellants. 
Two solid propellants were chosen for test: a double-base type 
consisting mainly of nitrocellulose and nitroglycerin, and a 
composite type consisting mainly of ammonium perchlorate 
and a styrene-base fuel binder. The driver gas was helium, 
the channel gas was air. Both the composite and double base 
types were tested with the tube first in the dump tank and 
: then in the cap configurations. The cap configuration pro- 
duced heating conditions sufficient to ignite both types. 


(a) (b) (c) (d) 


_ Fig. 4 (a) Second frame—Propellant burning vigorously in a swirling reverse flow which sweeps the flame forward. (b) Eleventh 
- frame—Flame being moved about during flow breakdown. Surface regression can be seen. (c) Nineteenth frame—Propellant burning 
in quiescent gas. Pronounced surface regression. (d) Thirty-second frame—Propellant burned into recess of hollow sting with its 
bushy flame burning in the motionless surrounding gas as yea pas 
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There are indications that the double base type first ignites 
and is then extinguished by the helium flow when in the dump 
tank configuration, while the results for the composite are still 
inconclusive. (The extinguishment by cold helium deserves 
study just as well as the ignition process.) 

A typical test showing photographic proof of ignition and 
burning was made with a movie camera taking 40 frames per 
second using Kodak TRI-X exposed only by the light emitted 
by the propellant within the tube. A model of ammonium 
perchlorate composite propellant was mounted about 7 in. 
from the cap end, a distance chosen to secure the longest heat- 
ing time with no chance of quenching by cold helium. It was 
a 3 in. long, 20 deg cone mounted in the tip of a hollowed stain- 
less steel tapered sting. This is indicated in the drawing (Fig. 
2) of the mounted model as seen through the window. The 
tube operating condition was calculated to produce a shock 
wave of Mach 4.0, with a flow of air behind the shock cor- 
responding to M = 1.5, P = 65 psia, T = 1100 K. 

Fig. 3 is the first exposed frame of the series showing the 
glowing propellant cone with an incipient flame on the surface. 
Fig. 4 shows four frames of the propellant combustion process. 
Fig. 4(a) is the second frame in the series and shows the pro- 
pellant vigorously burning with a flame which is being swept 
forward of the model by the air following the shock which has 
now reflected from the capped end of the tube and propa- 
gated upstream. The next two frames are the eleventh and 
nineteenth which show the propellant cone burning and the 
surface regression. The final frame, the thirty-second, shows 
the propellant surface burning down in the hollow of the 
stainless steel sting. 

It is clear that the method is successful in producing solid 
propellant ignition. The domain of ignition for several pro- 
pellants is currently being established with various types of 
gas in the channel. Light and temperature sensitive instru- 
mentation is being developed to obtain quantitative measure- 
ments of propellant ignition times and heat transfer rates. It 
is expected that this information will be used to interpret the 
thermal ignition theories and to determine the chemical reac- 
tion rate parameters for solid propellants. 
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Trajectories: 
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a Purdue University, Lafayette, Ind. 


Introduction 


N AN interesting technical note (1)* a minimal criterion 

was supplied for the two cases of vertical flight and level 
flight. 

With the present paper, a generalization is offered in con- 
nection with arbitrarily inclined rectilinear paths. The re- 
sults derived here include, as a particular case, those of (1). 
It is stressed that the minimal properties of functional forms 
of linear type can be usefully investigated by means of 
Green’s theorem. This circumstance was already shown in 
(4), in an article apparently not noticed by the author of (1). 
In connection with level flight, Green’s theorem supplies a 
complete proof of the necessary and sufficient conditions for 
the extremum, while the Legendre-Clebsch condition alone 
is only indicative of the local behavior of the solution w ith 
respect to weak variations. ae 


Euler-Lagrange Equations 7 


A rocket-powered vehicle moving along an arbitrarily in- 
clined rectilinear trajectory is considered. Its motion along 
the tangent to the flight path is represented‘ by 


g sin 8 
J 
= mexp 
D exp 7 


where m is the mass, g the acceleration of gravity (assumed 
constant), V. the equivalent exit velocity of the rocket engine 
(assumed constant), V the flight velocity, @ the path incli- 
nation with respect to a horizontal plane and D(h, V, L) the 
drag. The symbols ¢ and y denote modified mass and mod- 
ified drag function, respectively; A is the flight altitude and L 
the lift. Notice that h = A(@) and L = L(g, V); more spe- 
cifically 


Bex 
L = gy cos exp 


After accounting for Equations [4, 5], the y-function takes 
the form y = ¥(6, ¢, V), so that Equation [1] can be re- 
garded as a first order differential equation, with dependent 
variables V and g, and independent variable @. 

Among all sets of functions V(8) and ¢(@) satisfying Equa- 
tion [1] the Mayer problem seeks that particular set® such 
that AG = G; — G; is a minimum, where G = G(8, ¢). 
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4 Primed quantities denote derivative with respect to 8, 8 bein 
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The Euler-Lagrange equations associated with the afore- 
mentioned problem are written as ieee ant 


(6) 
ag dB 


where F = )J is the so-called fundamental function and 
(8) is a variable Lagrange multiplier. In view of the fact 
that y = ¥(6, ¢, V), Equations [6] lead to 


te 
, V, L), Equation [8] can he 


After considering that D = D(h 
transformed (2) into 


DV — V.) + VV ( 
g¢ exp V. 


oD 
sin + V cos 0 =Q...{9] 


Notice that Equation [9] contains, as a particular case, pre- 
vious results obtained by Hibbs, Cicala and Miele for level 
flight (@ = 0) and by Tsien, Evans, Leitmann, Ward and 


Miele for vertical flight (@ = 7/2). 4 


The Legendre-Clebsch condition is a necessary requirement 
of a local nature to be satisfied at all points of the extremal 
are. It is to be noted that, for the problem formulated in 
the (8, ¢, V) space, the Legendre condition does not supply 
any information on the minimal or maximal character of the 
: extremal solution, if employed in the form indicated by Bliss 
in Chapter 8 of (8). 

The above difficulty, however, can be quickly offset. by 
= introducing a new variable (y) such that y’ = V and by study- 
& 


Legendre-Clebsch Condition 


ing the problem in the (8, g, y) space. The y-function as- 
sumes the form (8, ¢, y’) so that Equation [1] is rewritten as 


J=¢'+" + ¢,y') = 
= - Since no condition is prescribed on the y-variable at one of 
the two end-points, the Euler equations for ¢g and y in com- 
bination with the transversality condition yield once more 
Equations [6-9]. 

With regard to the second variation,® the 
condition leads to 


Legendre-Clebsch 


(d¢ .{11] 


by’ 
og” 
where AJ, with ¢’, y’) supplied by Equation 
[10]. After transferring the problem back into the (8, ¢, V) 
space the above condition implies that 


[12] 


where F = with V) supplied by Equation [1]. 

For the problem of maximizing the distance flown (G = 
_ —) the multiplier \ is positive everywhere. As a conse- 
- quence, the Legendre-Clebsch condition is satisfied if 
oy 2 oy 

at all points of the extremal path. The relation [13] general- 
izes the two inequalities appearing in (1) to the case of an 
_ arbitrary path inclination 6. 


6 The variations dy’ and é6¢’ must be consistent with Equation 
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V 
Fig. 1 Comparison between the Eulerian path /PF and the ar- 
bitrary control path 7QF (maximum range problem; level flight) 


Level Flight 


For the particular case of level flight, the constraining 
equation, the optimizing condition and the Legendre condi- 
tion are, 


respectively, transformed into 


Cy 
Bp [15 


where Y = W(V, ¢) in view of the constancy of the altitude. 
Analysis by Green’s Theorem 


For problems of the so-called linear type, certain char- 
acteristic difficulties associated with the use of indirect vari- 
ational procedures can be by-passed by means of a new 
technique, based on Green’s theorem. Such a technique, 
originally developed by the writer (3) in 1950 has been re- 
cently applied to problems of optimum burning program (4, 
5). 

For the sake of brevity, the maximum range problem in 
level flight is now treated under the following simplifying 
assumptions: The drag polar is parabolic with constant coeffi- 
cients; the initial coordinates V;, g; are prescribed and con- 
sistent with Equation [15]; the final coordinates V;, g; are 
prescribed and consistent with Equation [15]. 

As a first step the properties of the function OW/OV are in- 
vestigated. The analysis shows that the above function is 
negative (Fig. 1) in the region A located on the left of the 
curve OW/OV = 0 and positive in the region B located on the 
right of the curvedy/OV = 0. 

As a second step the range Af is computed along the ex- 
tremal path 7PF and along an arbitrary control path /QF. 
The following difference is formed 


dy 
(AB)ipr — (AB)ier = WV, [17] 


By applying Green’s theorem the above cyclic integral is 
transformed into a surface integral associated with the area S 
bounded by the counterclockwise circuit JQFPI. 


In view of the fact that 0 /OV < 0 at all points of the area S 
(belonging to the region A) it is concluded’ that 

7A similar proof can be developed for control trajectories /QF 
which are flown entirely in the regicn B. An analogous remark 
applies to control paths which are flown partially in the region 
A and partially in the region B.— 
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Thus, the minimal properties of the Eulerian trajectory de- 
fined by Equation [15] are rigorously and completely proved, 
while the Legendre criterion alone gives only an indication of 
a local nature about the behavior of the solution. 

In closing, the writer emphasizes that several other appli- 
cations of Green’s theorem can be found in (6) for particular 
cases of vertical flight. Furthermore, the interested reader 
san find a general theory of linear problems in (7). 


Remarks 

in the restricted variational 
formulation used in the section on Euler-Lagrange equations, 
no condition can be imposed on the non-derivated variable 
(i.e., the velocity) at the end-points. 

It should also be noticed that, at least for level flight, the 
results of the Legendre-Clebsch condition can be obtained as a 
particular case of the approach based on Green’s theorem. 

To this effect, consider a small displacement from the region 
A to the region B (Fig. 1) along a line ¢ = const; and notice 
that along this line the partial derivative 0y/OV is, in the 
immediate neighborhood of the extremal arc, a monotonically 
increasing function of the velocity. It is consequently con- 
cluded that, at points of the extremal arc, the inequality [16] 
must be satisfied. 


It should be noticed that, 
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Some Effects of Charge Configuration 
in Solid Propellant Combustion! 


LEON GREEN Jr.? 


Aerojet-General Corp., Azusa, Calif. 


Introduction 


T HAS long been known that the stability of the reaction of 
a lateral-burning, solid propellant charge can be pro- 
foundly influenced by the shape of the flow channel or (in the 


1958 


case of internal-burning charges) the charge cavity. It has 
also been shown that unstable reaction is associated with the 
occurrence of transverse acoustical oscillations within the 
charge cavity, a phenomenon which may be called ‘“‘sonant” 
burning (1)* and which appears to be a condition necessary 
but not sufficient to produce “resonant” burning, or the 
occurrence of gross irregularities or ‘secondary peaks” in the 
mean chamber pressure vs. time history (2). In charges of 
small size, a cavity of irregular cross section (too irregular, 
presumably, to be compatible with a highly conservative 
mode of transverse wave motion) is sometimes sufficient to 


- suppress sonant and resonant burning (3), but this beneficial 
influence is reduced as the scale of the charge is increased. As 


for internal-burning charges with cavities of simple, sym- 
metrical cross section, the general development experience at 
this laboratory has been that charges with annular cavities 
are more difficult to “stabilize” than those with cylindrical 
‘avities (4). 

Most theoretical studies of unstable combustion of solid 
propellants have been limited to the sonance problem, i.e., 
have been concerned with the conditions under which the 
amplitude of a weak acoustical oscillation can grow to a finite 
value. A recent analysis by Cheng (5) indicated that the con- 
ditions for growth of a small pressure perturbation are more 
favorable in the case of a cylindrical cavity than in an annular 
cavity, a trend in apparent contradiction to the static test 
experience mentioned above, and it has been suggested that 
this discrepancy may be at least partly attributable to a lack 
of distinction between the sonance and resonance phenomena 
(2). Another possible source of disagreement is that the de- 
velopment experience was derived from tests of conventional 
propellant charges, with cavities of relatively high length-to- 
diameter ratio. Although the prime excitation of sonance is 
provided by the gas-phase combustion reaction (1), it has 
been suggested that an appreciable contribution may be af- 
forded by secondary viscosity effects (sound generation) in 
the subsonic jet leaving the aft end of the charge (4). This 
suggestion was based on the experience of numerous investiga- 
tors that resonant burning became more pronounced as the 
aft-end gas velocity increases (with increasing length-diameter 
or throat-port ratio), as well as on the observation that the 
configuration of the flow channel at the nozzle entrance (i.e., 
the shape and erosion resistance of the charge and restriction) 
exerts a critical effect upon reaction stability in cases where the 
gas velocity at this point is fairly high (1). This observed 
dependence of resonance severity upon gas velocity is qualita- 
tively consistent with a dimensional analysis of aerodynamic 
sound production which shows that the intensity of sound 
radiated by a subsonic jet is roughly proportiona! to the eighth 
power of the stream velocity and the square of the characteris- 
tic length involved, and that the efficiency of sound produc- 
tion is proportional to the fifth power of the Mach number (6). 

In order to obtain a more valid test of the Cheng theory, 
which does not consider viscosity effects, and also is restricted 
to charge cavities of a small length-to-diameter ratio (i.e., less 
than about 15), it was desired to compare the behavior of 
short, internal-burning charges using annular and cylindrical 
‘avities, employing the same propellant. In addition, tests of 
external-burning cylindrical charges using high frequency in- 
strumentation were considered desirable in order to determine 
the mode of any pressure oscillations prevailing, since ex- 
periments by Huffington (7) have indicated that the normal 
burning of cordite rods can show a fine periodic character, 
even though no instances of grossly irregular reaction of this 
type of charge are known. Finally, the reaction of a star-per- 
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RESTRICTION 


ROD SHELL 


PROPELLANT 0.D. 4 
PROPELLANT 1.0. 3.00 


Fig. 1 


forated charge was thought worthy of investigaton with high 
frequency response instrumentation, since the acoustics of 
irregular cavities are too complex to permit theoretical study. 


Equipment and Procedure 


For the purpose of this program, charges of the type de- 
scribed above were fabricated of a high energy composite 
propellant with ballistic properties given in Table 1. The 


Table 1 Some ballistic properties of propellant employed 


in the static testing 


Strand burning rate at 1000 psi and 60 F, in./sec.. 0.34 
Theoretical specific impulse at 1000 psi, Ibf-sec/Ibm. 239 

Characteristic gas velocity, c*, fps............... 4800 


propellant grains employed in this study were fabricated by 
casting a slurry of ammonium perchlorate oxidizer particles 
and thermosetting, resinous fuel in molds, followed by curing 
at an elevated temperature. The dimensions of the charges 
are shown in Fig. 1. Most of the experimental information 
obtained in the present study consisted of pressure measure- 
ments recorded by two parallel instrumentation systems of 
different frequency response. A block diagram of this dual 
system is shown in Fig. 2. The details of the high frequency 
response system and its limitations are described in (8). For 
the purpose of the present discussion, the pressure-measuring 
technique may be described as follows: The output of the 
variable-capacitance pressure pickup is recorded on magnetic 
tape and by a galvanometer oscillograph, and is also dis- 
played upon a cathode-ray oscilloscope viewed by a recording 
camera. In addition, the output of a variable-reluctance 
pressure pickup connected to the pressure tap in the nozzle 
entrance section is recorded by the galvanometer oscillograph. 
This record, of relatively low frequency response, provides 
the basis of the standard firing curves (pressure vs. time) 
which are then interpreted in terms of the high frequency 
response data. These latter data are available in two forms: 
First, as a frequency spectrogram of the observed pressure 
history, which is obtained by displaying the magnetic tape 
484 
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record (using repeated play-through) upon a frequency 
analyzer which reduces any complex oscillation into its Fourier 
components. This technique provides quantitative informa- 
tion on the frequencies prevailing during the test, as well as a 
semiquantitative measure of the amplitude of the various fre- 
quency components revealed. Second, the film record of the 
oscilloscope display is examined visually. This record yields 
quantitative amplitude data, but only of the total, syn- 
thesized pressure oscillation, and not necessarily of its in- 
dividual frequency components. 

The high frequency response pressure pickups employed in 
these tests were variable-capacitance units of an Aerojet- 
General design which employ a water-cooled shield as protec- 
tion from the hot combustion gases. The shielded pickups 
were mounted in the orifices at the fore-end of the chamber so 
as to be approximately flush with the inner surface of the for- 
ward head. Two pickup locations, spaced 120 deg apart, were 
provided. In a few tests where measurements of phase rela- 
tionships in the recorded pressures were attempted (with in- 
conclusive results), both locations were employed. In most of 
the tests, however, only one pickup was used, and the second 
orifice was plugged. 

All the tests comprising the present investigation were per- 
formed at 170 F, a temperature which had been found 
favorable for the occurrence of irregular reaction of the pro- 
pellant employed. The average chamber pressures (as de- 
termined by the ratio of burning surface area to nozzle throat 
area) ranged from about 500 to about 2000 psi, with the 
majority of the tests being conducted at nominal pressures in 
the vicinity of 1200 psi. The igniter consisted of a 50-gm 
charge of coarse (Grade FFG) black powder, loosely contained 
in a Koroseal bag, and initiated by a DuPont S-67 electric 
squib. 


Experimental Results 


The details of the experiments and their results have been 
reported in (9) and may be summarized as follows: As ex- 
pected, the external-burning rod grain yielded the most stable 
reaction; its three firing curves were uniformly smooth, and 
the frequency spectrogram data indicated the occurrence of 
only very low amplitude oscillations, almost exclusively in the 
longitudinal modes. The tubular charge ranked second in re- 
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action stability from both the gross and fine-time-scale view- 
points; only two of the seven tests using this configuration meacemeut 
— showed secondary pressure peaks (as distinguished from the 
; initial ignition peak), and only in the case of the most pro- 
nounced irregularity were oscillations (longitudinal) of ampli- vanits rare 
tude greater than a few psi detected. Irregular reaction was tw Re 
encountered frequently with the rod-in-shell charges; in only 
one of nine tests was a smooth firing curve of the normal shape - 
obtained. In one irregular test the film record revealed os- 
cillations of amplitude greater than 200 psi and with fre- 
. quency-vs.-time characteristics apparently corresponding to 
the first radial mode, as indicated by comparison with the apidimaes 
acoustical calculations of (10). 
Contrary to expectations, the star-perforated charge fur- 
nished the most exaggerated examples of resonant burning; of 
a in only one of six tests was a normal firing curve obtained, and - 
: the mean pressure irregularities observed were more pro- 
ae nounced than those yielded by the charges with cylindrical or 
7 annular cavities. In two of these tests, moreover, oscillations sale @ 
_ of high amplitude (up to 375 psi) were detected at frequencies 
corresponding to a high-order transverse mode. In the case of 
the complex, star-shaped cavity, no theoretical information eh 
concerning natural transverse modes of oscillation is avail- 
able, but photographic viewing experiments (using a chamber 
with a transparent fore-end closure) suggested one type of ah O 
transverse mode with a natural frequency corresponding 
closely to the frequency observed by the two pickups. This fe) 
mode is envisaged as a tangential motion (depicted in Fig. 3 
simply as a traveling front) in each of the star points, once 
they have become broadened by partial consumption of the 
ency grain. The calculated frequency of this motion is presented in | és 
irier Fig. 3 as a function of cavity width, together with the re- a Oo : * 
rma- corded frequencies measured at various time intervals of the = (ome) a 
as a two tests corresponding to the same range of cavity dimension. te 
fre- A comparison of the mean burning rates of the propellant c ws 
the in the various charge configurations considered is made in re 
elds Fig. 4. It may be seen that the propellant showed sig- = 1 
syn- nificantly different burning characteristics in each configura- nic a 
| in- tion. The internal-burning, star-perforated and tubular 
grains exhibited the same pressure exponent of burning rate, CALCULATED FREQUENCY OF - 
d in but, at a given pressure, the star-perforated grains burned 
yjet- more rapidly. On the basis of only three closely spaced datum 7 
tec- points, the external-burning rod grains appeared to show a 
ups higher pressure exponent, and the rod-in-shell charges showed ie i Oo 
r so an intermediate value. No satisfactory explanation for this 7 
for- behavior is apparent at present. 
in- Fig. 3 Frequency of high-amplitude oscillations in star-shaped 
t of 1 Smith, R. P., and Spencer, D. F., ‘Combustion Instability cavities vs. width of starpoint a 
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posium (International) on Combustion, Williams ard Wilkins, 10 : 
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2 Green, L., Jr., “Observations on the Irregular Reaction of | 
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oat Rockets,’’ McGraw-Hill, New York, 1950, p. 129. = 60 
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I. General Theory,’’ Proceedings of the Royal Society of London, | 
Series A, vol. 211, 1952, p. 564. ” ; 
7 Huffington, J. D., Unsteady Burning of Cordite,’’ ROD 
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_ Simplified Equations for Transient Heat- 
Conduction to Insulated Metal Slab 


JAMES H. WIEGAND! 


— 


_ Aerojet-General Corp., Sacramento, Calif. 


The transient heat-conduction to an insulated slab of 
infinite thermal conductivity can be described by Y 
(T, — T)/(T, — T;) = a exp (—6bX) for values of X> 0.5 and 
~Y< 0.9. The values of a and b have been empirically evalu- 
ated from the graphical solution of Grover and Holter as 


a =1+40.35/(1 + m)(1 + n); 1/b = (1+ m)(1 +n) — 0.60. 


Nomenclature 


intercept on log Y vs. X graph 
numerical slope on log Y vs. X graph 
heat capacity of insulation, Btu/hr °R 
heat capacity of metal slab, Btu/hr °R 

= heat transfer coefficient, Btu/hr °F ft? 

thermal conductivity of insulation, Btu/hr ft °R 
thickness of insulation, ft 
thickness of metal slab, ft 

k/hl = reciprocal of Biot number 
PmC ml m/pCl 

temperature of metal-insulation interface at time 0, °R 
temperature of hot gas, °R 

ao /|? 

unsatisfied temperature ratio (T, — T)/(T, — Ti) 
k/pC = thermal diffusivity of insulation 

time from start of heating, hr 
density of insulation, lbm/ft* 
density of metal, lbm/ft® 


3 


s 


> 


Introduction 


HE transient heat conduction in a multiple wall has been 

treated by Jaeger (5),? who pointed out the great complex- 
ity of the problem. He suggested that such a wall might be 
treated as an equivalent homogeneous wall having the same 
transient behavior at a specific point as the composite wall. 
He gave several solutions for such equivalent walls and 
showed that the analytical relations between the equivalent 
wall and the composite wall were different for different heating 
conditions. Unfortunately, he did not give the equivalent 
wall solution for the composite wall heated by a hot gas of 
constant temperature through a constant heat transfer co- 
efficient, the case of greatest interest in rocket work. 

The transient equations for a composite slab have been pre- 
sented by Mayer (6), and the numerical solution of the case 
of the second slab having infinite thermal conductivity was 
presented by Grover and Holter (3). Dimensionless plots 
similar to those of Gurney and Lurie (4) were presented by 
Grover and Holter, which correspond to the mid-plane lines 
on a graph for a homogeneous slab. The lines on such a graph 
are straight except for small values of time, and in homo- 
geneous solids they represent the solution when all terms ex- 
cept the first term have become negligible. The use of these 
straight lines was pointed out by Chung and Jackson (1) for 
a solid cylinder, and by Geckler (2) for a hollow cylinder. 
The constants in the simplified relations for a cylinder, a 
sphere and an infinitely wide plate were given by Wiegand (7). 

The temperatures at any point in a composite slab change 
in this same simple way with time, easily expressed as 

Y = (T, — T)/(T, — T;) = a exp (—bX) 
with the slope 6 the same at all points in the slab, and the in- 
tercept a varying with position. The possibility of using 
empirical relations for b and a has been explored and the pre- 
liminary results are shown in this paper. 
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Empirical Evaluation 


The tabulation of the constants for the mid-plane of an in- 
finitely wide slab given in my earlier paper (7) were used to 
derive the empirical equations. 


a= 1+ 0.29(1 + 
1/6 = (1 + m) — 0:60................ 8B] 


These do not, of course, predict the values derived from 
numerical calculation precisely, the value of 6 being 0.1 per 
cent low at m = 10, 1.5 per cent low at m = 2, and 1.5 per cent 
high at m = 0, and the error of a being less than 0.1 per cent 
for m > 0.05, and a maximum of 1.3 per cent at m = 0. 

In considering the composite slab graphs of Grover and Hol- 
ter (3), the graphs available in the publication were used, so 
the slopes and intercepts could not be read with the precision 
of the numerical tabulations for deriving Equations [2, 3] 
above. In an effort to find a simple means of correlation, the 
terms 1 + mand n were considered. It will be noted that 


and 


(m+n+mn) [4] 


where the uw is that used by Grover and Holter. Considering 
the value of dimensionless time X, we can write 


X = aO/l? = 0/(l/k) (IpC) 


If one were to write the value of X for an equivalent homo- 
geneous wall composed of a metal and an insulation, one 
could write 


Xs = = 


and if we visualize resistance and capacitance in usual terms 


Equation [8] can be written in terms of the ratio n 
(9] 
Recalling that k,, is infinite, Equation [6] becomes 
= + n) = + n)...... [10] 


Arguing then from the type of relations observed for the 
slope and exponent for a simple slab, and the line of reasoning 
on a composite wall, it appeared reasonable that terms 1 + m 
and 1 + n might correlate these composite wall data. 

The values of a and 6 were calculated from Figs. 1 to 5 of 
Grover and Holter (3), the intersection of the straight lines 
with X = 0 giving a = Yo directly. The slope of the straight 
lines is equal to —6/2.303, the figures being plotted to the 
logarithmic base 10. Table 1 presents the values of a and 6 
observed. When 1/6(1 + 7) was plotted vs. m on arithmetic 
graph paper, parallel lines were obtained of slope 0.96, whose 
intercepts in turn gave a straight line against 1/(1 + n). The 
resulting equation is 


1/b = (1 + n) (0.96 + m) — 0.67.......... 


If a slope of one were used for the lines, somewhat greater 
scatter of the intercepts was observed to give 


1/b = (1+ n)(1 + m) — 0.60........... 


The data for intercepts were much more variable than the 
slopes, due in part to the errors in reading from small figures. 
The values of a gave straight lines through the intercept on a 
graph of a vs. 1/(1 + n), with the slopes of these lines giving 
a reasonable correlation against 1/(1 + m), to give the rela- 
tion 


a=1+4 0.35/(1 + + n) 


The relation was verified by plotting a vs. 1/(1 + m) first and 
then the slopes vs. 1/(1 + n). 
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Table 1 Observed and calculated slopes and intercepts oe 

m bobs tous Eq. 1] Eq. [12] 
0 1.33 wal 2.56 3.45 2.50 
1.22 1.30 1.11 0.3 

= 1.10 O74 0.80 0.71 0.3 

1.10 0.48 0.45 0.42 0.3 

1.04 0.161 0.165 0.156 0.3 

0.5 «1.24 1.27 0.25 
0.64 0.65 0.61 0.3 

1.12 0.44 0.44 0.42 0.4 

1.08 0.265 0.270 0.256 0.5 

1.03 0.102 0.105 0.101 0.5 

1 1.14 0.72 0.78 0.71 0.3 
1.14 (0.44 0.44 0.42 0.3 
1.10 60.31 0.31 0.29 0.3 

1.06 0.190 ().192 0.185 0.3 

1.02 0.073 0.077 0.075 0.3 

2 1.07 0.44 0.44 0.42 0.3 
05 0.265 0.266 0.256 0.4 

1.04 0.189 0.190 0.185 0.4 

1.03 0.120. 0.122 0.119 0.4 

1.02 0.049 0.050 0.049 0.5 

6 1.06 0.162 0.159 0.156 0.4 
1.05 0.103 0.102 0.101 0.4 

1.04 0.074 0.075 0.075 0.5 

1.03 0.050 0.050 0.049 0.5 

1.01 0.0213 (0.0208 0.0207 0.5 


Values of a and 6 calculated from Equations [11, 12, 13] are 
given also in Table 1. It can be seen that Equation [13] 
vives good agreement with the values observed, and Equation 
{12] gives a smaller average deviation from the observed than 
does Equation [11]. 

Inspection of the Figs. 1 to 5 of Grover and Holter (3) in- 
dicated the values of X above which a straight line applied and 
the value of Y below which it applied. These values are also 
tabulated in Table 1 
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On the Application of Van Driest’s 
Method to a Highly Cooled Partially 
Dissociated ‘Tesbalem Boundary Layer' 


HENRY HIDALGO? 


AVCO Research Laboratory, Everett, Mass. 


The purpose of this note is to present an extension of the 


Van Driest method for the calculation of the shear stress : 


Juty 1958 


Industrial and Engineering © 


JET PROPULSION, vol. 25, 1955, pp. — 
I 


Solution of the Transient — 
Heat Conduction Equation for an Insulated, Infinite Metal Slab,”’ z 
JET PROPULSION, vol. 27, 1957, pp. 1249-1252. _ 


in a compressible turbulent boundary layer to include 
cases where the molecules in the boundary layer on a blunt 


body are partially dissociated. 


hese mean continuity, momentum and energy equations 
for a highly cooled, partially dissociated boundary layer 
on bodies of revolution follow from the general equations in 
(1)%as4 


ar) + 0 [1] 
or P oy 
_ Ou + ou Op (2) 
= — + ( 2 
__ OH OH H (3) 
i ( = 
pi + pi € 5 


it should be noted that the total enthalpy H includes the 
enthalpy of dissociation in addition to the perfect gas static 
enthalpy, such that H = hides + Auissociation + (U?/2). 


_ Equations [1—3] are valid under the following assumptions: 


- (a) The turbulent and laminar Lewis numbers are unity. 
~ (b) The turbulent and laminar Prandtl numbers are unity. 

(c) The density fluctuation in the boundary layer is 

As shown in (1), provided assumptions (a) and (b) are 
satisfied, then even without assumption (c) Equations [2 and 
3] show that the familiar relation between enthalpy and ve- 
locity (the Crocco integral) still exists in the case of a disso- 
ciated boundary layer with zero pressure gradient for either 
an insulated body or a constant surface enthalpy. Thus the 
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_ Fig. 1 Local turbulent skin friction coefficient for a highly 
cooled, partially dissociated turbulent boundary layer 


NORMALIZED SKIN FRICTION COEFFICIENT, 


static enthalpy, including the enthalpy of dissociation, is 
given by® 


h, hy be hy 2h, \uq 

~ where we have now dropped the bars. 
Furthermore if, as is also shown in (1), a 4th power velocity 
distribution is used in Equation [4], then the change in the 
flow properties through the boundary layer in thermodyna- 
~ mic equilibrium on a blunt body is found to be concentrated 
- in a region next to the wall. With the exception of a thin 
_ region near the wall, the density can be represented closely 

the expression 


p 

Po 

p 


5 The subscript © refers to conditions external to the boundary 
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When Equation [5] is used in the momentum integral equa- 
tion for the calculation of the wall shear stress, the value of 
the integrand near the wall will be in error. However, this 
error can be shown to be small because the momentum con- 
tribution in this region is negligible. When the calculation is 
carried out in a manner similar to the method of (2), the ex- 
pression obtained for the wall shear is the same as the cor- 
responding Van Driest formula for the case of zero dissoci- 
ation, with the exception that all temperatures are replaced 
by the corresponding enthalpies. More specifically, the 
formula for a highly cooled, partially dissociated boundary 


layer is a 
0.242 
Th (sin~' + sin~' B) = 0.41 + 
Mw 
where 
{2 = ?/2h B _ 
— B 
B 


(B? + 4A%)'/2 ~ (B2 + 4A42)'/2 


The results of (3) indicate that the value of D depends on 
whether the Prandtl mixing length or von Karman similarity 
is used in the analysis. The result is then 


D 
D 


(hy/h,,)~'/* (Prandtl mixing length) 
1 (von Kérmén similarity) 


It should be mentioned that the first term in the right-hand 
side of Equation [7] is a weak function of C, as defined in 
Equation [6]. As the value of C is increased from 1.0 to 
2.5, this term only increases from 0.41 to 0.45. It is also im- 
portant to mention that the approximation of low cooling 
made in the analysis (2) is not necessary in order to arrive 
at the corresponding Equation [7] for zero dissociation. 

The results of Equation [7] have been applied to the case of 
a hemisphere-cylinder body in hypersonic flight, by assuming 
that the two-dimensional formula is valid on the cylindrical 
section of the body where there is no pressure gradient. Be- 
cause of the entropy increase across the bow wave in front 
of the body, the Mach number of the inviscid flow over this 
cylindrical section is in the low and moderate supersonic 
range. In fact, it can be shown that for this case (u.,2/2)/H. 
< 1/2; therefore the parameter H/ ,,/h,, is taken to vary in the 
range 1.2 < H,/h, < 2.0. Ifa constant specific heat ratio 
y = 1.2 is used, this range would be equivalent to a free 
stream Mach number variation of 1.4 <M, < 3.2. Fig. | 
shows the ratio of c;../c,; calculated from Equation [7] as a 
function of Reynolds number (R,,) and H,/h,., for two values 
of Tw/T.,.. The compressible skin friction coefficient is based 
on the von Karman similarity and is normalized by using the 
incompressible skin friction coefficient c,;, which is propor- 
tional toR,~'”*. 

Fig. 1 indicates for the case under consideration a decrease 
in the skin friction coefficient of about 20 per cent when 1.2 
<H < 20. 

It should be mentioned also that Equation [7] allows for 
the calculation of the heat transfer on the cylinder through the 
use of Reynolds analogy (1). 
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Propellant Explosives Classification and 
the Effect on Field Handling of Missiles 


W. F. HAITE! 


Thiokol Chemical Corp., Huntsville, Ala. 


RESENT solid propellants fall into two main categories: 
Those in which the fuel itself contains sufficient chemi- 
cally bound oxygen for combustion, and those comprising an 
inert fuel and crystalline inorganic oxidizers, principally ni- 
trates and perchlorates. The latter, formulated from any of 
« number of inert fuels, fall into a general group known as 
composite propellants. Composite propellants utilizing am- 
monium perchlorate as oxidizer have the widest use in missile 
application; therefore, they are of particular concern today. 

Manufacture, storage and field use of solid propellant en- 
gines are controlled by the military, and are necessarily ac- 
complished under rigid safety rules. Composite propellants 
are higher energy materials than standard explosives and 
single and double base propellants which are generally known 
to be sensitive and readily detonated. Limited experience 
with respect to sensitivity caused composite propellants to be 
placed in the same general group and to be subject to similar 
hazards classification. 

Before proceeding, it is desirable to define one of the terms 
and some of the common classifications frequently referred 
to in discussions of propellant utilization: 

Quantity-Distance. The distance from a certain location, 
such as an inhabited building, public railway, magazine or 
operating building, required to protect that location against 
substantial structural damage from the ignition or explosion 
of a definite quantity of a specific class of explosives as defined 
below. 

Class IJ. Fire in these materials produces intense heat 
which is dangerous to personnel and equipment in the vicinity. 

Class IX. These materials can be expected to detonate 
when involved in a fire, and are subject to mass detonation by 
a detonating initiator. 

Principal reasons for the effectiveness of solid propellants 
are that they contain sufficient oxygen for combustion of the 
fuel and they can be ignited readily with a simple system. 
These inherent factors, although desirable, along with the 
fact that experience with propellants was very limited less 
than a decade ago, led to the establishment of these materials 
as Class LX explosives. Because of this high hazard classifica- 
tion, construction of facilities was quite costly, and produc- 
tion of any reasonable quantities necessitated the use of large 
acreage. 

Since the late 1940’s, continuing effort has been expended 
to obtain proper safety classifications for cured and uncured 
composite propellants. The drop test and friction test 
normally employed to catalog explosives were found to be 
inadequate, since the sensitivity of oxidizer-fuel combinations 
as measured in contact with high explosive charges was 
directly proportional to the oxidizer concentration; whereas 
the standard impact tests generally gave results which showed 
the inverse to be true. 

When mixed with a common oxidizer, different binders 
produced propellants with different sensitivities at the same 
oxidizer loading, and it was soon shown that per cent by 
weight of available oxygen was not the proper criterion for 
classification. 

If we look at these binders, we find that their densities 
differ, and it became readily apparent that sensitivity must be 
associated with amount of binder separating the oxidizer 
crystals. 

Presented at the ASME-ARS Aviation Conference, Dallas, 
Texas, March 17-20, 1958. 

1 Chief, Pocess Development Department. 
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These observations led to the tentative establishment of a 
relationship between the volume ratio of oxidizer to fuel and 
the detonating potential of the propellants. 

A definition of safety hazard associated with the production 
and use of polysulfide-perchlorate propellant based on the 
volume ratio of oxidizer to fuel was evolved from extensive 
tests of compositions comprising 57 to 88 per cent ammonium 
perchlorate, by weight, in quantities up to 2300 lb. Samples 
were subjected to heat, low energy electric sparks and impact, 
and initiated by electric squibs and explosive boosters. 

Negative results with large and small samples correlated 
and demonstrated that a mass effect, usually associated with 
explosives, is absent. The correlation between large and 
small samples is of particular importance because it obviates 
the testing of large masses to prove explosive characteristics 
for large missile applications. Tests to determine the hazards 
classification included: 

1 One pound samples of compositions similar to those 
which might exist during any phase of processing were tested 
by initiation with explosive boosters. Boosters comprised a 
No. 8 electric detonator and the detonator plus a 5-gram tetryl 
booster. Detonation was defined by noting spalling of con- 
crete or disintegration of cement blocks on which the samples 
were placed. 

2 Propellant containing 85 per cent perchlorate was con- 
fined in 12-in. sections of 4-in. pipe and initiated with 500 
grams of composition C-4 and a No. 8 electric detonator. 
Only partial detonation occurred in the section on which the 
explosive was placed, and there was no propagation to identi- 
cal samples butted up against the end of the first. 

3 Propellant containing 74 per cent perchlorate, confined 
in 4-ft and 1-ft sections of 8-in. pipe, was initiated as in item 2 
above, but with 2500 grams of composition C-4. Again, only 
partial detonation occurred in the 4-ft sections initiated on 
one end by the explosive, and there was no propagation to 
either of the two 1-ft sections butted in line to the other end. 

4 Attempts to cause transition from deflagration to 
detonation in large masses of uncured propellant were un- 
successful. Quantities up to 500 lb of standard propellants 
were initiated by squibs buried in the mass, and, as ex- 
pected, only rapid burning occurred. Negative results were 
obtained also with 50 grams of tetryl and a No. 8 detonator. 

5 Prototype engines containing approximately 2400 lb of 
propellant were destroyed by setting off the warhead. In 
these tests, the propellant failed to detonate; in fact, the in- 
tensity of the explosion was actually lessened due to the 
damping action of the propellant. 

Data typifying the tests with 1-lb samples are presented 
graphically in Fig. 1. The two distinct groupings are ration- 
alized by the difference in densities of the binders. It is 
obvious that a lower volume ratio of oxidizer-to-fuel will be 
obtained with a lower density binder. When composite pro- 
pellants with volume ratios in the nondetonating range are in 
the form of a cohesive mass, either in the cured or uncured 
state, the greatest hazard potential was found to be fire. The 
rubbery binder attenuates the detonation wave from the 
boosters and inhibits propagation of the wave. 

The observed results of the tests complement a theory on 
compression wave propagation based on the great difference 
in the moduli of the oxidizer and binder.? By application of 
acoustical theory, it was shown to be virtually impossible to 
transmit a wave for any significant distance through a hetero- 
geneous material composed of a low-modulus resin binder and 
a relatively high-modulus inorganic salt. 

There are a number of basic binders which can be grouped 
according to chemical structure and density. So long as they 
are inert and produce low-modulus materials, the hazard po- 
tential is independent of their chemical structure. The ap- 


2 Ritchey, H. W., ‘Explanation of Non-Detonating Character 
of Polysulfide-Ammonium Perchlorate Propellants Based on 
Acoustic Wave Transmission Velocity,” Twelfth Army-Navy-Air 
Force Solid Propellant Meeting (Confidential). 
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proximate range of density is from 1.0 to 1.3 grams/ce, and it 
is the density that is of primary consideration in applying the 
oxidizer-to-fuel volume ratio. 

Basically, when polymerized, these binders have rubbery 
properties which can be altered by varying the cure mecha- 
nism. We will consider that in all cases the propellants made 
from these binders have sufficient elongation and strength to 
withstand the temperature conditioning imposed and other 
requirements placed on the rocket engine. Because of the 
difference in density of the binders, a safe volume ratio of 
oxidizer to fuel can be maintained in spite of higher oxidizer 
contents. 

Generally, those propellants containing a higher per cent 
oxidizer have proved to be higher performance propellants. 
As shown in Fig. 1, the use of lower density binders permits 
higher performance without increasing the sensitivity to 
detonating initiators. 

As a result of the extensive tests which demonstrated the 
apparent insensitiveness of the compositions tested, the Chief 
of Ordnance Safety has established polysulfide-perchlorate 
propellants containing 74 per cent ammonium perchlorate or 
less as Class II materials. A formal recommendation, with 
high level endorsements, has been made to the Chief of Ord- 
nance Safety to revise the Ordnance Safety Manual to provide 
a proper definition of safety hazards associated with com- 
posite propellants and realistic safety classifications. The 
proposed revision is based on the conservative volume ratio 
of 2.85 (Fig. 1). The principal provisions of the proposed re- 
vions are that: 

1 Composite propellants made of an essentially inert fuel 
and in which the volume ratio of oxidizer to fuel is 2.85 or less 


shall be regarded as Class II material. 


2 When the explosive in the warhead is separated from the 


~ Class II propellant by more than 6 in., quantity-distance shall 


be assessed on the weight of the warhead alone. 

In the construction and utilization of buildings for the pro- 
duction and use of propulsion units for missiles, rigid restric- 
tions have been placed on the quantity of explosives that may 
be located in a single area or areas adjacent to administrative 
or operating buildings, which depend on the distances from 


1 Launching battery—three launchers, two missiles /launcher 

2 Control room—barricaded, Class IX propellant 

3 Control room—barricaded, Class IT propellant 7 

4 Ready storage—30 missiles Ca 
Perimeter—Class IX propellant 

— Perimeter—Class II propellant 


Fig. 3 Typical missile launching site land area 
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and it adjacent buildings, public property or operating areas. 
ig the These restrictions have been given definite numerical values 3. It was assumed that the missile contained 6000 Ib of pro- 
from which quantity-distance tables have been established. pellant and carried a 1000-lb high explosive warhead. It is 
bbery The data from these tables are shown graphically in Fig. 2, evident that the land area is greatly reduced when the engines 
\echa- intraline distance being that distance between operating build- contain Class II propellant. The required areas, based on 
made ings and inhabited building distance that distance between quantity-distance assessments for Class LX propellants is 1200 
sth to operating buildings and administrative buildings. As a acres; a reduction of 68 per cent of this area is realized with 
other general rule, the inhabited building distances apply to public Class II propellant. 
of the highways and railroads. The graphical presentation of the Using a very conservative figure of $300 per acre, $360,000 
tio of data shows that a considerable saving in area alone could be is saved on each site. Added costs of providing services, 
idizer realized by utilizing a propellant with a hazards classification transportation and facilities could easily increase the figure to 
of Class IT. one-half million dollars. When the over-all economics of 
- cent Missile lauching sites which are established on the perimeter missile lauching sites are considered, the dollars saved may be 
lants. of strategic areas and major cities could be severely handi- insignificant—but where, around strategic areas with high 
rmits capped because of area required by the use of Class LX pro- population densities, can 1200 acres of open land be found? 
ty to pellants. According to the existing ground rules, the weight In the foreseeable future, missile systems containing 250,000 
of the propellant mass and explosive in the warhead would be lb of propellant might be used. Such a mass, if Class LX 
d the combined and considered as Class IX. The proposed re- propellant, requires a minimum land area of 1340 acres. To 
Chief vision to the Ordnance Safety Manual would relieve this have only one dozen of these missiles available for ready use, 
lorate situation—normally, the distance between the high explosive 16,000 acres are needed. These areas can be reduced by a 
ite or warhead and the propellant is 6 in. or greater, and, with Class factor of four if proper barricades are provided. With Class 
with {I propellant, quantity-distance would be assessed only on the II propellants, the maximum areas required are 15 acres and 
Ord- weight of Class [X material in the warhead. 156 acres, respectively. 
ovide 
com- 
The i 
ratio 
New Patents 
t fuel 
r less _ aa George F. McLaughlin, Contributor 
n the Temperature compensated regulator for Baller and C. A. Johnstone, Indianapolis, Gyroscope with direct current torquing 
shall fluid supply lines (2,829,492). H. Klein- Ind., assignor to U. S. Navy. (2,824,451). T.O.Summers Jr., Sherman 
man, New York, N. Y., assignor to Reac- Means for intercepting the movement of Oaks, Calif., assignor to Summers Gyro- 
tion Motors, Ine. the missile launching tube to its vertical scope Co. 
? pro- Thermostatic control of a pressure bal- loading position so the tube and ejecting Cylindrical armature constructed of 
stric- anced valve in the supply line to a rocket means may be held in inclined alignment, non-magnetic material so as to have no 
-may engine. A threaded valve stem, rotated permitting the ejector to be extended into inherent. residual magnetism remaining 
ative in response to changes in ambient temper- the tube to eject the dud. after energization of the armature by 
from phn ne goed variation in the valve cross- — Primary zone for gas turbine combustor _—gravity sensitive means. 
secuional area. (2,829,494). W. L. Christensen, New- Exercise head for torpedoes or other 


A typical example of a missile lauching site is given in Fig. 


town, Conn., assignor to the U.S. Navy. 

Ceramic target mounted in the mixing 
chamber so hot gases from the primary 
combustion zone will pass between the 
target and a shield in the combustor. 
Incoming fuel is vaporized, and air passing 
through a venturi draws the mixture into 
the combustion zone. 


underwater, surface, or air missiles, ships, 
or craft (2,824,537). T. C. Boyle, Pasa- 
dena, Calif., assignor to the U. S. Navy. 
A movable wall adapted to reduce the 
volume of a chamber containing liquid 
ballast. A discharge valve opens auto- 
matically in response to pressure applied 
to the ballast by movement of the wall. 


Aerodynamic bomb ejector mechanism 
(2,829,561). E. H. Granfelt, Huntington 
Valley, Pa., assignor to the U. 8. Navy. 

Arms connected to an actuator flap 

under the wing of an aircraft arranged to 
exert a substantial downward force on the 
rocket or bomb immediately after it is re- 
leased. 
Apparatus for regulating fuel flow in gas 
turbine engines (2,829,662). F.G. Carey, 
Cheltenham, England, assignor to Dowty 
Fuel Systems Ltd. 

Power operated means moving a differ- 
ential member in the manual control ac- 
cording to difference between a selected 
speed signal and the engine speed signal. 
Apparatus for machining radome walls to 
controlled electrical thickness (2,824,413). 
A. J. Whitehill, Seattle, Wash., assignor to 
Boeing Airplane Co. 

A positional pick-up controls an actu- 
ator to move a machine tool toward the 
radome wall in response to predetermined 
displacement of a feeler element on a tra- 
versing head. 


2.831.921 


Variable —_ exhaust nozzle (2,831,321). 
G. Laucher, Canoga Park, ( valif., 
assignor to Marquardt Aircraft Co. 
In a turbojet engine, a circumferential 
force exerted against segments of a multi- 
segment nozzle varies the exit area upon 


7 Dud jettisoning device for rocket launcher 
: (2,831,399). M. H. Meekins, M. H. 


Eprrors Note: Patents listed above were selected from the Official Gazette of the U.S. 


- Patent Office. Printed copies of patents may be obtained from the Commissioner of fore-and-aft movement between the nozzle 
Patents, Washington 25, D. C., at a cost of 25 cents each; design patents, 10 cents. and a guide encircling the segments. 
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Book Reviews 


Transport Properties in Gases, Proceed- 
ings of the Second Biennial Gas Dy- 
namics Symposium at Northwestern 
University, edited by A. B. Cambel and 
John B. Fenn, Northwestern University 
Press, Evanston, 1958, 182 pp. $7.50 


Reviewed by A. K. OppENHEIM 
University of California 


It has become a custom over the recent 
years to hold special meetings devoted to 
some special branches of scientific activity 
and to publish their proceedings in sepa- 
rate volumes. A valuable addition to the 
literature resulting from such an activity 
last year is the present volume of the 
Proceedings of the Second Biennial Gas 
Dynamics Symposium co-sponsored by 
the AmerIcAN Rocket and 
Northwestern University. 

The importance of the knowledge of 
transport properties in gases to modern 
developments in gas dynamics hardly 
needs emphasis. Great strides made over 
the last fifty years in the understanding of 
the macroscopic mechanics of continuum 
are impressive. In equations formulating 
this approach, transport properties ap- 
pear as coefficients. Their magnitude and 
dependence on other parameters of the 
problem are often of critical importance, 
but their true character cannot be revealed 
by a macroscopic study, characteristic of 


Pressure Vessel 
Technology 
tanium... 


and now 


Ali Bulent Cambel, Northwestern University, Associate Editor 


the gas dynamic method of analysis. 
Instead microscopic inquiries based on 
other disciplines usually have to be made. 
However, the role transport parameters 
play in modern fluid dynamics is so pro- 
found that their study cannot be neglected 
by a gas dynamicist. 

The volume starts with an excellent 
review of the current status of the theories 
of gas transport properties presented by 
the distinguished trio of the University of 
Wisconsin: Curtis, Hirschfelder and 
Bird. After Touloukian’s description of 
the collection and correlation of data on 
thermophysical properties made currently 
at Purdue, there follows a section on 
experimental methods for the determina- 
tion of transport properties. 

The subject is opened by J. Kestin in a 
critical review of the problem of the tech- 
niques for the measurement of the vis- 
cosity of gases at elevated pressures and 
temperatures. From a general consider- 
ation of criteria for a successful instrument 
it appears that at present the best choice 
is offered by the capillary viscometer. 
Other alternatives, such as the oscillating 
body (on the development of which the 
author spent considerable effort), the 
rotating cylinder and the falling body 
viscometer, seem still to require an inten- 
sive period of development even before 
their significance could be accurately 
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BARBER 
COLMAN 


Controle 


Spheres of titanium alloy represent 
the optimum configuration and 


TEMPERATURE 


material of maximum strength/weight & POSITIONING 
for the storage of gases and liquids CONTROLS 
under extremes of pressure and 
temperature. 

Rheem Aircraft, presently supplying 
pressure units to current ICBM pro- as 


grams, is directing its years of experi- 
ence to the development and pro- 
duction of all types of vessels and 
propellant chambers. 


AIR VALVES 


ACTUATORS 
RHEEM MANUFACTURING CO. 


AIRCRAFT DIVISION GROUND TEST 
1 woodruff avenue, downey, EQUIPMENT 
Write for complete data 
a or consult the Barber- 


Colman engineering sales 
office nearest you: 

Los Angeles, Seattle, 
Fort Worth, New York, 
Baltimore, Montreal, 
Rockford. 


BARBER-COLMAN COMPANY 
Dept.G— 1470 Rock Street, Rockford, Illinois. 


Engineers: Join Rheem in challenging technical work. 


ght vehicles 


AERONUTRONIC SYSTEMS, INC., 2 subsidiary 
of Ford Motor Company, is undertaking expanded 
military and commercial programs involving the 
most advanced research, development, experimen- 
tation and prototype production at plants in Glen- 
dale and Van Nuys, California, and at modern, new 
facilities overlooking the Pacific Ocean at Newport 
Beach, California. The following positions are 
open: 


SENIOR STRUCTURAL ENGINEERS, with 
graduate degrees for design and analysis with ap- 
proximately 8 to 10 years’ experience and 3 to 4 
years’ supervisory experience in the missile field. 
Will be required to apply knowledge of high 
temperature materials and methods, thermal stress, 
dynamics, etc. to advanced hypersonic vehicles and 
re-entry bodies. 


PhD AND MS RESEARCH SPECIALISTS with 5 
to 7 years’ experience in heat transfer and fluid 
mechanics, thermodynamics, combustion and chemi- 
cal kinetics, thermoelasticity or metallurgy for as- 
signment on theoretical and experimental programs 
related to re-entry technology and advanced rocket 
propulsion. Specific assignments are open in re- 
entry body design, high temperature materials 
studies, boundary layer heat transfer with chemical 
reactions, thermal stress analysis, high temperature 
thermodynamics, studies of rocket chamber cooling 
methods, and liquid and solid rocket engine design. 
For assignment at our new Newport Beach facility. 


PHYSICAL CHEMISTS OR PHYSICISTS, Senior 
level. Advanced degree required with at least 5 
years’ experience in the analysis of combustion 
processes in high energy fuel systems. For assign- 
ment to work on combustion program of wide scope 
located in our new laboratory at Newport Beach, 
California. 


Qualified scientists and engineers are invited to 
contact Mr. L. T. Williams, Aeronutronic Systems, 
Inc., Building 19, 1234 Air Way, Glendale, 
California. 
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Why wait for the answers? © 


... Observe and evaluate phenomena as they 
occur—with new Kodak Linagraph Direct Print Paper! — 

This new photorecording paper gives immediate readout in 
moving-mirror galvanometer oscillographs—with no photochemical 
developing. Gives sharp, legible traces over a wide range of pian 4 
recording frequencies—from 0 cps to 3000 cps. 

Linagraph Direct Print records can be used in recording 
system readers, accept pen and pencil notations readily. They 
have adequate permanence and may be stabilized 
for extended archival storage. 

Now available in 5” x 200’, 6” x 100’, 7” x 200’, and 12” x 200’ rolls. 
Other sizes on request. For complete details, write: 


EASTMAN KODAK COMPANY 
Graphic Reproduction Division 
Rochester 4, N. Y. 
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we have 
the “Know-How!” 


We have developed new tech- 
niques, new methods, new 
processes that effect production 
economy so necessary to a success- 
ful missile program. 


Here at Newbrook you will find 
men with experience gained 


pecializing in 


@ Motor Cases from doing ...a modern plant 
Solid and Liquid % with up-to-date equipment... 
propellants Pe, precision inspection to meet your 

® Jato Cases . most exacting quality control 

requirements. 
Nozzles 


And most important, Newbrook 
specialization results in strict 
reliability ! Let us help you with 
om Missile Hardware problems. 


Plenum Chambers 
Blast Tubes 


Fuel Injectors 


NEWBROOK MACHINE CORPORATION. 
es 20 Mechanic Street Phone 45 
ae SILVER CREEK, NEW YORK 


assessed. The main difficulties encoun- 
tered in their use, the author points out, 
are connected with complexities of the 
flow field which they create. 

A résumé of recent measurements of 
thermal conductivities, especially in rela- 
tion to viscosity, is offered by Keyes, 
while Guildner describes recent improve- 
ments made at MIT on the use of the 
coaxial cylinder method for the measure- 
ment of thermal ‘conductivity in gases. 
As a demonstration of its virtue, he pre- 
sents data on thermal conductivity for 
carbon dioxide in the region of the critica! 
point and within a pressure range of from 
1 to 100 atm. 

Experimental methods for the measure- 
ment of diffusion coefficients in gases and 
liquids at elevated pressures are reviewed 
by Reamer and Sege. They then describe 
the equipment used by them at the Cali- 
fornia Institute of Technology which 
operates at pressures up to 10,000 lb per 
square inch in the temperature range from 
40 to 460 F. Typical results for some 
hydrocarbon liquids are presented. 

The fundamental aspects of thermal dif- 
fusion in gases are discussed by Chap- 
man. Since, as he points out, thermal 
diffusion in binary gas mixtures depends 
more than most gas phenomena on the 
precise nature of the forces between the 
unlike molecules, it presents a useful tool 
for the investigation of such forces. 

Basic features of the relaxation times in 
gases are elegantly summarized by Bou- 
dart. This is followed by Estermann’s 
description of molecular beam applications 
to the measurement of transport proper- 
ties in gases at low densities, Lin’s review 
of shock tube applications for the study of 
transport properties in gases at high 
temperatures, and Baulknight’s descrip- 
tion of the calculation of transport proper- 
ties at elevated temperatures. 

The determination of equilibrium gas 
emissivities from spectroscopic data is 
given a comprehensive treatment by Pen- 
ner and Thomson. They set down the 
physical principles involved in the evalu- 
ation of gas emissivities, that is, the black 
body radiation laws, the basic laws for 
distributed radiation, the Einstein coeffi- 
cients for induced absorption and emis- 
sivity and for spontaneous emissions, the 
absorption oscillator strengths and the 
spectral line profiles including the various 
line-broadening phenomena, especially 
those due to the Doppler and to the col- 
lison effects. They then consider the 
spectral distribution of the centers of 
vibration-rotation bands. 

The importance of the knowledge of 
transport properties in various specific 
applications is explored in a number of 
articles forming the third part of the book. 
Bartz considers it from the viewpoint of 
jet engines and rockets; Bromberg, with 
reference to hypersonic aerodynamics; 
Friend, with respect to chemical and 
process industries; Sherman, in low den- 
sity gases; and Rosenbluth, in relation to 
magnetohydrodynamics. The volume 
closes with a reference list compiled by 
Hilsenrath. 

There is no doubt in the reviewer’s mind 
that the volume forms a distinct stepping 
stone in our comprehension of the intri- 
‘acies of transport properties and their 
evaluation. 
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No other missile test 
lab can do all the | 
things we can do 


_ Take our 100-foot test tower. With it we can 
simulate the pressure and aerodynamic loadings 
that occur during the actual launching and flight 
of a missile. 

This is just one of the many specialized facil- 
ities and tools that we have built up in the last 
ten years. During this decade—because of our 
complete responsibility for the development of 
a major weapon system — we have also piled up 
a wealth of experience in many related fields. 

The upshot of this is that we can handle any 
project from far-out basic research to static and 
dynamic testing of complete airframes. 

For example, we can simulate aerodynamic 
heating and loading at the same time...on a 
complete wing or a large section of fuselage. 
We’re old hands at applying heat gradients from 
leading to trailing edge and spanwise. And we 
are equipped to compute temperature differen- 
tials, and to predict specific area temperatures. 

We’re completely qualified to work with any 
kind of metals or nonmetals; materials and 
processes; electrical, hydraulic and pneumatic 
components and systems; structures; reliability; 
environmental testing (including re-entry) ; and 
instrumentation. 

Now we'd like to solve your problems. The 
tougher the better. 

Please write to E. R. Schenkel, Missile Divi- 
sion Test Lab Manager, North American Avia- 
tion, 12214 Lakewood Blvd., Downey, Calif. 


MISSILE DEVELOPMENT DIVISION 


North American Aviation, Inc. 
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Jet and Rocket 
Engines 


Experimental Investigation of an Axial- 
flow Compressor Inlet State Operating at 
Transonic Relative Inlet Mach numbers, 
III: Blade-row Performance of Stage 
with Transonic Rotor and Subsonic Stator 
at Corrected Tip Speeds of 800 and 1000 
Feet Per Second, by Francis C. Schwenk, 
Seymour Lieblein and George W. Lewis 
Jr.. NACA RM E53G17, Sept. 1953, 69 
pp. (Declassified from Confidential by 
authority of NACA Res. Abs. 121, p. 17, 
11/4/57.) 

Diffusion Factor for Estimating Losses 
and Limiting Blade Loadings in Axial- 
flow-compressor Blade Elements, by 
Seymour Lieblein, Francis C. Schwenk and 
Robert L. Broderick, NACA RM E53D01, 
June 1953, 43 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abs. 
121, p. 17, 11/4/57.) 

Investigation of an Axial-flow-compres- 
sor Rotor with Circular-arc Blades Oper- 
ating up to a Rotor-inlet Relative Mach 
Number of 1.22, by William H. Robbins 
and Frederick W. Glaser, NACA RM 
E53D24, July 1953, 63 pp. (Declassified 


Ep1Tor’s Note: Contributions from E. R. 
G. Eckert, J. P. Hartnett, T. F. Irvine Jr. and 
P. J. Schneieer of the Heat Transfer Labora- 
tory, University of Minnesota, are gratefully 
acknowledged. 
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Self-Aligning Bearings 
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All Word Rights Reserved 


CHARACTERISTICS 


ANALYSIS 


Stainless Steel 
Ball and Race 


Chrome Alloy ff types operating under a radial 
ultimate loads (3000-893,000 Ibs.). 


Steel Ball and Race 


Bronze Race and 
Chrome Steel Ball 


Thousands in use. Backed by years of service life. Wide variety 
of Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in 
similar size range with externally or internally threaded shanks. 
Our Engineers welcome an opportunity of studying individual 
requirements and prescribing a type or types which will serve 
under your demanding conditions. Southwest can design special 
types to fit individual specifications. As a result of thorough 
study of different operating conditions, various steel alloys 
have been used to meet specific needs. Write for revised Engi- 
neering Manual describing complete line. Dept. JP-58. 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


from Confidential by authority of NACA 
Res. Abs. 121, p. 17, 11/4/57.) 


General Considerations of Mach Num- 
ber Effects on Compressor-blade Design, 
by John F. Klapproth, NACA RM E53- 
L23a, Apr. 1954, 24pp. (Declassified from 
Confidential by authority of NACA Res. 
Abs. 121, p. 18, 11/4/57.) 


Investigation of am Axial-flow Compres- 
sor Rotor Having NACA High-Speed 
Blade Sections (A:Isb Series) at Mean 
Radius Relative Inlet Mach Numbers up 
to 1.13, by Melvyn Savage, John R. Er- 


win and Robert P. Whitley, NACA RM 
L53G02, Nov. 1953, 43 pp. (Declassified 


from Confidential by authority of NACA 
Res. Abs. 121, p. 21, 11/4/57.) 


Low-speed Cascade Tests of Two 45° 
Swept Compressor Blades with Constant 
Spanwise Loading, by Loren A. Beatty, 
Melvyn Savage and James C. Emery, 
NACA RM 153107, March 1954, 65 pp. 
(Declassified from Confidential by ’author- 
ity of NACA Res. Abs. 121, p. 23, 11/4/ 


57.) 


Experimental Investigation of an Axial- 
flow Supersonic Compressor Having Sharp 
Leading-edge Blades with an 8-percent 
Mean Thickness-chord Ratio and of the 
Effect of a Leading-edge Radius, by 
Theodore J. Goldberg, NACA RM L54- 
K16, Feb. 1955, 33 pp. (Declassified 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


from Confidential by authority of NACA 
Res. Abs. 121, p. 25, 11/4/57.) 


Investigation of Turbines for Driving 
Supersonic Compressors, I: Design and 
Performance of First Configuration, by 
Warner L. Stewart, Harold J. Schum and 
Warren J. Whitney, NACA RM E52 
C25, June 1952, 27 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abs. 121, p. 16, 11/4/57.) © 


Aerodynamics of Jet 
Propelled Vehicles 


Experimental Aerodynamic Forces and 
Moments at Low Speed of a Missile Model 
During Simulated Launching from the 
Midsemispan Location of a 45° Sweptback 
Wing-fuselage Combination, by Williain 
J. Alford Jr., H. Norman Silvers and 
Thomas J. King Jr., NACA RM L54K11:, 
Feb. 1955, 36 pp. (Declassified froin 
Confidential by authority of NACA Res. 
Abs. 122, p. 138, 12/3/57.) 


Effect of Large Deflections of a Canard 
Control and Deflections of a Wing-tip 
Control on the Static-stability and In- 
duced-roll Characteristics of a Crucifor 
Canard Missile at a Mach Number of 
2.01, by M. Leroy Spearman, NACA RM 
L53K03, Dec. 1953, 20 pp. (Declassi- 
fied from Confidential by authority of 
NACA Res. Abs. 122, p. 12, 12/3/57.) 


ROD END 
TYPES 


Several years’ 


ENGINEER, ME 


FLUTTER & VIBRATION 


ANALYST 


SENIOR POSITION WITH 


GENERAL ELECTRIC 


experience in 


RECOMMENDED USE 


For types operating under high temper- 
ature (800-1200 degrees F.). 


For types operating under normal loads 
with minimum friction requirements. 


vibration analysis of turbo ma- 
chinery assemblies — plus a 
strong background in funda- 
mental concepts and applica- 
tions in this field—are essential 
for this advanced assignment. 
M.S. desirable. 

The work concerns the de- 
velopment of nuclear power 
systems for aircraft at General 
Electric. It deals specifically 
with high-temperature reactor 
core components. 

Nuclear experience is not es- 
sential. In-plant courses and 
graduate study—on full tuition 
refund basis—provide essential 
nuclear theory and technology. 


Please write in confidence, stating 
salary requirements to: 


Mr. P. W. Christos, Div. 34-MS 
AIRCRAFT NUCLEAR PROPULSION DEPT. 


GENERAL @@ ELECTRIC 


P. O. Box 182 Cincinnati 15, Ohio 
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Wing-Tunnel Investigations at Low and 
ape Transonic Speeds of the Feasibility of | 


Self-actuating Spoilers as a Lateral- 
control Device + a Missile, by Harleth 
ontribut G. Wiley and William C. Hayes Jr. 
ontributor RM L53K27, Jan. 1954, 24 pp. 
sity diagrs., tab. (Declassified from Confi- 
NAC dential by authority of NACA Res. Abs. 
of NACA 12, 12/3/57.) 
> Dee Molecular Approach to Problems of 
esien ari} High Altitude, High Speed Flight, by 
g G. N. Patterson, Univ. Toronto, Inst. 


Acrophys., Rev. 10, May 1957, 42 pp., 21 


RM 
ole Drag Measurements in Flight on the 10 


of NACA Percent Thick and 8 Percent Thick Wing 
: X-1 Airplanes, by John J. Gardner, NACA 
kM L8K05, Nov. 1948, 17 pp. (De- 
classified from Confidential by authority 
; 0! NACA Res. Abs. 122, p. 10, 12/3/57.) 
A Study of Several Parameters Con- 
trolling the Trajectories of a Supersonic 
wces and Antiaircraft Missile Powered with Solid 
le Model or Liquid Fuel Rockets, by Ralph F. 
rem the Huntsberger, NACA RM _ A6G22, Apr. 
westback 1947, 39 pp. (Declassified by authority 
Willian § °! NACA Res. Abs. 122, p. 1, 12/3/57.) 
vers and Theory of Bluff Bodies, by Isao Iami, 
54K1 1: Maryland Univ., Inst. Fluid Dynamics 
ed froin Appl. Math., TN BN-104 (ASTIA AD 
CA Res. 128412), June 1957, 89 pp. 


The Effects of Liquefaction of Air on the 
1 Canard Pressure Distribution and Forces on a 
Wing-tip Cone Cylinder Body at M-5.0, by W. T. 
and In- Strike and W. N. MacDermott, Arnold 
Crucifor Engng. Dev. Center, TN 57-14 (ASTIA 
mber of 1 D 123509), June 1957, 44 pp. 
CA RM Measurements of the Upwash on a 
Declassi- Cone Cylinder Model, by Oscar C. Hold- 
ority of erer, Army Ballistic Missile Agency, Dev. 
3/57.) Operations Div., Rep. DA-R-11, June 
1957, 23 pp. 

The Aerodynamics of Bodies in Non- 


CESSNA T-37A, 


dade 


JET MENTOR 


uniform Flow, by Enrico Pistolesi and r 
Marino Marini, Univ. Pisa, Inst. Aeron., 4 
(AFOSR TR-56-37), 1956, 38 pp., 11 
fig. @ 
Supersonic Flow Past Slender Bodies : $ 
of Arbitrary Shape, (Practical Numerical 
Methods for Electronic Computers), by ls 
Robert Sauer, Munchen, Technische Hoch- : [3 
schule. TN 1 (AFOSR-TN-57-162; AS- 
TIA AD 126454), 1957, 15 pp. Zz 8 
Zero-lift Drag of Several Conical and adi ma - @ 
Blunt Nose Shapes Obtained in Free 4 
Flight at Mach Numbers of 0.7 to 1.3, 7 4 
by Robert O. Piland and Leonard W. 
Putland, NACA RM 154427, March Asthe pioneerAmerican builderofsmallgas ¢ 
1954, 14 pp. | (Declassified from Confi-' — turbines, Continental Aviation and Engi- ® 
dential by authority of NACA Res. Abs. a 
121, p. 23, 11/4/57.) neering Corporation occupies a unique . 
Transonic Drag Measurements of Eight ona oon in this relatively new field of power. §$ 
Body-Nose Shapes, by William E. Stoney ot only is it the largest producer of engines ° 
RM 14 pp. of this type and size; it is also the most 
classifiec rom onfidential Dy au- 2 
thority of NACA. Res. Abs. 121, p. 22, thoroughly qualified by experience, andthe $ 
11/4/57.) most completely equipped. Those having H 
present or prospective need for gasturbine @ 
= 
Heat Transfer and Fluid power—for applications on land, at sea, or § 
Siew in the air—are invited to lay their problems ° 
before C. A. E. engineers. 
Experimental Study of Stagnation Tem- 
7 perature in Free Molecular Flow, by F. M. = $ 
7 Devienne, J. Aeron. Sci., vol. 24, June $ MODEL 141 
: 1957, pp. 403-406, 412. @ AIR COMPRESSOR 
Turbulent Transport of Heat and Mo- 
mentum from an Infinite Rough Plane, 


by T. H. Ellison, J. Fluid Mech., vol. 2, 


CONTINENTAL AVIATION & ENGINEERING CORPORATION 


; adiabatic Rotational Free Stream, by 
_ Franklin D. Hains, J. Aeron. Sci., mee | 24, 12700 KERCHEVAL AVENUE, DETROIT 15, MICHIGAN 
—" Aug. 1957, pp. 623-624 


, pp. 62: 
Note on the Boundary-layer Instability 
on a Flat Plate Stopped Suddenly, by SUBSIDIARY OF CONTINENTAL MOTORS CORPORATION 
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Challenge in Electronics 


Electronics activities are broad and fast-growing at 
_ Chance Vought. Projects involve advanced guidance and 


control and fire control systems for missiles and high- 
performance manned aircraft. They begin with investiga- 
Wi tions and theory and progress through systemization and 


packaging to detailed hardware design. Key responsibili- 
ae ties await additional men who are qualified in these areas. 
Advanced degrees are preferred. 
~» 
Stability and Control Engineer. E.E., M.E., or A.E. with emphasis 

on flight stability and control problems or dynamics. (Special 
consideration given graduate study or extensive experience in transients 
or closed loop stability analysis.) To assist in design of autopilot and 
control systems for high-performance missiles and aircraft. 


Antenna Design Engineer. E.E. or Physics Degree with demonstrated 
aptitude for antenna design. To join active projects involving design of 
flush-mounted, recessed and external antennas at all frequencies for 
very high-performance aircraft and missiles. 


Fire Control and Microwave Systems Engineer. Requires E.E. or Physics 
Degree; at least 2 years experience in radar, data link, or fire 
control systems; and strong ability in this work. 


Test Equipment Engineer. Requires E.E. or Physics Degree and at least 
2 years experience in this or related field. (Desirable: broad background in 
electronics design with emphasis on digital computers or microwave 
systems.) To join in the design of complete checkout systems for 
missiles and associated subsystems. 

Guidance Design Engineer. E. E. or Physics Degree, plus 2 or more years 
experience. To design various active and self-contained missile guidance 
systems, and to design and develop radar beacons. 

Reliability Analyst. Requires M.E., Physics, E.E., or Math Degree; 
broad knowledge of electronic and mechanical systems; experience in 
operations research or reliability. Helpful: statistical methods experience. 

Electronic Packaging Engineer. M.E., E.E., or equivalent packaging design 
experience. To help design ground, airborne and shipboard 
electronic equipment for use in severe environments. Involves 
consideration of heat transfer, shock, vibration and other factors. 

Field Service Engineer. B.S. in E.E. or M.E. and experience in airborne 
electronics equipment. To assist and consult with Navy personnel 
on operation and maintenance of autopilot and fire control systems 
for high performance aircraft. 


To arrange for a personal interview, or for more ——— a 
tion on these cr other current openings, return coupon to: 


C. A. Besio 
2 Supervisor, Engineering Personnel 
JP-1 


| 

- CHANCE VOUGHT AIRCRAFT, Dept. 


Dallas, Texas 


wn | Iama Engineer, 
interested in the opening for. 
| Name__— 
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Francis R. Hama, J. Aeron. Sci., vol. 24, 
June 1957, p. 471. 

Theory of the Interaction Between 
Oblique Shock Waves and Laminar 
Boundary Layers, by M. Honda, Univer- 
sity Tohoku, Rep. Inst. High Speed Mech., 
vol. 8, 1957, pp. 109-130. 

Dissociation Effects in Hypersonic Vis- 
cous Flow, by Y. H. Kuo, J. Aeron. Sci., 
vol. 24, May 1957, pp. 345-350. 

Charts Adapted from Van Driest’s 
Turbulent Flat-plate Theory for Deter- 
mining Values of Turbulent Aerodynamic 
Friction and Heat-Transfer Coefficients, 
by D. B. Lee and M. A. Faget, NACA TN 
3811, Oct. 1956, 16 pp. 

Average Properties of Compressible 
Laminar Boundary Layer on a Flat Plate 
with Unsteady Flight, by Franklin K. 
Moore and Simon Ostrach, NACA 7'N 
3886, Dec. 1956, 35 pp. 

Stability of Laminary Boundary Layer 
Near a Stagnation Point over an Imper- 
meable Wall and a Wall Cooled by Normal 
Fluid Injection, by Morris Morduchow, 
Richard G. Grape and Richard P. Shaw, 
NACA TN 4037, Aug. 1957, 56 pp. 

Simplified Method for Estimating Com- 
pressible Laminar Heat Transfer with 
Pressure Gradient, by Eli Reshotko. 
NACA TN 3888, Dec. 1956, 16 pp. 

Heat Transfer in a Laminar Boundary 
Layer at Mach 2.5 from a Surface Having a 
Temperature Distribution, by B. N. Prid- 
more-Brown, Univ. Toronto, Inst. Aero- 
phys., UTIA Rep. 45, Feb. 1957, 61 pp. 

The Effect of Fluid Injection on the 
Compressible Turbulent Boundary-layer, 
Preliminary Tests on Transpiration Cool- 
ing on a Flat Plate at M-2.7 with Air as the 
Injected Gas, by Morris W. Rubesin, 
Constantine C. Pappas and Arthur F. 
Okuno, NACA RM A55119, Dec. 1955. 

Heat Transfer to Laminar Boundary 
Layers with Variable Free-stream Veloc- 
ity, by R. A. Seban, Trans. ASME, 
vol. 79, no. 7, Oct 1957, p. 1545. 

Nonsteady Surface Temperature Ef- 
fects on Forced Convection Heat Transfer, 
by E. M. Sparrow and J. L. Gregg, J. 
Aeron. Sci., vol. 24, Oct. 1957, pp. 776- 

Summary of Low Prandtl Number Heat 
Transfer Results for Forced Convection 
on a Flat Plate, by E. M. Sparrow and 
J. L. Gregg, J. Aeron. Sci., vol. 24, Nov. 
1957, pp. 852-853. 


Combustion, Fuels and 
Propellants 


Volumetric and Phase Behavior in Part 
of the Nitric Acid-Nitrogen Dioxide-Water 
System, by H. H. Reamer and B. H. Sage, 
Calif. Inst. Techn., Jet. Propulsion Lab., 
Progr. Rep. 20-328, July 1957, 34 pp. 

Liquid Gas Containers Designed to 
Insure Liquefaction Temperature, Avia- 
tion Age, vol. 28, Dec. 1957, pp. 82-86. 

High-energy Fuels, by Donald R. 
Geehrig, Ordnance, vol. 42, Nov.—Deec. 
1957, pp. 563-566. 

The Hydrazine Flame, by M. Gilbert, 
Calif. Inst. Techn., Jet Propulsion Lab., 
Progress Rep. 20-318, May 1957, 13 pp. 

Gaseous Detonations X, Study of Reac- 
tion Zone, by John P. Chesnick and G. B. 
Kistiakowsky, Harvard University Dept. 
Chem., Oct. 1957, 14 pp. 

Combustion Instability in Liquid Pro- 
pellant Rocket Motors, Princeton Univ., 
Dept. of Aeron. Engng., 21st Quarterly 
Progr. Rep., 1 May-31 July 1957, Rep. 
216-u, Sept. 1957, 21 pp., 3 fig. 
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HE WAS too weak to speak... 

this child of eight. But the words 
were plain to see in her eyes: 

“Can’t you make me well 
again, Doctor?” 

It’s terribly hard ... even fora 
doctor who sees tragedy enough... 
to watch a child fade from the sun- 
light of life—a victim of cancer. 

We had succeeded in prolonging 
her life by many months—thanks 


a child die 
of cancer 


to recent advances in the treatment 
of leukemia. 


But that’s not enough! Cancer is 
a disease that ranks today as the 
Number 1 disease-killer of children. 
We can... we must... find ways 
to battle it, and win over it. 

Research, supported by the 
American Cancer Society, is striv- 
ing tirelessly towards that goal. 
Let’s give... boldly, generously 


to the American Cancer Society 
Crusade... and help eliminate this 
mortal enemy which will take the 
lives of more than 250,000 Ameri- 
cans this year alone. 


Send your gift to CANCER in 
care of your local post office. 


CANCER 
SOCIETY 
we 
j 
‘ 
JuLy 1958 


1,000 TO 100,000 MILES 


Cislunar Space 


300 TO 1,000 MILES 


Exosphere 
= 
2 
= 
2 
Mesosphere 
= 
Stratosphere 


EARTH TO 10 MILES 


Troposphere 


in confidence, to 


Send complete resume, 
MR. GEORGE R. HICKMAN 


Engineering Employment Manager 
MEPURBLIG AVIATION 


FARMINGDALE, LONG ISLAND, NEW YORK 


* SCIENTI 


Diversity of 
vital R&D 
PROJECTS AT 
EVERY ALTITUDE 
from Troposphere 
to Cislunar 
space with 
Republic Aviation 


THEORETICAL FLUID DYNAM- 


ICS—AE or Physics, MS or 
PhD. To conduct basic re- 
search in fluid dynamics related 
to hypervelocities of flight in 
rarefied atmosphere at orbital 
speeds. Requires strong aero- 
nautics, understanding kinetic 
theory of gases applied to field, 
Reynolds number and dissocia- 
tion effects. 


— DEVELOP- 


ENT—aAr. To perform para- 
metric studies in preliminary 
design stage of aircraft, missile 
or space ship projects. Needs 
imaginative cast of mind plus 
solid background in propulsion, 
aerodynamics, stability and con- 
trol, trajectory and wind tunnel 
testing. 


PROPULSION SYSTEMS—rc- 
quires strong theoretical back- 
ground in order to evaluate new 
propulsion systems for specific 
applications. Also responsible 
for general power plant and 
ejector performance analyses, 
formulation of engine control 
requirements. 


HEAT TRANSFER — AE or 
ME. Heavy experience in heat 
transfer, thermodynamics. To 
work in area of aerodynamic 
heating of aircraft or missiles 
and re-entry problems of space 
vehicles. 


ENGINE AIR-INLET & EX- 


HAUST—AE, ME. Experience 
in analysis of internal and ex- 
ternal aerodynamics at  super- 
sonic speeds. Responsible for 
coordinating inlet design with 
airframe and engine configura- 
tion. 


WEAPONS 
EVALUATION 


ea? 


- 


and 

MISSILE 

FEASIBILITY 
STUDIES 


ARMOUR RESEARCH FOUN- 
DATION offers the creative 
and imaginative individual 
the opportunity to work in 
small project groups on chal- 
lenging Operations Research 
and Weapons Systems pro- 
grams. These programs en- 
compass research on all forms 


of air and ground warfare. 


Knowledge of probability 
theory, statistics, and physics 
Experience and 


with 


required. 

familiarity modern 
weapons systems concepts are 
very desirable. Must be able 
to integrate results from many 
scientific disciplines into sys- 


tems studies. 


ARMOUR RESEARCH FOUN- 
DATION _ is 
metropolitan area offering 


located in a 


outstanding cultural and edu- 
cation advantages. For fur- 
ther information please write 


to: 


E. P. Bloch 
ARMOUR RESEARCH 
FOUNDATION 
of Illinois Institute 
of Technology 
10 West 35th Street, 
Chicago 16, Illinois. 
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COMPLETE 


TESTING 


FACILITIES 


* Qualification Tests 
* Evaluation Tests 

* Performance Tests 
* Environmental Tests 


Manufacturers of Metal Boss 
Seals to Military Specifications 


AIRCRAFT 
EQUIPMENT 


TESTING CO. 
1812 Fleet St., Baltimore 31, Md. 


ORleans 5-8337 ORleans 5-2222 


This is your chance to examine 
FREE the first complete chronicle of 
the men and events that launched 
the space age: What’s Going on in 
Space? by Commander David C. 
Holmes. 


The author, a rocket and missile ex- 
pert and member of the American 
Rocket Society, tells you the inside- 
Washington reaction to Sputnik .. . 
where the U.S. may still hold a lead 
in the rocket race... the a 
hour count down of E xplorer I . 
the X-15 project latest spac .- 
medicine discoveries . . . equipment 
and space-suit problems . . . up- 
coming manned satellites and moon 
shots... ahistory of rocket develop- 
ment since Goddard's earliest experi- 
ments . .. PLUS dozens of other 
vital topics. 
F or your FREE examination copy of 
“What's Going on in _ Space?” 
simply mail this ad to Funk & Wag- 
nalls, 153 East 24th St., New York 
10, N. Y. Then, after ten days you 
may either send $3.95 plus postage 
or return book without obligation. 
Enclose check for $3.95 now, and 
we'll pay postage. Refund and re- 
turn privilege guaranteed. 


Name 
Address 


1958 


INVITATION | 


| to the readers of | 
| JET PROPULSION | 
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Engineers — Scientists 


The Bell Aircraft Rockets Division forges ahead with new 
types of rocket engines and propellents to provide the 
higher thrust and greater efficiency needed to push mis- 
siles, satellites and manned space vehicles through the 
earth’s atmosphere into outer space. 

These programs are the outgrowth of over a decade of 
Bell experience in rocketry, beginning with the record- 
breaking and X-2 


and continuing with the development of rocket power 


supersonic rocket-powered aircraft, 


plants for projects like the Rascal air-to-surface missile. 


Continued growth and expansion in the Rockets Divi- 
sion have opened a number of select positions in the 


following fields: 


Structural Analysis 
Instrumentation 
Rotating Machinery 
Controls Development 


Propulsion Systems 
Development 


Combustion and Fuels 
Research 


Systems Installation Laboratory Testing 


To learn more about the personal opportunities and unex- 
celled benefits now available to you as a member of our 
Rockets Division engineering team, send resume of your 
qualifications to: Supervisor of Engineering Employment, 
Dept. K-32, Bell Aircraft Corporation, P.Q. Box One, 
Buffalo 5, New York. 
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pressure tight 
thin wall casting 


Here is a cast aluminum access door 
for an aircraft fuel cell. Of special 
interest is the fact that in its flat area 
the thin section is held to 0.080" + 
.010", yet it is pressure tight. Grooves 
for sealing gaskets are cast in to elimi- 
nate machining; the fuel cell access 
door is used as cast. 

The unusual foundry methods of 
Morris Bean & Co. readily lend them- 
selves to casting large areas of pres- 
sure tight thin walls, in compound 
curves or flat planes. Casting such sec- 
tions for your missile, aircraft, wave 
guide or electronics parts, may save 
you time and money. May we send 
you technical information? 
_ Morris Bean & Company 

Yellow Springs 6, Ohio 7 


the talented 


engineer C scientist 


APL OFFERS 
GREATER 
FREEDOM OF 

ACTIVITY 


APL has responsibility for 
the technical direction of 
much of the guided missile 
program of the Navy Bureau 
of Ordnance. As a result 
staff members participate in 
assignments of challenging 
scope that range from basic 
research to prototype test- 
ing of weapons and weapons 
systems. 


A high degree of freedom 
of action enables APL staff 
members to give free rein 
to their talents and ideas. 
Thus, professional advance- 
ment and opportunities to 
accept program responsi- 
bility come rapidly. Promo- 
tion is rapid, too, because 
of our policy of placing 
professional technical men 
at all levels of supervision. 


APL’s past accomplish- 
ments include: the first ram- 
jet engine, the Aerobee high 
altitude rocket, the super- 
sonic Terrier, Tartar, and 
Talos missiles. Presently the 
Laboratory is engaged in 
solving complex and ad- 
vanced problems leading to 
future weapons and weap- 
ons systems vital to the 
national security. Interested 
engineers and physicists are 
invited to address inquiries 
to: 


Professional Staff 


Appointments 


The Johns Hopkins University 
Applied Physics Laboratory 


8617 Georgia Avenue, Silver Spring, Md. 
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Donald W. Douglas, Jr. 
_ President, Douglas Aircraft Company a 
Lp m sure you’ve heard about Douglas projects like For you engineers who can help us move forward, 
Th Nike-Ajax, Nike-Hercules, Nike-Zeus, opportunities are almost as limitless as space 
Honest John, Genie and Sparrow. While these itself. 
are among the most important defense programs If you thrive on tough problems — and there 
in our nation today, future planning is moving are many — we'd like to discuss a future at 
into even more stimulating areas. Douglas with you. 
Working as we are on the problems of space Please write to Mr. C. C. LaVene 
flight and at the very borderline of the unknown, Douglas Aircraft Company, Box J-620 
engineering excellence in all fields is essential. Santa Monica, California 
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AMERICAN 
LATEX 


PRODUCTS CORPORATION 


of the new 


FREEDLANDER 


Research & 
Development 
Laboratory 


DEVOTED 
EXCLUSIVELY TO 
THE UNLIMITED 

APPLICATIONS OF 


STAFOAM 


STAFOAM. a unique type of 
plastic, offers an unequalled 
combination of properties for 
industrial design, such as 
foaming-in-place, adhesion, 
strength, thermal insulation, 
energy absorption, resistance 
to chemicals and solvents, 
water impermeability, flame 
resistance, and wide tempera- 
ture service. 

INQUIRE NOW! Your 
investigation may uncover an 
entirely new dimension in the 
practical use of foamed 
plastics. 

Revealing brochures and 
technical data on new miracle 
STAFOAM will be mailed 
upon request. 


ATEX 


PRODUCTS CORPORATION 


3341 West El Segundo Boulevard 
Hawthorne, California 
ORegon 8-5021 * OSborne 6-0141 


A DIVISION OF THE DAYTON RUBBER CO. 


atk Index to Advertisers 


AEROJET-GENERAL CORPORATION... 
D'Arcy Advertising Co., Los Angeles, Calif. 


AERONUTRONIC Systems, INc 
Honig, Cooper & Miner Adv., Los Angeles, Calif. 


AIRCRAFT EQuIPMENT TESTING COMPANY 
Mahool Advertising, Inc., Baltimore, Md. 


AMERICAN LATEX Propucts CoRPORATION 
Raymond E. Wallace Adv., Inc., Palos Verdes Estates, Calif. 


APPLIED Puysics LABORATORY, THE JOHNS Hopkins UNIVERSITY 
M. Belmont ver Standig, Inc., Washington, D. C. 


Armour FounpaTION OF ILLINOIS INSTITUTE OF TECHNOLOGY. . 


AsTRODYNE, INC 
Batten, Barton, Durstine & Osborn, Inc., Los Angeles, Calif. 


BarBER-CoLMAN CoMPANY 
Howard H. Monk & Associates, Inc., Rockford, Ill. 


BELL Arrcrart CoRPORATION 
Baldwin Bowers & Strachan, Inc., Buffalo, N. Y. 


CHANCE VouGut ArRcrAFT, INC 
Tracy-Locke Company, Inc., Dallas, Texas 


Coast MANuFactuRING & SuppLy CoMPANY 
Peter Hurst Adv., Inc., San Francisco, Calif. 


CoNTINENTAL AVIATION & ENGINEERING COMPANY 
The Hopkins Agency, Detroit, Mich. 


DIvERSEY ENGINEERING COMPANY 
Roark & Colby Advertising, Chicago, Ill. 


Dovuetas AtrcraFT Company, INc 
J. Walter Thompson Co., Los Angeles, Calif. 


EastMAN Kopak ComMpaNy 
The Rumrill Co., Rochester, N. Y. 


Funk & WaAGNALLS 
Albert Frank-Guenther Law, Inc., New York, N. ¥ 


GENERAL ELEctTRIC COMPANY 
AIRcRAFT NUCLEAR PROPULSION DEPT 
Deutsch & Shea, Inc., New York, N. Y. 


LUKENS STEEL CoMPANY 
J. M. Mathes, Inc., New York, N. 


Morris Bean & CoMPpAaNy 
Odiorne Industrial Adv., Inc., Yellow Springs, Ohio 


NEwBROOK MacuHiNE CorPORATION 
Melvin F. Hall Adv. Agency, Inc., Buffalo, N. 


NITROGEN Division oF ALLIED CHEMICAL CORPORATION 
G. M. Basford Co., New York, N. Y. 


Nortu AMERICAN Aviation, INC 
Batten, Barton, Durstine & Osborn, Inc., Los Angeles, Calif. 


OLDSMOBILE Division, GENERAL Motors 
D. P. Brother & Co., Detroit, Mich. 


The McCarty Co., Los Angeles, Calif. 


AVIATION CORPORATION 
Deutsch & Shea, Inc., New York, N. Y. 


RHEEM MANUFACTURING COMPANY 
Getz & Sandborg, Inc., Beverly Hills, Calif. 


SouTHWEsT Propucts CoMPANY 
O. K. Fagan Adv. Agency, Los Angeles, Calif. 


THIOKOL CHEMICAL CORPORATION. 
Dancer-Fitzgerald-Sample, Inc., New York, N. Y. 


U.S. Arr Force 
Erwin Wasey, Ruthrauff & Ryan, Inc., New York, N. 
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